Дарфис заңы көп фазалы ағынға - Википедия - Darcys law for multiphase flow

Моррис Мускат т.б.[1][2] дамыды басқарушы теңдеулер көпфазалы ағын үшін (бір вектор әрбір сұйықтық үшін теңдеу фаза ) кеуекті медиа Дарси теңдеуін қорыту ретінде (немесе Дарси заңы ) кеуекті ортадағы су ағыны үшін. Кеуекті орталар әдетте болады шөгінді жыныстар сияқты крастикалық жыныстар (негізінен құмтас ) немесе карбонатты жыныстар.

Номенклатура
ТаңбаСипаттамаSI бірліктері
индекс: а фазасы, векторлық компонент
векторлық компонент екі тор ұяшығының арасындағы бағытты жанасу бетініңм2
жанасу беті екі (әдетте көрші) тор ұяшықтары арасындам2
3-ось бойындағы бірлік вектор (z - бұл жерде еске салғыш: 3 - z бағыты)1
ауырлық күшінің үдеуіХаным2
ауырлық күшінің үдеуі бағытпенХаным2
абсолютті өткізгіштік 3x3 ретінде тензорм2
a = w, o, g фазасының салыстырмалы өткізгіштігібөлшек
бағытты салыстырмалы өткізгіштік (яғни 3х3 тензор)бөлшек
қысымПа
торлы жасушаның жанасу беті арқылы көлемдік ағын (Дарси жылдамдығы)Ханым
тор ұяшығының жанасу беті арқылы көлемдік ағыным3/ с
ағынның жылдамдығы (сұйықтық)Ханым
Дарси (сұйықтық) осі бойынша жылдамдық Ханым
Дарси (сұйықтық) жылдамдығыХаным
градиент операторм−1
динамикалық тұтқырлықПа с
масса тығыздықкг / м3
 Мұндағы a = w, o, g

Сұйықтықтың қазіргі фазалары су, мұнай және газ болып табылады және олар сәйкесінше a = w, o, g индексімен ұсынылған. Гравитациялық үдеу бағытымен бейнеленген немесе немесе . Мұнай инженериясында кеңістіктік координаттар жүйесі z осі төмен қаратып оңға бағытталғанына назар аударыңыз. Үш сұйықтық фазасының ағын теңдеулерін байланыстыратын физикалық қасиет мынада салыстырмалы өткізгіштік әрбір сұйықтық фазасы мен қысымы. Сұйық-жыныстық жүйенің бұл қасиеті (яғни су-мұнай-газ-тас жүйесі) негізінен сұйықтықтың функциясы болып табылады қанықтылық және ол байланысты капиллярлық қысым ағынды процесс, оған бағынышты екенін білдіреді гистерезис әсер.

1940 жылы б.з.б. Леверетт[3] қосу үшін деп атап көрсетті капиллярлық қысым ағын теңдеуіндегі әсер, қысым фазаға тәуелді болуы керек. Ағын теңдеуі содан кейін болады

 Мұндағы a = w, o, g

Леверетт сонымен қатар капиллярлық қысымның айтарлықтай екенін көрсетті гистерезис әсерлер. Бұл дренаж процесіне арналған капиллярлық қысымның an капиллярлық қысымнан өзгеше екенін білдіреді имбиция бірдей сұйықтық фазаларымен жүреді. Гистерезис басқарушы ағын теңдеуінің формасын өзгертпейді, бірақ гистерезиске қатысатын қасиеттер үшін құрылтай теңдеулерінің санын көбейтеді (көбейтеді).

1951-1970 жылдары коммерциялық компьютерлер ғылыми және инженерлік есептеулер мен модельдік модельдеу сахнасына шықты. Компьютерлік модельдеу мұнай қабаттарының динамикалық мінез-құлқы көп ұзамай мұнай өнеркәсібінің мақсатына айналды, бірақ ол кезде есептеу қуаты өте әлсіз болды.

