【摘要】：裂縫內(nèi)可壓縮流動(dòng)廣泛存在于油氣開采、CO_2埋存、地下水污染處理等工業(yè)生產(chǎn)及自然界中,是能源開發(fā)、環(huán)境治理等領(lǐng)域重點(diǎn)關(guān)注的問題。例如,在天然氣開采中的高壓氣井井筒附近,氣體的可壓縮性對(duì)流動(dòng)的影響顯著。本文將基于氣體動(dòng)理學(xué)格式(Gas Kinetic Scheme,GKS)在孔隙尺度上模擬裂縫內(nèi)可壓縮流動(dòng),并分析流動(dòng)達(dá)到阻塞的影響因素,研究基質(zhì)巖塊系統(tǒng)與裂縫系統(tǒng)之間的竄流對(duì)出流流量的影響。以下是本文研究工作的主要內(nèi)容:1、結(jié)合直角網(wǎng)格技術(shù),開發(fā)了針對(duì)多孔介質(zhì)內(nèi)流動(dòng)的GKS模擬方法。檢驗(yàn)了該數(shù)值方法模擬復(fù)雜區(qū)域的不可壓低速滲流和可壓縮流動(dòng)的可行性與有效性。測(cè)試結(jié)果表明:開發(fā)的GKS程序能夠有效地處理復(fù)雜邊界形狀;相比于單松弛格子玻爾茲曼方法,GKS能夠精確實(shí)現(xiàn)壁面無(wú)滑移邊界條件,從而正確反映滲透率與黏性無(wú)關(guān)的特性;對(duì)于Berea砂巖切片結(jié)構(gòu)中的復(fù)雜流動(dòng),模擬結(jié)果與實(shí)驗(yàn)吻合較好,能較精確地計(jì)算出滲透率。2、模擬了曲張管道內(nèi)的可壓縮流動(dòng),分析了管道的各項(xiàng)特征參數(shù)(如管道長(zhǎng)度,凸起物的幅度,凸起物個(gè)數(shù)等)及管道截面形狀對(duì)管內(nèi)可壓縮流動(dòng)的影響,總結(jié)了管內(nèi)流動(dòng)發(fā)生阻塞的相關(guān)規(guī)律:對(duì)于該類曲張管道模型,管內(nèi)阻塞的位置均是位于喉道處,且管道越短或凸起幅度越大或凸起物越多,管道越容易發(fā)生阻塞;臨界壓力比與管道長(zhǎng)度存在線性關(guān)系;各項(xiàng)參數(shù)相同時(shí),圓管道比平行通道更容易發(fā)生阻塞。3、提出了與氣藏儲(chǔ)層的平均壓力、氣藏儲(chǔ)層的滲透率、裂縫尺寸、氣體動(dòng)力學(xué)黏度等相關(guān)的基質(zhì)-裂縫竄流邊條件,并將其應(yīng)用到直圓管內(nèi)可壓縮流動(dòng)的研究中。模擬結(jié)果表明:流入管內(nèi)的氣體總量是氣藏儲(chǔ)層的平均壓力與井筒壓力比值的非線性函數(shù),且存在飽和現(xiàn)象;當(dāng)氣藏儲(chǔ)層的平均壓力與井筒壓力比值較大時(shí),井筒中的大部分氣體都產(chǎn)自裂縫與井筒連接處附近的儲(chǔ)層區(qū)域�？傊�,本文結(jié)合直角網(wǎng)格技術(shù),開發(fā)了針對(duì)多孔介質(zhì)流動(dòng)的GKS模擬方法,檢驗(yàn)了該數(shù)值方法模擬復(fù)雜區(qū)域的不可壓低速滲流和可壓縮流動(dòng)的可行性與有效性,為進(jìn)一步研究多孔介質(zhì)內(nèi)的輸運(yùn)問題提供了有效的工具;此外,本文還提出了一種基質(zhì)-裂縫竄流邊條件,分析了井筒附近的流道四周存在竄流時(shí),對(duì)流道內(nèi)可壓縮流動(dòng)的影響,為以后的相關(guān)研究奠定了基礎(chǔ)。
【圖文】：

華 中 科 技 大 學(xué) 碩 士 學(xué) 位 論 文函數(shù)求矩可得： +1= 1 Δ + . (2.8)因此，為了使單元的宏觀守恒量在每個(gè)離散時(shí)間步中得到更新，我們首先必須構(gòu)造出各個(gè)界面處的數(shù)值通量 。

