【摘要】：非均勻多孔介質(zhì)中的單相穩(wěn)態(tài)滲流可以用類拉普拉斯方程來描述,準(zhǔn)確求解該方程是油藏?cái)?shù)值模擬最基礎(chǔ)的問題之一。地下油氣儲(chǔ)層往往具有強(qiáng)非均質(zhì)性,在計(jì)算強(qiáng)非均質(zhì)地層中的單相穩(wěn)態(tài)滲流時(shí),如果滲流率為標(biāo)量,傳統(tǒng)數(shù)值方法將網(wǎng)格界面平均滲透率取為相鄰網(wǎng)格滲透率的各種平均值,例如調(diào)和平均值、幾何平均值或算術(shù)平均值等,但是由于對(duì)其中數(shù)學(xué)問題機(jī)制的理解不清晰,這些傳統(tǒng)數(shù)值方法都不具備通用性。該問題的本質(zhì)在于不同滲透率網(wǎng)格界面交點(diǎn)處滲流速度和壓力梯度存在發(fā)散現(xiàn)象。對(duì)于滲透率為張量形式的情形,傳統(tǒng)數(shù)值方法同樣基于調(diào)和平均的思想來處理非均質(zhì)滲流問題,因此同樣面臨精度差的困境。對(duì)于具有張量形式滲透率的單相穩(wěn)態(tài)滲流,在不同滲透率網(wǎng)格界面的交點(diǎn)處同樣存在滲流速度和壓力梯度發(fā)散的現(xiàn)象。無論是在標(biāo)量滲透率還是張量滲透率的情形下,上述發(fā)散現(xiàn)象的相關(guān)研究工作都是在規(guī)則的矩形網(wǎng)格中進(jìn)行的,本文將研究非規(guī)則網(wǎng)格中的相關(guān)發(fā)散現(xiàn)象,并在二維非結(jié)構(gòu)化網(wǎng)格中建立用于計(jì)算具有張量形式滲透率的非均勻介質(zhì)單相穩(wěn)態(tài)滲流的有限分析數(shù)值格式。本文首先推導(dǎo)出具有張量形式滲透率的二維類拉普拉斯方程在任意形狀角域的局部?jī)缏山馕鼋?該解析解的冪指數(shù)由角域的滲透率分布決定,與邊界條件無關(guān),具有內(nèi)稟特性。通過理論推導(dǎo)發(fā)現(xiàn),角域流動(dòng)存在三種基本流動(dòng)模式:冪律流動(dòng),線性流動(dòng)以及滯止流。并且在同一角域,不同流動(dòng)模式可以共存,角域冪律解析解則是這些流動(dòng)模式所對(duì)應(yīng)基本解的線性組合。從任意形狀角域的冪律解析解出發(fā),建立了二維非結(jié)構(gòu)化網(wǎng)格中求解具有張量形式滲透率的單相穩(wěn)態(tài)滲流的有限分析數(shù)值格式。數(shù)值算例表明,該數(shù)值方法精度很高,而且隨著網(wǎng)格加密,其向真值的收斂速度比傳統(tǒng)方法快得多,在實(shí)際應(yīng)用中,推薦使用細(xì)化參數(shù)為n = 2或n = 3的網(wǎng)格進(jìn)行計(jì)算,此時(shí)等效滲透率的誤差將低于4%,并且與介質(zhì)非均質(zhì)性的強(qiáng)弱無關(guān)。相比之下,在計(jì)算強(qiáng)非均質(zhì)問題時(shí),傳統(tǒng)算法要得到準(zhǔn)確的結(jié)果,則需要將網(wǎng)格充分細(xì)分;且隨著介質(zhì)非均質(zhì)性強(qiáng)度的增加,所需要網(wǎng)格細(xì)分程度將會(huì)相應(yīng)劇增。另外,本文還得到了非均勻多孔介質(zhì)中二維單相穩(wěn)態(tài)滲流方程的解析解,該解析解可以表示為一個(gè)無窮冪級(jí)數(shù)。針對(duì)每個(gè)奇點(diǎn)角域都能計(jì)算得到一組固有的冪指數(shù),結(jié)合給定的邊界條件,每個(gè)冪級(jí)數(shù)項(xiàng)前的系數(shù)可以通過數(shù)值方法確定,從而得到相應(yīng)的精確解。本文構(gòu)造的有限分析數(shù)值格式和得到的冪級(jí)數(shù)解析解是針對(duì)二維非均勻多孔介質(zhì)中的單相穩(wěn)態(tài)滲流提出的。由于相應(yīng)的二維類拉普拉斯方程同樣適用于非均勻介質(zhì)中的其他擴(kuò)散問題,如熱傳導(dǎo)、靜電場(chǎng)等,因此本文提出的有限分析數(shù)值格式和得到的解析解同樣可應(yīng)用于相關(guān)領(lǐng)域的研究,具有廣闊的應(yīng)用前景。
【圖文】：

２．２邐二維任意形狀角域中類拉普拉斯方程的冪律解析解逡逑２．２．１角點(diǎn)附近壓力分布的局域性態(tài)逡逑如圖２．１所示的一個(gè)任意形狀角域，整個(gè)平面劃分為了邋７Ｖ個(gè)子角域，，各個(gè)子邐、逡逑角域的分界面在角點(diǎn)處交匯。在每一個(gè)子角域中，壓力方程，即類拉普拉斯方程逡逑（２．４）