Есептеу қабілетінің әлсіздігінде су қоймасының модельдері сәйкесінше өрескел болды, бірақ ұлғайту статикалық параметрлері өте қарапайым болды және өрескелдікке ішінара өтелді. Өзекше тығын масштабында алынған жыныстық қисықтардан (оны көбінесе микро шкала деп атайды) резервуар модельдерінің өрескел тор торларына дейін (көбінесе макро шкала деп аталады) салыстырмалы өткізгіштік қисықтарын жоғарылату туралы мәселе әлдеқайда қиын және ол әлі күнге дейін жалғасып келе жатқан маңызды зерттеу саласына айналды. Бірақ жоғарылатудағы прогресс баяу болды және 1990-2000 жылдарға дейін салыстырмалы өткізгіштік пен тензор ұсыну қажеттілігінің бағытқа тәуелділігі айқын көрсетіле алмады,[4][5] кем дегенде бір қабілетті әдіс болса да[6] 1975 жылы жасалған болатын. Осындай көтерілу жағдайларының бірі көлденең қозғалысқа қосымша мұнайға қатысты су (және газ) тігінен бөлінетін көлбеу қойма болып табылады. Тор ұяшығының тік өлшемі, сонымен қатар, тор ұяшығының көлденеңінен әлдеқайда аз, сәйкесінше кіші және үлкен ағын аймақтарын жасайды. Мұның бәрі х және z бағыттары үшін әр түрлі салыстырмалы өткізгіштік қисықтарын қажет етеді. Геологиялық гетерогенділік сияқты су қоймаларында ламиналар немесе көлденең өткізгіштік құрылымдары сонымен қатар бағытты салыстырмалы өткізгіштікті тудырады. Бұл бізге салыстырмалы өткізгіштік, ең көп жағдайда, тензормен ұсынылуы керек екенін айтады. Ағын теңдеулері болады

 Мұндағы a = w, o, g

Жоғарыда аталған жағдай көрініс тапты түсіру су айдау (немесе жаңартылған газ айдау) немесе қысымның төмендеуі арқылы өндіру. Егер сіз белгілі бір уақыт аралығында суды жаңартуды (немесе газдың төмендеуін) айдайтын болсаңыз, онда x + және x- бағыттары бойынша әр түрлі салыстырмалы өткізгіштік қисықтары пайда болады. Бұл дәстүрлі мағынада гистерезис процесі емес және оны дәстүрлі тензормен ұсынуға болмайды. Бұл бағдарламалық жасақтама кодында IF-операторымен ұсынылуы мүмкін және бұл кейбір коммерциялық резервуарлық тренажерлерде кездеседі. Процесс (дәлірек айтқанда, процестердің реттілігі) кен орнының қалпына келуіне арналған резервтік жоспарға байланысты болуы мүмкін немесе айдалатын сұйықтық су қоймасының күтпеген ашық бөлігінің салдарынан басқа қабат жыныстарының қабатына ағуы мүмкін. Кінә немесе а қаптаманың артында тығыздалмайтын цемент айдау ұңғымасының. Салыстырмалы өткізгіштік нұсқасы сирек қолданылады, және біз оның басқарушы теңдеудің өзгермейтінін (аналитикалық формасы) емес, көбейетінін (әдетте екі еселенетінін) атап өтеміз. құрылтай теңдеулері қатысты қасиеттер үшін.

Жоғарыда келтірілген теңдеу кеуекті ортадағы сұйықтық ағынының ең жалпы теңдеуінің векторлық формасы болып табылады және ол оқырманға қатысатын терминдер мен шамаларға жақсы шолу жасайды. Алдыға ауыспас бұрын дифференциалдық теңдеу ішіне айырымдық теңдеулер, компьютерлерде пайдалану үшін ағын теңдеуін компонент түрінде жазу керек. Компонент түріндегі ағын теңдеуі (қолдану арқылы) жиынтық конвенция ) болып табылады

 Мұндағы a = w, o, g  қайда = 1,2,3

Дарси жылдамдығы бұл сұйықтық бөлшегінің жылдамдығы емес, көлемдік ағын (жиі таңбамен ұсынылады) ) сұйықтық ағынының Кеуектердегі сұйықтық жылдамдығы (немесе қысқа, бірақ дұрыс емес деп аталады кеуек жылдамдығы) қатынасымен Darcy жылдамдығымен байланысты

 Мұндағы a = w, o, g

