【正文】 6 180m 處煤層壓力隨時(shí)間分布Fig 56 Pressure in the 180m by times表 54 180m 處隨時(shí)間變化的最大壓力值時(shí)間（s） 0 1 1—30壓力值（MPa） 圖（a1）開切眼處壓力梯度（0s）：深部煤層群采動(dòng)影響下瓦斯運(yùn)移規(guī)律的數(shù)值模擬研究32圖（a3）開切眼處壓力梯度（30s）圖（b1）全面開挖 180m 壓力梯度（0s）圖（b2）全面開挖 180m 壓力梯度（30s）圖 57 被卸壓煤層壓力梯度分布Fig 57 Pressure gradient in the coal從圖 56 可以明顯發(fā)現(xiàn)，由于工作面開采位置的延伸，由于應(yīng)力狀態(tài)的不同，透氣性的差異，開切眼附近的瓦斯壓力降低要快。隨著工作面的延伸，被卸瓦斯壓力也逐步降低。隨著工作面的推進(jìn)，壓力是呈動(dòng)態(tài)變化的，最后趨于一個(gè)固定值。工作面過開切眼后，瓦斯壓力逐漸下降，待巖層移動(dòng)基本穩(wěn)定時(shí)。采空區(qū)中部的裂隙被上覆巖層壓實(shí)，壓力也趨于一個(gè)固定值。隨著壓力的變化和工作面的推進(jìn)。覆巖活動(dòng)中形成穩(wěn)定的結(jié)構(gòu)，（論文）33從而使冒落位移增大。煤體所承受的壓力降低，瓦斯流量增加位移增大，壓力降低發(fā)生變形。開采過程中被卸壓煤層開切眼周圍的瓦斯壓力梯度的變化情況,在開切眼周圍形成了明顯的瓦斯壓力梯度分布，并且隨著煤層的開采、瓦斯的卸壓，瓦斯壓力梯度分布范圍逐漸在擴(kuò)大，這是由于瓦斯壓力值逐漸減小造成的。瓦斯壓力梯度越小，滲透率也越小，這是由于瓦斯壓力對(duì)煤的孔隙和裂隙有促進(jìn)和擴(kuò)展作用。一般情況下，煤層滲透率隨壓力（深度）的增加而減少，壓力越大，孔隙裂隙趨于閉合使?jié)B透率減小。X 方向壓力梯度： ecpPx)(???? y 方向壓力梯度： (54)nypN)(由公式可以明顯看出：煤層煤體滲透率與壓力梯度成正比。煤層開采后，覆巖與開采層的采空區(qū)之間形成瓦斯壓力梯度。 6 結(jié)論與展望 結(jié)論 本文針對(duì)高瓦斯、低透氣性的深部煤層群煤層瓦斯的賦存條件，運(yùn)用理論分析、數(shù)值模擬等手段，系統(tǒng)的研究了卸壓開采上覆巖層采動(dòng)裂隙的卸壓煤層變形、采動(dòng)應(yīng)力變化、壓力梯度變化及瓦斯運(yùn)移規(guī)律等。根據(jù)開采后位移分布、應(yīng)力分布的不同，煤層開采后，沿煤層走向得出煤巖體裂隙分布規(guī)律。豎向依次為煤壁支撐影響區(qū)、離層區(qū)、重新壓實(shí)區(qū)。橫向依次為冒落帶、裂隙帶、彎曲下沉帶。在停采煤柱、開切眼附近形成應(yīng)力增高閉合區(qū)；遠(yuǎn)離開切眼、工作面應(yīng)力增加平緩，形成應(yīng)力平緩裂隙未變區(qū)；采空區(qū)上方形成充分卸壓裂隙發(fā)育區(qū)，該區(qū)可分為卸壓波動(dòng)區(qū)、卸壓增大區(qū)以及卸壓緩慢變化區(qū)。開采過程中，卸壓瓦斯壓力值逐漸變小到一定程度趨于一定值，壓力梯度逐漸變大，瓦斯?jié)B透率增大。：深部煤層群采動(dòng)影響下瓦斯運(yùn)移規(guī)律的數(shù)值模擬研究34 展望 煤與瓦斯共采技術(shù)，在開采高瓦斯煤層的同時(shí)，利用巖層運(yùn)動(dòng)的特點(diǎn)將煤層瓦斯開采出來將是煤層瓦斯的重要途徑。在井下因采動(dòng)影響地層壓力發(fā)生變化，由于開采卸壓，煤層中的瓦斯壓力升高，煤中的原來的孔隙系統(tǒng)反而降低。有效地進(jìn)行井下瓦斯抽放和地面煤層氣的開發(fā)，是煤與瓦斯共采的關(guān)鍵技術(shù)問題，巖層運(yùn)動(dòng)的關(guān)鍵層理論。因此在提高瓦斯率研究的同時(shí)，也應(yīng)該注重將井下瓦斯抽取，形成一套屬于煤與瓦斯共采的基礎(chǔ)理論和技術(shù)體系，從而真正實(shí)現(xiàn)了煤與瓦斯共采，更好的為煤礦安全生產(chǎn)服務(wù)。本文不足在于由于時(shí)間問題和本人學(xué)習(xí)能力有限，未能數(shù)值模擬出煤巖體中的離層裂隙以及沒有在三維中更加生動(dòng)的模擬分析。致謝經(jīng)過一學(xué)期的忙碌，大學(xué)最后的任務(wù)——畢業(yè)論文接近尾聲，回首和論文一起“走”過的日子，感慨良多。做論文的過程中有對(duì)新知識(shí)汲取的新鮮感，也有研究不懂時(shí)候的煩惱感。這個(gè)過程中如果沒有老師和同學(xué)們的鼓勵(lì)幫助，恐怕很難完成本次畢業(yè)論文。由于本論文題目中涉及到大量的信息量與知識(shí)量、知識(shí)點(diǎn)。此中我的論文導(dǎo)師梁冰老師給予了我最大的鼓勵(lì)與幫助，不厭其煩的耐心講解每一個(gè)細(xì)節(jié)，有問必答，為我確定完成論文的目標(biāo)與進(jìn)度。完成過程中，我見識(shí)到了他的學(xué)識(shí)與很強(qiáng)的理論基礎(chǔ)與應(yīng)用實(shí)際能力。在這里我懷著無比感激的心對(duì)老師道一聲感謝，謝謝您！同時(shí)我要感謝我的父母，一路的求學(xué)之路，如若不是你們強(qiáng)烈的后盾支持與鼓勵(lì)就沒有現(xiàn)在的我，能本次完成論文也離不開父母給我的關(guān)心。最后，畢業(yè)來臨之際，我想感謝我的母校，遼寧工程技術(shù)大學(xué)力學(xué)與工程學(xué)院對(duì)我的栽培。（論文）35 潘妍 2022 年 6 月 9 日參考文獻(xiàn)[1]煤炭科學(xué)院北京開采所，煤礦地表移動(dòng)與覆巖破斷規(guī)律及其應(yīng)用[M]北京：煤炭工業(yè)出版社，1981[2][J].煤炭學(xué)報(bào),1995,20(1):15[3]錢銘高,[M].北京:煤炭工業(yè)出版社,[4]、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程[M].北京:煤炭工業(yè)[5][J].煤炭學(xué)報(bào)，1996,21(6):576581[6]楊倫,于廣明，[J].煤炭學(xué)報(bào),1997,22(5):477480[7]“四帶”模型與動(dòng)態(tài)位移及分析[J],煤炭學(xué)報(bào),1996,20(1):5153[8]滕云海,[J],煤炭學(xué)報(bào),1999,24(1):2627[9][J].煤炭學(xué)報(bào),1995,20(1):3942：深部煤層群采動(dòng)影響下瓦斯運(yùn)移規(guī)律的數(shù)值模擬研究36[10]王金莊，康建榮，[J].中國(guó)礦業(yè)大學(xué)學(xué)報(bào),(4):331334[11]錢鳴高,繆協(xié)興,[J].煤炭學(xué)報(bào),1996,21(3):225230[12][D].徐州:中國(guó)礦業(yè)大學(xué),1999[13]錢鳴高,繆協(xié)興,[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,2022[14][M]徐州:中國(guó)礦業(yè)大學(xué)出版社,2022[15]程遠(yuǎn)平,俞啟香,[J]中國(guó)礦業(yè)大學(xué)學(xué)報(bào)2022,8,22(4):483486[16]石必明,俞啟香, 2022(33)3:3942[17][D].中國(guó)科學(xué)技術(shù)大學(xué),2022,11:8892[18]蘇現(xiàn)波,陳江峰，[M].科學(xué)出版社,2022,2[19]張新民,莊軍,[M].科學(xué)出版社,[20]Mavor M J,Pratt T J and C R Nelson. Quantitative Evaluation of Coal Seam Gas Content Estimate Accuracy[A].Paper SPE 29577 Presented at the 1995 SPE Joint Rocky Mountain Regional/LowPermeability Reservoirs Symp[C].Denver,1995:2022[21]王建博，王紅巖，[M].地質(zhì)出版社，[22]Wilfred Rupert Bodden. Coalbed methane Production and Desorption Testing in the Black Warrior Basin of Alabama[D].geology University of Alabama,1976[23][J].重慶大學(xué)學(xué)報(bào),2022,23(sl):15[24]梁冰,[M].科學(xué)出版社，[25][D].遼寧工程技術(shù)大學(xué)博士學(xué)位論文。2022[26][D].東北大學(xué)碩士論文 2022[27]陳炎光,[M].徐州:中國(guó)礦業(yè)大學(xué)出版社,1994[28]何國(guó)清,[M].中國(guó)礦業(yè)大學(xué)出版社,1991[29]“砌體梁”結(jié)構(gòu)的關(guān)鍵塊分析[J]煤炭學(xué)報(bào),1994(6)[30][D]中國(guó)礦業(yè)大學(xué)博士后研究工作報(bào)告 1990,4[31][M]煤炭工業(yè)出版社,1981（論文）37附錄 A例子描繪了流體在流動(dòng)域中的傳遞過程：流體在多孔介質(zhì)中流動(dòng)，由低速開始流動(dòng)，然后速度加快，最終由井筒套管進(jìn)入井內(nèi)。Darcy 定律描述了距離井口較遠(yuǎn)的流體的流動(dòng)速度；緊接著通過由 Brinkman 方程控制的流域；NavierStokes 方程描述了流體在井中的運(yùn)動(dòng)狀態(tài)。開始時(shí)流體的流動(dòng)模型是 Darcy 定律與 Brinkman 方程耦合來共同控制的，然后在井口是由 NavierStokes 方程控制的，所以流體在基質(zhì)的流動(dòng)有三個(gè)不同的方程來控制。：深部煤層群采動(dòng)影響下瓦斯運(yùn)移規(guī)律的數(shù)值模擬研究38welnsp?inletsu?brnss dlbru?brdlp?pumdl?Darcy0 m m 1 m 4 m圖 幾何邊界條件是用于 COMSOL Multiphysics 中建立Darcy、Brinkman 和 NavierStokes 耦合模型接下來的章節(jié)將會(huì)描述這個(gè)模型。第一節(jié)介紹 Darcy 定律、Brinkman 定律、NavierStokes 方程和所用的邊界條件。第二節(jié)利用表格說明各個(gè)步驟應(yīng)用到的數(shù)據(jù)。最后一節(jié)顯示結(jié)果，然后給出了用 COMSOL Multiphysics 建立模型的各個(gè)步驟的詳細(xì)解釋。模型的描述圖 敘述了幾何和邊界條件的模型，流體從以 r=4m 的圓柱區(qū)域內(nèi)開始，最終流入井中。流體在模型中的流動(dòng)是在一個(gè)可滲透的儲(chǔ)水層中實(shí)現(xiàn)的。假設(shè)流體的流動(dòng)在的區(qū)域?yàn)檫_(dá)西流動(dòng)，在 的區(qū)域?yàn)?Brinkman 流動(dòng)。在井中 ?r .??r ?r的流動(dòng)為 Brinkman 流動(dòng)。所以達(dá)西流動(dòng)和 Brinkman 流動(dòng)的交界面在 米處，?rBrinkman 流動(dòng)和 Navierstokes 流動(dòng)的交界面在 米處井壁套管的射流處，在套管.?的其他地方是沒有流體流動(dòng)的。設(shè)油層中的流體是均勻的和各向同性的，密度和粘度都是常數(shù)。假設(shè)流動(dòng)是穩(wěn)態(tài)的，油井出口處的流量和壓力都是已知的。這個(gè)問題是個(gè)二維軸對(duì)稱問題，用柱坐標(biāo) 表示，坐標(biāo)原點(diǎn)在井眼中心和上覆層??zr,交界處，如圖 1 所示，設(shè)井眼半徑為 米，滲透層的厚度為 米，?rw 1?b處的半徑為 米。假設(shè)滲透層和油藏相鄰，所以油藏的內(nèi)半徑就是 米。?res 孔在油層套管的中部，高度為 米。上下覆層很薄，且是均勻，各向同性的，具有不滲透性。BrinkmanNavierstokes（論文）39在這個(gè)例子中，用 COMSOL Multiphysics 解決流體在射孔的流動(dòng)壓力來代替在井頂端層面的固定的壓力值。而且知道抽水率和油井出路的壓力值，就可以解決這種穩(wěn)態(tài)流場(chǎng)的流動(dòng)問題。Darcy 定律達(dá)西定律描述了多孔介質(zhì)中由壓力梯度引起的流體流動(dòng)，因?yàn)榱鲌?chǎng)的模型是非常薄的，所以可以不考慮重力勢(shì)能的影響， 對(duì)于各向同性介質(zhì)中的穩(wěn)態(tài)流動(dòng)問題，由達(dá)西定律并應(yīng)用連續(xù)方程可得到控制方程： sdlQp????????????式中： —滲透率張量， ； —?jiǎng)诱扯认禂?shù)， ；?2m1?kgm—流體的密度， ； —重力加速度；f?3?kg—表示單位時(shí)間內(nèi)單位體積的油藏源產(chǎn)出的流體體積， ；sQ 1?s—在達(dá)西流動(dòng)區(qū)域中流體的壓力。dlp由于達(dá)西定律和布林克曼方程的交界面處應(yīng)滿足壓力和速度的連續(xù)性。因?yàn)樵谶_(dá)西定律中壓力是可變的，所以在達(dá)西布林克曼界面壓力的邊界條件是受到制約的，所以要滿足壓力連續(xù)條件： brdlp?式中：下標(biāo) 和 分別代表達(dá)西定律和布林克曼方程。dlbr其余的邊界條件包括，在不可滲透層對(duì)進(jìn)水口速度的限制。在進(jìn)水口，應(yīng)滿足適合的流動(dòng)速率。沒有流動(dòng)能夠通過圍層，邊界條件為達(dá)西定律問題： 在滲透層邊界處 brWpesdl???2??????????n 在 Darcy—Brinkman 交界處 （4）brdl 在上下覆層 0????????dlp??n式中： —滲透層邊界處的法線；下標(biāo) 和 分別代表達(dá)西定律和布林克曼方程。 dlbr：深部煤層群采動(dòng)影響下瓦斯運(yùn)移規(guī)律的數(shù)值模擬研究40Brinkman 方程Brinkman 方程描述了介于達(dá)西流動(dòng)和 Navierstokes 流動(dòng)之間的一種流動(dòng)情況。達(dá)西流動(dòng)描述的是多孔介質(zhì)中流體的流動(dòng)速度比較小的情形，這時(shí)由剪應(yīng)力引起的動(dòng)量傳遞比較小，所以達(dá)西流動(dòng)的控制方程中不考慮剪應(yīng)力效應(yīng)。當(dāng)流動(dòng)速度比較大，剪應(yīng)力效應(yīng)不能忽略時(shí)，就需要采用 Brinkman 方程來描述多孔介質(zhì)中的流動(dòng)：????0???????????????FpkbrrTbrbruu??0??br式中： —流體的密度， ； —?jiǎng)诱扯认禂?shù)，?3?kgm1?skgm—流動(dòng)速度矢量， ； —壓力， ；u1sp1?—滲透率， 。 k?k該布林克曼模式制約了速度在達(dá)西布林克曼界面：布林克曼速度等于達(dá)西速度對(duì)邊界。除為穿孔，這個(gè)例子定義了所有其他的界限，在流模型與無滑移條件，這是一份聲明中關(guān)于速度。對(duì)穿孔，該模型要求的條件，代表的連續(xù)性的壓力。邊界條件對(duì)布林克曼流區(qū)。 在達(dá)西流動(dòng)和 Brinkman 流動(dòng)的交界面上，速度是連續(xù)的，除了射孔以外，布林克曼速度應(yīng)等于達(dá)西速度，這個(gè)例子定義了所有其他的邊界條件，對(duì)于速度定義了在模型的流動(dòng)中無滑移的條件；對(duì)于射孔，應(yīng)滿足壓力連續(xù)的條件?？傊珺rinkman 流動(dòng)域的邊界條件為： 在 Darcy—Brinkman 交界處dlbru? 在上下覆層 0r 在套管處br 在射孔處??????0?????nuITbrbrrp?NavierStokes 方程井眼中流體