النفاذية النسبية

في التدفق متعدد الأطوار في الأوساط المسامية ، تُعدّ النفاذية النسبية لطورٍ ما مقياسًا لا بُعديًا للنفاذية الفعّالة لهذا الطور. وهي نسبة النفاذية الفعّالة لهذا الطور إلى النفاذية المطلقة. ويمكن اعتبارها تكييفًا لقانون دارسي ليناسب التدفق متعدد الأطوار.

التركيبة

بالنسبة للتدفق ثنائي الطور في وسط مسامي في ظل ظروف الحالة المستقرة ، يمكننا كتابة

qأنا=-كأناμأناPأنالأنا=1،2{\displaystyle q_{i}=-{\frac {k_{i}}{\mu _{i}}}\nabla P_{i}\qquad {\text{for}}\quad i=1,2}

أينqأنا{\displaystyle q_{i}}هو التدفق ،Pأنا{\displaystyle \nabla P_{i}}هو انخفاض الضغط ،μأنا{\displaystyle \mu _{i}}هي اللزوجة . الرمز السفليأنا{\displaystyle i}يشير ذلك إلى أن المعلمات خاصة بالطورأنا{\displaystyle i}.

كأنا{\displaystyle k_{i}}هنا نفاذية الطور (أي النفاذية الفعالة للطور)أنا{\displaystyle i}), كما لوحظ من خلال المعادلة أعلاه.

النفاذية النسبية ،كرأنا{\displaystyle k_{\mathit {ri}}}، للمرحلةأنا{\displaystyle i}ثم يتم تعريفها منكأنا=كرأناك{\displaystyle k_{i}=k_{\mathit {ri}}k}، مثل

كرأنا=كأنا/ك{\displaystyle k_{\mathit {ri}}=k_{i}/k}

أينك{\displaystyle k}تمثل نفاذية الوسط المسامي في التدفق أحادي الطور، أي النفاذية المطلقة . يجب أن تكون النفاذية النسبية بين الصفر والواحد.

في التطبيقات، غالبًا ما يتم تمثيل النفاذية النسبية كدالة لتشبع الماء ؛ ومع ذلك، نظرًا للتخلف الشعري، غالبًا ما يلجأ المرء إلى دالة أو منحنى تم قياسه أثناء التصريف وآخر تم قياسه أثناء الامتصاص .

وفقًا لهذا النهج، يُعيق وجود الأطوار الأخرى تدفق كل طور. وبالتالي، يكون مجموع النفاذية النسبية لجميع الأطوار أقل من 1. مع ذلك، تم الحصول على نفاذية نسبية ظاهرية أكبر من 1 لأن نهج دارسي يتجاهل تأثيرات الاقتران اللزج الناتج عن انتقال الزخم بين الأطوار (انظر الافتراضات أدناه). قد يُعزز هذا الاقتران التدفق بدلًا من إعاقته. وقد لوحظ ذلك في مكامن النفط الثقيل عندما يتدفق الطور الغازي على شكل فقاعات أو بقع (غير متصلة). [ 1 ]

افتراضات النمذجة

يُطلق على الصيغة المذكورة أعلاه لقانون دارسي أحيانًا اسم قانون دارسي الموسع، وهو مُصاغ لتدفق متعدد الأطوار أفقي أحادي البعد وغير قابل للامتزاج في وسط مسامي متجانس ومتناحٍ . يُهمل هذا النموذج التفاعلات بين السوائل، لذا يفترض أن الوسط المسامي الصلب والسوائل الأخرى تُشكل مصفوفة مسامية جديدة يمكن أن يتدفق عبرها طور ما، مما يعني أن أسطح التماس بين السوائل تظل ثابتة في حالة التدفق المستقر، وهو أمر غير صحيح، ولكن هذا التقريب أثبت جدواه على أي حال.

يجب أن يكون كل تشبع من تشبعات الطور أكبر من التشبع غير القابل للاختزال، ويفترض أن كل طور مستمر داخل الوسط المسامي.