Көлемдік ағын - интенсивті шама, сондықтан бір уақытта қанша сұйықтық келетінін сипаттауда жақсы емес. Мұны түсіну үшін қолайлы көлем - бұл ағынның көлемділігі деп аталады, ол берілген уақыт ішінде берілген аймақтан қанша сұйықтықтың шығатынын (немесе шығатындығын) айтады және бұл қатынас бойынша Darcy жылдамдығымен байланысты.

 Мұндағы a = w, o, g

Ағынның көлемдік жылдамдығын байқаймыз скаляр шама болып табылады және бағытты беттің (ауданның) қалыпты векторы және көлемдік ағын (қамсыздандыру жылдамдығы) күтеді.

Резервуар моделінде геометриялық көлем торлы ұяшықтарға бөлінген, ал қазір қызығушылық аймағы екі көршілес ұяшықтардың қиылысу аймағы болып табылады. Егер бұл шынымен көршілес ұяшықтар болса, онда бұл аймақ жалпы бүйір беті болып табылады, ал егер ақаулар екі ұяшықты бөлсе, онда қиылысу ауданы әдетте екі көршілес ұяшықтардың толық бүйір бетінен аз болады. Осылайша, көп фазалы ағын теңдеуінің нұсқасы, оны резервуарлық тренажерларда қолданғанға дейін және оны қолданғанға дейін қолданылмайды

 Мұндағы a = w, o, g

Ол кеңейтілген (компоненттік) формада болады

 Мұндағы a = w, o, g

(Алғашқы) гидростатикалық қысым z тереңдікте (немесе деңгейде) z анықтамалық тереңдікте (немесе төменде)0 арқылы есептеледі

 Мұндағы a = w, o, g

Гидростатикалық қысымның есептеулері жүргізілген кезде, әдетте, фазалық индекс қолданылмайды, бірақ формуланы / шаманы нақты тереңдікте қандай фаза байқалатынына сәйкес ауыстырады, бірақ біз анықтық пен дәйектілік үшін фазалық индексті енгіздік. Алайда, гидростатикалық қысымды есептеу кезінде дәлдікті арттыру үшін тереңдікте өзгеретін ауырлық күшінің үдеуін қолдануға болады. Егер мұндай жоғары дәлдік қажет болмаса, онда ауырлық күшінің үдеуі тұрақты болып, есептелген қысым деп аталады артық қысым. Резервуар модельдеуінде мұндай жоғары дәлдіктің қажеті жоқ, сондықтан ауырлық күшінің үдеуі осы пікірталаста тұрақты болып саналады. Қойма моделіндегі бастапқы қысым (бастапқы) формуласы бойынша есептеледі артық қысым қайсысы

 Мұндағы a = w, o, g

Ағын теңдеуінің жақша ішіндегі шарттарды жеңілдету үшін ағын потенциалын енгізуге болады -потенциал, айқындалатын пси-потенциал, ол анықталады

 Мұндағы a = w, o, g

Ол екі шарттан тұрады, олар абсолютті қысым және ауырлық күші. Есептеу уақытын үнемдеу үшін интегралды бастапқыда есептеуге және кесте ретінде сақтауға болады, оны есептеуге қарағанда арзан кесте іздеуде қолдану керек. Енгізу -потенциал мұны білдіреді

 Мұндағы a = w, o, g

ПСИ-потенциалды көбінесе «сандық қысым» деп те атайды, өйткені функция қабаттағы жазықтыққа / тереңдікке ауысқаннан кейін қабаттағы кез-келген нүктедегі қысымды көрсетеді0. Тәжірибелік инженерлік жұмыста ұңғымаларда өлшенген қысымды санау деңгейіне жіберу немесе қабат қысымының таралуын картаға түсіру өте пайдалы. Осылайша, резервуардағы сұйықтықтың қозғалысының бағытын бір қарағанда көруге болады, өйткені қысым қысымының таралуы потенциалды үлестіруге тең. Екі қарапайым мысал мұны анықтайды. Су қоймасы тығыз тақтатас қабаттарымен бөлінген бірнеше ағынды блоктардан тұруы мүмкін. Бір су қоймасынан немесе ағынды блоктан шыққан сұйықтық ақаулыққа бір тереңдікте түсіп, екінші қабаттағы ағыннан немесе басқа тереңдіктегі ақаулықтан шығуы мүмкін. Сол сияқты сұйықтық бір ағынды қондырғыдағы өндірістік ұңғымаға түсіп, екінші ағынды блоктағы немесе резервуардағы өндірістік ұңғымадан шыға алады.