= 2+2 2 1 2 2+12 3 3 4 , 12( 2+ 2+ ( + 2)/2 )、 2=12( 3+ 2+ ( + 4) /2 )、 4=14( 4+ 4+ 2 2 2+ ( )，則將(2.13)、(2.14)、(2.15)、(2.16)、(2.17的值。偏導(dǎo)數(shù)均可根據(jù)宏觀量的偏導(dǎo)數(shù)求得，而宏觀量(2.11)中的 0可根據(jù)公式(2.4)通過界面處的宏可壓縮問題，界面處的宏觀量可通過界面左右線性插值與中心差分的模板如圖 2.2 所示，所考
【學(xué)位授予單位】：華中科技大學(xué)
【學(xué)位級(jí)別】：碩士
【學(xué)位授予年份】：2019
【分類號(hào)】：TE37;O35

【相似文獻(xiàn)】

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

1 祁大同,王尚錦;奇點(diǎn)法在葉輪機(jī)械任意回轉(zhuǎn)面亞音速可壓縮流動(dòng)計(jì)算中的應(yīng)用[J];應(yīng)用力學(xué)學(xué)報(bào);1989年02期

2 王爾智;陶瑞民;;求解非定�？蓧嚎s流動(dòng)問題的邊界元方法的研究[J];中國(guó)電機(jī)工程學(xué)報(bào);2006年24期

3 馬東軍,孫德軍,尹協(xié)遠(yuǎn);一維多介質(zhì)可壓縮流動(dòng)數(shù)值方法[J];計(jì)算物理;2003年02期

4 徐麗;姜明洋;蔡靜靜;吳碩;;無(wú)粘可壓縮流動(dòng)的改進(jìn)高精度方法[J];計(jì)算力學(xué)學(xué)報(bào);2019年06期

5 馬東軍,蔡勇,孫德軍,尹協(xié)遠(yuǎn);兩種多介質(zhì)流體可壓縮流動(dòng)界面捕捉方法的數(shù)值研究(英文)[J];中國(guó)科學(xué)技術(shù)大學(xué)學(xué)報(bào);2002年02期

6 杜東興,李志信,過增元,岳敏;微細(xì)光滑管內(nèi)可壓縮流動(dòng)換熱特性的數(shù)值研究[J];清華大學(xué)學(xué)報(bào)(自然科學(xué)版);2000年11期

7 丁文靜,何雅玲,陶文銓;基于脈管制冷機(jī)的二維可壓縮流動(dòng)數(shù)值計(jì)算程序的開發(fā)與應(yīng)用[J];工程熱物理學(xué)報(bào);2004年06期

8 馬東軍,孫德軍,尹協(xié)遠(yuǎn);高密度比多介質(zhì)可壓縮流動(dòng)的PPM方法[J];計(jì)算物理;2001年06期

9 張鐳,袁禮;用改進(jìn)的耦合型Level Set方法計(jì)算一維雙介質(zhì)可壓縮流動(dòng)[J];計(jì)算物理;2001年06期

10 孟慶昌;張志宏;顧建農(nóng);劉巨斌;;回轉(zhuǎn)橢球體可壓縮流動(dòng)數(shù)值計(jì)算[J];武漢理工大學(xué)學(xué)報(bào)(交通科學(xué)與工程版);2008年06期

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

1 羅喜勝;;可壓縮流動(dòng)中水蒸汽的同質(zhì)凝結(jié)的理論分析與數(shù)值模擬[A];第七屆中國(guó)空氣動(dòng)力學(xué)學(xué)會(huì)物理氣體動(dòng)力學(xué)專業(yè)委員會(huì)換屆選舉大會(huì)暨學(xué)術(shù)交流會(huì)論文集[C];2016年

2 李萬(wàn)愛;任玉新;;非結(jié)構(gòu)無(wú)波動(dòng)間斷譜元法在可壓縮流動(dòng)中的應(yīng)用[A];中國(guó)力學(xué)大會(huì)-2015論文摘要集[C];2015年

3 袁湘江;陶建軍;沈清;涂國(guó)華;;可壓縮流動(dòng)中小激波形成的機(jī)理研究[A];慶祝中國(guó)力學(xué)學(xué)會(huì)成立50周年暨中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)’2007論文摘要集（下）[C];2007年

4 陸夕云;;可壓縮流動(dòng)的數(shù)值計(jì)算和分析[A];中國(guó)力學(xué)大會(huì)-2015論文摘要集[C];2015年

5 姜元勇;徐曾和;;局部非熱平衡條件下填充床中的可壓縮流動(dòng)與氣固反應(yīng)[A];慶祝中國(guó)力學(xué)學(xué)會(huì)成立50周年暨中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)’2007論文摘要集（下）[C];2007年