استنادًا إلى بيانات من تجارب مختبر تحليل اللب الخاص (SCAL)، [ 2 ] يمكن بناء نماذج مبسطة للنفاذية النسبية كدالة للتشبع (مثل تشبع الماء ). ستركز هذه المقالة على نظام النفط والماء.

مقياس التشبع

تشبع الماءSw{\displaystyle S_{\mathit {w}}}يمثل نسبة حجم المسام المملوءة بالماء، وينطبق الأمر نفسه على تشبع الزيت.So{\displaystyle S_{\mathit {o}}}وبالتالي، فإن التشبعات هي في حد ذاتها خصائص أو متغيرات مُقاسة. وهذا يُعطي القيد.

Sw+So=1Sw=1-So{\displaystyle S_{\mathit {w}}+S_{\mathit {o}}=1\Leftrightarrow S_{\mathit {w}}=1-S_{\mathit {o}}}

لذلك، عادةً ما تُكتب دوال النموذج أو العلاقات الخاصة بالنفاذية النسبية في نظام الزيت والماء كدوال لتشبع الماء فقط، وهذا ما يجعل من الطبيعي اختيار تشبع الماء كمحور أفقي في العروض البيانية.Swأنار{\displaystyle S_{\mathit {wir}}}(يشار إليه أيضًا بـSwج{\displaystyle S_{\mathit {wc}}}وأحياناًSwر{\displaystyle S_{\mathit {wr}}}ليكن ) تشبع الماء غير القابل للاختزال (أو الأدنى أو الملازم)، وليكنSoرw{\displaystyle S_{\mathit {orw}}}يمثل هذا الحد الأدنى من تشبع الزيت المتبقي بعد عملية حقن الماء (الامتصاص). وتكون نافذة تشبع الماء المتدفق في عملية غزو/حقن/امتصاص الماء محدودة بقيمة دنيا.Swأنار{\displaystyle S_{\mathit {wir}}}وقيمة قصوىSwoر=1-Soرw{\displaystyle S_{\mathit {wor}}=1-S_{\mathit {orw}}}رياضياً، تُكتب نافذة التشبع المتدفقة على النحو التالي:

SwأنارSwSwoر=1-Soرw{\displaystyle S_{\mathit {wir}}\leq S_{\mathit {w}}\leq S_{\mathit {wor}}=1-S_{\mathit {orw}}}
توحيد قيم تشبع الماء

من خلال تعديل تشبع الماء ليناسب نطاق التشبع المتدفق، نحصل على قيمة تشبع ماء مُعَيَّرة (جديدة أو أخرى).

Swن=Swن(Sw)=Sw-Swأنار1-Swأنار-Soرw=Sw-SwأنارSwoر-Swأنار{\displaystyle S_{\mathit {wn}}=S_{\mathit {wn}}(S_{w})={\frac {S_{w}-S_{\mathit {wir}}}{1-S_{\mathit {wir}}-S_{\mathit {orw}}}}={\frac {S_{w}-S_{\mathit {wir}}}{S_{\mathit {wor}}-S_{\mathit {wir}}}}}

وقيمة تشبع الزيت المعيارية

Soن=1-Swن=So-Soرw1-Swأنار-Soرw{\displaystyle S_{\mathit {on}}=1-S_{\mathit {wn}}={\frac {S_{o}-S_{\mathit {orw}}}{1-S_{\mathit {wir}}-S_{\mathit {orw}}}}}

نقاط النهاية