Кеуекті орта үшін көп фазалы ағын теңдеуі енді айналады

 Мұндағы a = w, o, g

Бұл көп фазалы ағын теңдеуі бағдарламалық жасақтамашы үшін дәстүрлі түрде мұнай өнеркәсібінде қолданылатын резервуарлық тренажерға арналған бағдарлама кодын жазу үшін теңдеуді дифференциалдық теңдеуден айырмашылық теңдеуге түрлендіре бастаған кезде бастама болды. Белгісіз тәуелді айнымалылар дәстүрлі түрде мұнай қысымы (мұнай кен орындары үшін) және тартылған сұйықтықтар үшін көлемдік шамалар болды, бірақ олардың біреуі мұнай қысымы және сұйықтық компоненттері үшін масса немесе моль шамалары үшін шешілетін модель теңдеулерінің жиынтығын қайта жазуы мүмкін.[7]

Жоғарыда келтірілген теңдеулер SI бірліктерінде жазылған және біз барлық материалдық қасиеттер SI бірліктерінде анықталған деп есептейміз. Бұдан шығатын нәтиже - теңдеулердің жоғарыда келтірілген нұсқаларында түрлендіруге арналған бірліктердің тұрақтылығы қажет емес. Мұнай өнеркәсібі әр түрлі қондырғыларды қолданады, олардың кем дегенде екеуі кең таралған. Егер сіз SI бірліктерінен басқа бірліктерді қолданғыңыз келсе, көпфазалы ағын теңдеулеріне арналған конверсияның тұрақты константаларын орнатуыңыз керек.

Бірліктерді түрлендіру

Жоғарыда келтірілген теңдеулер SI бірліктері (қысқа SI) деп басу бірлік D (дарси) үшін абсолютті өткізгіштік SI емес бірліктерде анықталады. Сондықтан бірлікке байланысты тұрақтылар жоқ. Мұнай өнеркәсібі SI қондырғыларын қолданбайды. Оның орнына олар SI қондырғыларының арнайы нұсқасын қолданады, оны біз қолданбалы SI бірліктері деп атаймыз немесе олар АҚШ пен Ұлыбританиядан шыққан өріс бірліктері деп аталатын басқа бірліктер жиынтығын пайдаланады. Температура теңдеулерге кірмейді, сондықтан біз оны қолдана аламыз фактор-белгі әдісі (егер бірлік-фактор әдісі деп те аталады), егер бізде H өлшем бірлігімен / параметрге ие болса, біз осы айнымалыны / параметрді конверсиялық тұрақтыға көбейтеміз, содан кейін айнымалы біз қалаған G бірлікке ие болады. Бұл дегеніміз біз H * C = G трансформациясын қолданамыз және өткізгіштік анықтамасының SI емес әсері өткізгіштік үшін конверсия коэффициентіне кіреді. H * C = G түрлендіруі барлық кеңістіктік өлшемдерге қатысты болады, сондықтан біз белгілерді ескермей, негізгі шарттарға назар аударамыз, содан кейін жақшаны гравитациялық мүшемен аяқтаймыз. Конверсияны бастамас бұрын, Дарсидің ағын теңдеуінің де (бір фазалы) және Мускат және басқалардың жалпыланған (немесе кеңейтілген) көп фазалы ағын теңдеулерін байқаймыз. резервуардың жылдамдығын (көлем ағыны), көлем жылдамдығын және тығыздықты қолданады. Осы шамалардың өлшем бірліктеріне с (немесе S) префиксін алатын стандартты беттік жағдайларда оларды аналогтарынан ажырату үшін r (немесе R) префиксі беріледі. Бұл теңдеулерді Өріс бірліктеріне айналдырған кезде өте маңызды. Бірліктерді түрлендіру қарапайым болып көрінетін тақырыпқа егжей-тегжейлі тоқталудың себебі, көптеген адамдар бірліктерді түрлендіру кезінде қателіктер жібереді.