２．２邐二維任意形狀角域中類拉普拉斯方程的冪律解析解逡逑２．２．１角點(diǎn)附近壓力分布的局域性態(tài)逡逑如圖２．１所示的一個(gè)任意形狀角域，整個(gè)平面劃分為了邋７Ｖ個(gè)子角域，各個(gè)子邐、逡逑角域的分界面在角點(diǎn)處交匯。在每一個(gè)子角域中，壓力方程，即類拉普拉斯方程逡逑（２．４）
【學(xué)位授予單位】：中國科學(xué)技術(shù)大學(xué)
【學(xué)位級(jí)別】：博士
【學(xué)位授予年份】：2018
【分類號(hào)】：O357.3

【相似文獻(xiàn)】

相關(guān)期刊論文 前10條

1 臥龍;;滲透率預(yù)示互聯(lián)網(wǎng)股大調(diào)整(下)[J];股市動(dòng)態(tài)分析;2018年14期

2 季小凱;郭建斌;邢同菊;譚希鵬;;煤系沉積巖應(yīng)力 應(yīng)變與應(yīng)變 滲透率特征[J];煤田地質(zhì)與勘探;2015年03期

3 楊曉晶;;高端有機(jī)成嬰幼兒配方奶粉發(fā)展趨勢(shì)[J];中國食品;2017年12期

4 盧紹群;龔根平;;評(píng)估電壓偏差與諧波約束下配網(wǎng)光伏最大滲透率[J];科學(xué)家;2017年02期

5 邵維志;解經(jīng)宇;遲秀榮;李俊國;吳淑琴;肖斐;;低孔隙度低滲透率巖石孔隙度與滲透率關(guān)系研究[J];測(cè)井技術(shù);2013年02期

6 車廷信;王學(xué)武;田文博;戰(zhàn)劍飛;韓雪;;滲流介質(zhì)對(duì)滲透率測(cè)試結(jié)果的影響[J];科技導(dǎo)報(bào);2012年35期

7 賀青玉;;3D電視增長(zhǎng)顯著[J];家用電器;2011年05期

8 張進(jìn)軍;昝軍才;;采空區(qū)自由堆積滲透率假設(shè)及應(yīng)用[J];陜西煤炭;2011年04期

9 朱春明;鄧金根;王利華;王希玲;韋龍貴;何寶生;曾祥林;;防砂管外堆積砂滲透率實(shí)驗(yàn)研究[J];石油鉆采工藝;2010年01期

10 劉全穩(wěn);;關(guān)于極差與滲透率極差及其他[J];西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版);2009年01期

相關(guān)會(huì)議論文 前10條

1 胡恒山;王軍;;基于動(dòng)電效應(yīng)的巖芯滲透率實(shí)驗(yàn)測(cè)量[A];中國地球物理2010——中國地球物理學(xué)會(huì)第二十六屆年會(huì)、中國地震學(xué)會(huì)第十三次學(xué)術(shù)大會(huì)論文集[C];2010年

2 閆國亮;孫建孟;;不同尺度下滲透率模型研究進(jìn)展[A];“地球·資源”全國博士生學(xué)術(shù)論壇會(huì)議論文摘要集[C];2011年

3 張德會(huì);李健康;;多孔介質(zhì)滲透率在金屬成礦學(xué)研究中的意義[A];固體地球系統(tǒng)復(fù)雜性與地質(zhì)過程動(dòng)力學(xué)學(xué)術(shù)討論會(huì)論文摘要集[C];2004年

4 常莉;謝馥勵(lì);張玉君;王克協(xié);;利用聲偶極測(cè)井資料求取各向異性儲(chǔ)層衰減與滲透率研究[A];中國地球物理學(xué)會(huì)第二十四屆年會(huì)論文集[C];2008年

5 盧穎忠;李銘;劉福利;昌慶珍;;大慶油田滲透率分布特征研究[A];1999年中國地球物理學(xué)會(huì)年刊——中國地球物理學(xué)會(huì)第十五屆年會(huì)論文集[C];1999年

6 張亞蒲;何應(yīng)付;楊正明;鮮保安;;煤層氣藏滲透率應(yīng)力敏感性研究[A];煤層氣勘探開發(fā)理論與實(shí)踐[C];2007年

7 魏斌;;利用流動(dòng)單元計(jì)算高含水油田滲透率[A];2001年中國地球物理學(xué)會(huì)年刊——中國地球物理學(xué)會(huì)第十七屆年會(huì)論文集[C];2001年

8 曹成;李天太;劉剛;高潮;王宇;;考慮吸附效應(yīng)的頁巖基質(zhì)滲透率計(jì)算模型[A];2015年全國沉積學(xué)大會(huì)沉積學(xué)與非常規(guī)資源論文摘要集[C];2015年

9 杜堯俊;李文延;;超高圍壓與超高滲透力作用下砂巖滲透率的實(shí)驗(yàn)研究[A];第五屆全國構(gòu)造物理、第三屆全國高溫高壓學(xué)術(shù)討論會(huì)論文摘要[C];1992年

10 李閩;;非線性有效壓力計(jì)算[A];中國力學(xué)大會(huì)——2013論文摘要集[C];2013年

相關(guān)重要報(bào)紙文章 前10條

1 每經(jīng)記者 金U

本文編號(hào)：2652703