يترككرow{\displaystyle K_{\mathit {row}}}لتكن النفاذية النسبية للنفط، ولتكنكرw{\displaystyle K_{\mathit {rw}}}لنفترض أن النفاذية النسبية للماء هي نفاذية الماء. هناك طريقتان لقياس نفاذية الطور (أي النفاذية الفعالة للطور). إذا قمنا بقياس نفاذية الطور بالنسبة إلى نفاذية الماء المطلقة (أيSw=1{\displaystyle S_{\mathit {w}}=1}نحصل على مُعامل نقطة نهاية لكل من النفاذية النسبية للنفط والماء. إذا قمنا بقياس نفاذية الطور بالنسبة لنفاذية النفط مع وجود تشبع غير قابل للاختزال بالماء، كرow{\displaystyle K_{\mathit {row}}}نقطة النهاية هي واحدة، ولم يتبق لدينا سوىكرw{\displaystyle K_{\mathit {rw}}}معلمات نقطة النهاية. لتحقيق كلا الخيارين في النموذج الرياضي، من الشائع استخدام رمزين لنقطة النهاية في نموذج النفاذية النسبية ثنائية الطور. معلمات نقطة النهاية للنفاذية النسبية للنفط والماء هي:

كرow(Swأنار)=كرoتوكرw(Swoر)=كرwر{\displaystyle {\begin{aligned}K_{\mathit {row}}(S_{\mathit {wir}})=K_{\mathit {rot}}&&{\text{and}}&&K_{\mathit {rw}}(S_{\mathit {wor}})=K_{\mathit {rwr}}\end{aligned}}}

لهذه الرموز مزاياها وحدودها. الرمزكرoت{\displaystyle K_{\mathit {rot}}}التأكيد على أنها تمثل أعلى نقطة فيكرow{\displaystyle K_{\mathit {row}}}يحدث ذلك عند تشبع الماء غير القابل للاختزال، وهو أكبر قيمة لـكرow{\displaystyle K_{\mathit {row}}}قد يحدث ذلك عند التشبع الأولي بالماء. رمز نقطة النهاية المنافسةكرoر{\displaystyle K_{\mathit {ror}}}يحدث ذلك في تدفق الامتصاص في أنظمة النفط والغاز. إذا كانت أساس النفاذية هو النفط مع وجود ماء غير قابل للاختزال، فإنكرoت=1{\displaystyle K_{\mathit {rot}}=1}الرمزكرwر{\displaystyle K_{\mathit {rwr}}}يؤكد ذلك على حدوثه عند تشبع الزيت المتبقي. رمز بديل لـكرwر{\displaystyle K_{\mathit {rwr}}}يكونكرwo{\displaystyle K_{\mathit {rw}}^{o}}مما يؤكد أن النفاذية المرجعية هي نفاذية النفط مع الماء غير القابل للاختزالSwأنار{\displaystyle S_{\mathit {wir}}}حاضر.

ثم تُكتب نماذج النفاذية النسبية للزيت والماء على النحو التالي

كرow=كرoتكرowن(Swن)وكرw=كرwركرwن(Swن){\displaystyle {\begin{aligned}K_{\mathit {row}}=K_{\mathit {rot}}\cdot K_{\mathit {rown}}(S_{\mathit {wn}})&&{\text{and}}&&K_{\mathit {rw}}=K_{\mathit {rwr}}\cdot K_{\mathit {rwn}}(S_{\mathit {wn}})\end{aligned}}}

الوظائفكرowن{\displaystyle K_{\mathit {rown}}}وكرwن{\displaystyle K_{\mathit {rwn}}}تُسمى هذه المعاملات بالنفاذية النسبية المعيارية أو دوال الشكل للنفط والماء، على التوالي. معلمات نقطة النهايةكرoت{\displaystyle K_{\mathit {rot}}}وكرwر{\displaystyle K_{\mathit {rwr}}}(وهو تبسيط لـكرwoر{\displaystyle K_{\mathit {rwor}}}) هي خصائص فيزيائية يتم الحصول عليها إما قبل أو مع تحسين معلمات الشكل الموجودة في دوال الشكل.

غالبًا ما تحتوي المقالات التي تناقش نماذج النفاذية النسبية على العديد من الرموز. ويتجنب عدد من محللي اللب ومهندسي المكامن والعلماء المشغولين استخدام الرموز السفلية المملة والمستهلكة للوقت، ويكتبون على سبيل المثال Krow بدلاً من ذلك.كرow{\displaystyle K_{\mathit {row}}}أوكرow{\displaystyle k_{\mathit {row}}}أو النفاذية النسبية للنفط أو النفاذية النسبية للغاز الخامل. لذلك، من المتوقع وجود مجموعة متنوعة من الرموز، وهي مقبولة طالما تم شرحها أو تعريفها.

تمت مناقشة تأثيرات شروط الانزلاق أو عدم الانزلاق على حدود تدفق المسام على معلمات نقطة النهاية من قبل بيرج وآخرون [ 3 ] [ 4 ] .

نموذج كوري

يُعدّ معامل كوري [ 5 ] [ 6 ] [ 7 ] تقريبًا شائع الاستخدام للنفاذية النسبية، وهو عبارة عن قانون أُسّي في حالة التشبع. وتكون معاملات كوري للنفاذية النسبية للنفط والماء كما يلي:

مثال على ارتباط كوري لتدفق الامتصاص معشمالo{\displaystyle N_{\mathit {o}}}=شمالw=2{\displaystyle N_{\mathit {w}}=2}وكرwر=0.6{\displaystyle K_{\mathit {rwr}}=0.6}.
كرow(Sw)=كرoت(1-Swن)شمالo{\displaystyle K_{\mathit {row}}(S_{w})=K{_{\mathit {rot}}}(1-S_{\mathit {wn}})^{N_{\mathit {o}}}}
كرw(Sw)=كرwرSwنشمالw{\displaystyle K_{\mathit {rw}}(S_{w})=K{_{\mathit {rwr}}}S_{\mathit {wn}}^{N_{\mathit {w}}}}

إذا كانت أساس النفاذية هو النفط العادي مع وجود ماء غير قابل للاختزال، فإنكرoت=1{\displaystyle K_{\mathit {rot}}=1}.

المعايير التجريبيةشمالo{\displaystyle N_{\mathit {o}}}وشمالw{\displaystyle N_{\mathit {w}}}تُسمى هذه المعلمات معلمات شكل المنحنى أو ببساطة معلمات الشكل، ويمكن الحصول عليها من البيانات المقاسة إما عن طريق التفسير التحليلي للبيانات المقاسة، أو عن طريق التحسين باستخدام محاكاة عددية لتدفق النواة لمطابقة التجربة (والتي تسمى غالبًا مطابقة التاريخ).شمالo=شمالw=2{\displaystyle N_{\mathit {o}}=N_{\mathit {w}}=2}يكون ذلك مناسبًا في بعض الأحيان. الخصائص الفيزيائيةكرoت{\displaystyle K_{\mathit {rot}}}وكرwر{\displaystyle K_{\mathit {rwr}}}يتم الحصول عليها إما قبل أو بالتزامن مع تحسينشمالo{\displaystyle N_{\mathit {o}}}وشمالw{\displaystyle N_{\mathit {w}}}.

في حالة نظام الغاز والماء أو نظام الغاز والنفط، توجد علاقات كوري مشابهة لعلاقات النفاذية النسبية للنفط والماء الموضحة أعلاه.

نموذج LET

لا يمتلك نموذج كوري، أو معادلة كوري، سوى درجة حرية واحدة لشكل كل منحنى نفاذية نسبية، وهي معامل الشكل N. أما معادلة LET [ 8 ] [ 9 ] فتضيف درجات حرية إضافية لاستيعاب شكل منحنيات النفاذية النسبية في تجارب SCAL [ 2 ] وفي نماذج الخزانات ثلاثية الأبعاد التي يتم تعديلها لتتوافق مع الإنتاج التاريخي. وتشمل هذه التعديلات عادةً منحنيات النفاذية النسبية ونقاط النهاية.

مثال على ارتباط LET لتدفق الامتصاص مع L وE وT جميعها تساوي 2 وكرwر=0.6{\displaystyle K_{\mathit {rwr}}=0.6}.

يتم وصف تقريب LET بثلاثة معلمات L و E و T. وبالتالي، فإن العلاقة بين النفاذية النسبية للماء والنفط مع حقن الماء هي

كرw=كرwرSwنلwSwنلw+هـw(1-Swن)تيw{\displaystyle K_{\mathit {rw}}={\frac {{K_{\mathit {rwr}}}S_{\mathit {wn}}^{L_{\mathit {w}}}}{{S_{\mathit {wn}}}^{L_{\mathit {w}}}+{E_{\mathit {w}}}{(1-S_{\mathit {wn}})}^{T_{\mathit {w}}}}}}

و

كرow=كرoت(1-Swن)لo(1-Swن)لo+هـoSwنتيo{\displaystyle K_{\mathit {row}}={\frac {K_{\mathit {rot}}(1-S_{\mathit {wn}})^{L_{o}}}{{(1-S_{\mathit {wn}})^{L_{o}}}+{E_{\mathit {o}}}S_{\mathit {wn}}^{T_{\mathit {o}}}}}}

مكتوب باستخدام نفسSw{\displaystyle S_{w}}التطبيع كما هو الحال بالنسبة لكوري.

فقطSwأنار{\displaystyle S_{\mathit {wir}}}،Soرw{\displaystyle S_{\mathit {orw}}}،كرoت{\displaystyle K_{\mathit {rot}}}، وكرwر{\displaystyle K_{\mathit {rwr}}}لها معنى فيزيائي مباشر، بينما تُعتبر المعاملات L و E و T تجريبية. يصف المعامل L الجزء السفلي من المنحنى، وبالتشابه والخبرة، يمكن مقارنة قيم L بمعامل كوري المناسب. يصف المعامل T الجزء العلوي من المنحنى بطريقة مشابهة لوصف المعامل L للجزء السفلي. يصف المعامل E موضع ميل المنحنى (أو ارتفاعه). القيمة 1 هي قيمة محايدة، ويتحكم المعاملان L و T في موضع الميل . زيادة قيمة المعامل E تدفع الميل نحو الطرف العلوي من المنحنى، بينما تقليل قيمة المعامل E يدفع الميل نحو الطرف السفلي. تشير الخبرة في استخدام معادلة LET إلى النطاقات المعقولة التالية للمعاملات L و E و T : L ≥ 0.1، E > 0، و T ≥ 0.1.

في حالة نظام الغاز والماء أو نظام الغاز والنفط، توجد علاقات LET مماثلة لعلاقات النفاذية النسبية للنفط والماء الموضحة أعلاه.

التقييمات

بعد أن وضع موريس موسكات وآخرون مفهوم النفاذية النسبية في أواخر ثلاثينيات القرن العشرين، ازداد عدد نماذج النفاذية النسبية بشكل مطرد. وقد أدى ذلك إلى الحاجة لتقييم أكثر النماذج شيوعًا في الوقت الراهن. ومن أحدث وأشمل هذه التقييمات (حتى عام ٢٠١٩) ما قام به مقدسي وآخرون [ ١٠ ] وسخائي وآخرون [١١]. قيّم مقدسي وآخرون [ ١٠ ] نماذج كوري ، وشيريسي ، وLET للنفاذية النسبية للنفط/الماء باستخدام طريقة متطورة تأخذ في الحسبان عدد معلمات النموذج غير المؤكدة. ووجدوا أن نموذج LET، الذي يحتوي على أكبر عدد من المعلمات غير المؤكدة (ثلاثة)، هو الأفضل بوضوح لكل من النفاذية النسبية للنفط والماء. قام Sakhaei et alios [ 11 ] بتقييم 10 علاقات نفاذية نسبية شائعة الاستخدام لأنظمة الغاز/النفط والغاز/المكثفات، ووجدوا أن LET أظهر أفضل توافق مع القيم التجريبية لكل من النفاذية النسبية للغاز والنفط/المكثفات.

انظر أيضاً

مراجع

  1. برافو، إم سي؛ أراوجو، إم. (2008). "تحليل السلوك غير التقليدي للنفاذية النسبية للنفط أثناء اختبارات استنزاف الزيوت الثقيلة المشبعة بالغاز". المجلة الدولية لتدفق متعدد الأطوار . 34 (5): 447-460 . Bibcode : 2008IJMF...34..447B . doi : 10.1016/j.ijmultiphaseflow.2007.11.003 .
  2. 1 2 ماكفي، سي.؛ ريد، ج.؛ زوبيزاريتا، آي. (2015). التحليل الأساسي: دليل لأفضل الممارسات . إلسيفير. ISBN 978-0-444-63533-4.
  3. بيرغ، إس.؛ سينس، إيه دبليو؛ هوفمان، جيه بي؛ سميتس، آر إم إم (2007). "التدفق في الأوساط المسامية مع شرط حدودي انزلاقي". ورقة SCA2007-13 مقدمة في الندوة الدولية لعام 2007 لجمعية SCA، كالجاري، كندا، 10-12 سبتمبر 2007 .
  4. بيرغ، س.؛ سينس، أ. و.؛ هوفمان، ج. ب.؛ سميتس، ر. م. م. (2008). "تدفق ثنائي الطور في وسط مسامي مع شرط حدودي انزلاقي". النقل في الأوساط المسامية . 74 (3): 275-292 . Bibcode : 2008TPMed..74..275B . doi : 10.1007/s11242-007-9194-4 . S2CID 37627662 . 
  5. جودا، إتش إم؛ بيرنبروش، بي. (2004). استخدام نموذج بروكس-كوري المعدل لدراسة النفاذية النسبية للنفط والماء في هياكل مسامية متنوعة . doi : 10.2118/88538-MS . ISBN 978-1-55563-979-2.{{cite book}}تم |journal=تجاهله ( مساعدة )
  6. بروكس، آر إتش؛ كوري، إيه تي (1964). "الخواص الهيدروليكية للوسط المسامي". أوراق هيدرولوجية . 3 .
  7. كوري، أ. ت. (نوفمبر 1954). "العلاقة المتبادلة بين النفاذية النسبية للغاز والنفط". مجلة الإنتاج الشهرية . 19 (1): 38-41 .
  8. لوميلاند، ف.؛ إيبلتوفت، إ.؛ توماس، و.هـ. (2005). "علاقة جديدة متعددة الاستخدامات للنفاذية النسبية" (ملف PDF) . وقائع الندوة الدولية لعام 2005 لجمعية الكيمياء التطبيقية، أبو ظبي، الإمارات العربية المتحدة، 31 أكتوبر - 2 نوفمبر 2005 .
  9. لوميلاند، ف. (2018). "نظرة عامة على عائلة LET من الارتباطات متعددة الاستخدامات لوظائف التدفق" (ملف PDF) . وقائع الندوة الدولية لعام 2018 لجمعية SCA، تروندهايم، النرويج، 27-30 أغسطس 2018 .
  10. 1 2 مقدسي، ل.؛ غوادانيني، أ.؛ إنزولي، ف.؛ بارتوسيك، م. (2015). "تفسير منحنيات النفاذية النسبية ثنائية الطور من خلال صيغ متعددة ومعايير جودة النموذج". مجلة علوم وهندسة البترول . 135 : 738-749 . Bibcode : 2015JPSE..135..738M . doi : 10.1016/j.petrol.2015.10.027 . hdl : 11311/968828 .
  11. 1 2 سخائي، ز.؛ عازين، ر.؛ عوسفوري، س. (2016). "تقييم علاقات النفاذية النسبية التجريبية/النظرية لأنظمة الغاز والنفط/المكثفات". ورقة بحثية قُدِّمت في المؤتمر الأول لثنائية النفط والغاز والبتروكيماويات في الخليج العربي، الذي عُقد في جامعة الخليج العربي في بوشهر، إيران، في 20 أبريل 2016 .