Енді біз конверсия жұмысын бастауға дайынбыз. Алдымен біз теңдеудің ағынды нұсқасын алып, оны қалай жазамыз

 Мұндағы a = w, o, g

Біз композиттік конверсия коэффициентін өткізгіштік параметрімен бірге орналастырғымыз келеді. Мұнда біздің теңдеуіміз SI бірліктерінде жазылғанын, ал оң жағындағы айнымалылар / параметрлер тобы (бұдан әрі қысқаша параметрлер деп аталады) өлшемсіз топты құрайтынын ескереміз. Енді біз әрбір параметрді түрлендіреміз және осы түрлендірулерді бір конверсия тұрақтысына жинаймыз. Енді біз конверсиялық тұрақтылар тізімі (С) қолданылатын бірліктен SI бірлікке ауысатындығын және мұндай конверсия тізімдері үшін өте кең таралғанын ескереміз. Сондықтан біздің параметрлер қолданбалы бірліктерге енгізіліп, оларды (артқа) SI бірліктеріне айналдырады деп есептейміз.

 Мұндағы a = w, o, g

Өлшемсіз параметр болып табылатын салыстырмалы өткізгіштікті жойғанымызға назар аударыңыз. Бұл композициялық конверсия коэффициенті ағынның теңдеуі үшін Дарси тұрақтысы деп аталады және солай болады

Біздің параметр тобымыз SI базалық өлшемдерінде өлшемсіз болғандықтан, екінші кестеде көріп отырғанымыздай, SI бірліктерін біздің құрамдық конверсия коэффициенті үшін бірліктерге қосудың қажеті жоқ. Әрі қарай біз теңдеудің жылдамдық нұсқасын алып, оны қайта жазамыз

 Мұндағы a = w, o, g

Енді біз әрбір параметрді түрлендіреміз және осы түрлендірулерді бір конверсия тұрақтысына жинаймыз.

 Мұндағы a = w, o, g

Өлшемсіз параметр болып табылатын салыстырмалы өткізгіштікті жойғанымызға назар аударыңыз. Бұл композициялық конверсия коэффициенті ағынның теңдеуі үшін Дарси тұрақтысы деп аталады және солай болады

Қысым градиенті және ауырлық күші ағыны мен жылдамдық теңдеулері үшін бірдей, сондықтан бір рет қана талқыланады. Мұндағы міндет - қысым градиенті үшін қолданылатын бірліктерге («H-бірліктер») сәйкес келетін гравитациялық терминге ие болу. Сондықтан біз конверсия коэффициентін гравитациялық параметрлермен бірге орналастыруымыз керек. SI бірліктерінде «жақшаны» былай жазамыз

 Мұндағы a = w, o, g

және оны қайта жазыңыз

 Мұндағы a = w, o, g

Енді біз әрбір параметрді түрлендіреміз және осы түрлендірулерді бір конверсия тұрақтысына жинаймыз. Біріншіден, біздің теңдеуіміз SI бірліктерінде жазылғанын, ал оң жағындағы параметрлер тобы өлшемсіз топты құрайтынын ескереміз. Сондықтан, біздің параметрлер қолданбалы бірліктерге енгізілген деп санаймыз және оларды (артқа) SI бірліктеріне айналдырамыз.

 Мұндағы a = w, o, g

Бұл консистенция-түрлендіру үшін композициялық түрлендіру коэффициентін береді

Біздің параметрлер тобымыз SI бірліктерінде өлшемсіз болғандықтан, екінші кестеде көріп отырғанымыздай, SI бірліктерін біздің құрамдық конверсия коэффициенті үшін бірліктерге қосудың қажеті жоқ.

Бұл аналитикалық теңдеулерге арналған, бірақ бағдарламашы ағын теңдеуін ақырлы айырым теңдеуіне, одан әрі сандық алгоритмге айналдырғанда, олар есептеу амалдарының санын барынша азайтуға асық. Мұнда термоядролық реакцияның көмегімен біріне дейін төмендетуге болатын екі тұрақтысы бар мысал келтірілген