6 范進(jìn)之;李樺;;基于Guauss-Lobatto節(jié)點(diǎn)的CPR方法在求解可壓縮流動(dòng)問題時(shí)計(jì)算精度的研究[A];中國(guó)力學(xué)大會(huì)-2015論文摘要集[C];2015年

7 姜元勇;徐曾和;;填充床中可壓縮流動(dòng)與非等溫氣固反應(yīng)[A];慶祝中國(guó)力學(xué)學(xué)會(huì)成立50周年暨中國(guó)力學(xué)學(xué)會(huì)學(xué)術(shù)大會(huì)’2007論文摘要集（下）[C];2007年

8 朱建州;楊焱;;氣動(dòng)聲學(xué)中不同模態(tài)的轉(zhuǎn)換以及關(guān)于聲模態(tài)和聲源的界定與分離的再思考[A];第九屆全國(guó)流體力學(xué)學(xué)術(shù)會(huì)議論文摘要集[C];2016年

9 任虎;朱建州;;渦輪機(jī)內(nèi)旋轉(zhuǎn)流動(dòng)相關(guān)問題研究[A];2018年全國(guó)工業(yè)流體力學(xué)會(huì)議摘要集[C];2018年

10 李偉鵬;劉洪;;高精度激波捕獲格式的研究與分析[A];高超聲速專題研討會(huì)暨第五屆全國(guó)高超聲速科學(xué)技術(shù)會(huì)議論文集[C];2012年

相關(guān)博士學(xué)位論文 前7條

1 田浩;三維可壓縮流動(dòng)的非定常分離及飛行器運(yùn)動(dòng)/流動(dòng)耦合的非線性動(dòng)態(tài)穩(wěn)定性[D];中國(guó)空氣動(dòng)力研究與發(fā)展中心;2015年

2 王芳;非緊致結(jié)構(gòu)氣動(dòng)噪聲數(shù)值積分預(yù)測(cè)方法研究[D];西北工業(yè)大學(xué);2015年

3 徐浩然;大型風(fēng)力機(jī)翼型氣動(dòng)性能優(yōu)化研究[D];揚(yáng)州大學(xué);2015年

4 林健宇;切割網(wǎng)格方法及激波與氣泡相互作用研究[D];中國(guó)科學(xué)技術(shù)大學(xué);2016年

5 程萬(wàn);伴隨凝結(jié)的無(wú)粘可壓縮流動(dòng)研究[D];中國(guó)科學(xué)技術(shù)大學(xué);2011年

6 姜元勇;多孔相變材料中的流動(dòng)與反應(yīng)[D];東北大學(xué);2009年

7 李德波;氣固兩相剪切流動(dòng)大規(guī)模并行直接數(shù)值模擬研究[D];浙江大學(xué);2011年

相關(guān)碩士學(xué)位論文 前8條

1 曾偉;裂縫內(nèi)可壓縮流動(dòng)的氣體動(dòng)理學(xué)數(shù)值模擬研究[D];華中科技大學(xué);2019年

2 王博;二維可壓縮流動(dòng)轉(zhuǎn)捩計(jì)算[D];西北工業(yè)大學(xué);2004年

3 左計(jì)學(xué);應(yīng)用場(chǎng)協(xié)同理論的蒸汽噴射器優(yōu)化設(shè)計(jì)研究[D];大連理工大學(xué);2012年

4 李震;基于有限差分方法的可壓縮流動(dòng)大渦數(shù)值模擬[D];哈爾濱工業(yè)大學(xué);2014年

5 王杰;高溫堆直接循環(huán)系統(tǒng)瞬態(tài)分析算法研究[D];哈爾濱工程大學(xué);2016年

6 劉宇;基于渦方法的渦運(yùn)動(dòng)研究[D];北京交通大學(xué);2011年

7 袁瑞峰;氣體動(dòng)理論格式的浸入邊界法研究[D];西北工業(yè)大學(xué);2015年

8 季鑫龍;可壓縮流動(dòng)的計(jì)算及AUSM+格式的應(yīng)用[D];華東理工大學(xué);2014年

，

本文編號(hào)：2590273