Салалық бірліктерді қолдана отырып, векторлық түрдегі ағын теңдеуінің ағындық нұсқасы болады

 Мұндағы a = w, o, g

және компонент түрінде ол болады

 Мұндағы a = w, o, g  қайда = 1,2,3

Салалық бірліктерді қолдана отырып, векторлық түрдегі ағын теңдеуінің жылдамдық нұсқасы болады

 Мұндағы a = w, o, g

және компонент түрінде ол болады

 Мұндағы a = w, o, g

Бөлшектерді түрлендіру - бұл тіпті техникалық мамандар үшін де сирек кездесетін қызмет, бірақ бұл адамдар оны қалай дұрыс жасау керектігін ұмытып кетудің себебі болып табылады.

Сондай-ақ қараңыз

Әдебиеттер тізімі

  1. ^ Muskat M. and Meres M.W. 1936. Кеуекті орта арқылы гетерогенді сұйықтық ағыны. Дж. Апплда жарияланған мақала. Физ. 1936, 7, 346–363 бб. https://dx.doi.org/10.1063/1.1745403
  2. ^ Muskat M. and Wyckoff RD, Botset H.G. and Meres M.W. 1937. Газ-сұйық қоспалардың құм арқылы өтуі. 1937 жылғы AIME операциялары, 123, 69–96 беттерінде жарияланған. SPE құжатының идентификаторы - SPE-937069-G. https://dx.doi.org/10.2118/937069-G
  3. ^ Леверетт 1941. Кеуекті қатты денелердегі капиллярлық тәртіп. AIME 1940 ж. Тулса жиналысында ұсынылған мақала. Транзакцияларда AIME 1941 ж., 142, 159–172 бб. SPE құжатының идентификаторы - SPE-941152-G. https://dx.doi.org/10.2118/941152-G
  4. ^ Pickup G.E және Sorbie K.S. 1996. Фазалық өткізгіштік тензорларын қолданатын кеуекті ортадағы екі фазалы ағынның масштабы. SPE-28586 қағазы алғаш рет АҚШ-тың Жаңа Орлеан қаласында өткен SPE Жыл сайынғы техникалық конференциясында және көрмесінде ұсынылды, 1994 ж. 25-28 қыркүйек. SPEJ-28586-PA Қағаз 1996 жылдың желтоқсан айының SPEJ-де жарық көрді. https://dx.doi.org/10.2118/28586-PA
  5. ^ Kumar A.T.A. және Джераулд Г.Р. 1996 ж. Масштабтың су қоймасындағы тастан блокқа дейінгі масштабтан сұйықтық ағынына әсері. SPE / DOE 35452 қағазы 1996 жылы 21-24 сәуірде АҚШ-тың Оклахома штатындағы Тулса қаласында өткен жақсартылған мұнай қалпына келтіру бойынша SPE / DOE 10-шы симпозиумында ұсынылды. SPE-35452-MS қағазы. https://dx.doi.org/10.2118/35452-MS
  6. ^ Kyte JR және Berry D.W. 1975. Сандық дисперсияны басқарудың жаңа жалған функциялары. SPE 5105 қағазы алғаш рет 1974 жылы 6-9 қазан аралығында Хьюстон қаласында өткен SPE-AIME 49-шы күзгі жиналысында ұсынылды. SPEJ 5105-PA қағазы SPEJ 1975 ж. Тамызында 269-275 беттерінде шығарылды. https://dx.doi.org/10.2118/5105-PA
  7. ^ Жас Л.С. және Stephenson R.E. 1983. Су қоймаларын модельдеуге арналған жалпыланған композициялық тәсіл. SPE 10516 қағазы алғаш рет Жаңа Орлеанда өткен SPE су қоймасының симпозиумында, 31 қаңтар - 3 ақпан 1982 ж. Ұсынылды. SPE-10516-PA қағазы алғаш SPEJ 1983 ж. Қазанында басылды, 23 том, № 05, 727-742 бет; AIME 275 транзакциясы; қазір ResearchGate арқылы; https://www.researchgate.net/publication/244956766 және арқылы https://dx.doi.org/10.2118/10516-PA