張先敏, 吳浩宇, 馮其紅, 劉 珊, 吳海明
(1.非常規(guī)油氣開發(fā)教育部重點(diǎn)實(shí)驗(yàn)室(中國石油大學(xué)(華東)),山東青島 266580; 2.中國石油華北油田公司勘探開發(fā)研究院,河北任丘 062552; 3.中國石油大學(xué)(華東)理學(xué)院,山東青島 266580; 4.河北華北石油天成實(shí)業(yè)集團(tuán)有限公司,河北任丘 062552)
中國主要含煤盆地普遍存在多煤層疊置現(xiàn)象,多煤層的發(fā)育使煤儲(chǔ)層總厚度較大、資源豐度較高,但由于受單煤層厚度制約,對(duì)單一煤儲(chǔ)層開發(fā)通常難以達(dá)到預(yù)期產(chǎn)氣效果,經(jīng)濟(jì)效益也較差,因此對(duì)多煤層發(fā)育地區(qū)實(shí)施多層合采工藝技術(shù)是降低煤層氣開發(fā)成本、提高產(chǎn)能的重要舉措之一[1]。近年來,諸多學(xué)者已就煤層氣藏多層合采的相關(guān)問題開展了大量理論研究[2-7]。自從1961年Lefkovits等[8]提出不考慮層間竄流的多層合采油藏不穩(wěn)定滲流模型以來,國內(nèi)外許多學(xué)者對(duì)該類油氣藏進(jìn)行了大量研究。Gao等[9]在不考慮吸附效應(yīng)的基礎(chǔ)上,構(gòu)建了多層氣藏壓力動(dòng)態(tài)分析的近似解析數(shù)學(xué)模型;之后,又進(jìn)一步將該模型推廣到了考慮吸附效應(yīng)的多層氣藏的壓力動(dòng)態(tài)分析[10]。李成勇等[11]建立了多層合采氣井的井底壓力動(dòng)態(tài)響應(yīng)數(shù)學(xué)模型,并分析了影響其動(dòng)態(tài)響應(yīng)特征的具體因素。劉啟國等[12]基于考慮井筒儲(chǔ)集和各層不同表皮情況下的層間無竄流多層氣藏滲流數(shù)學(xué)模型,分析了多層合采氣井分層產(chǎn)量貢獻(xiàn)情況。顧岱鴻等[13]開展了層間非均質(zhì)性致密氣藏多層合采產(chǎn)量的變化規(guī)律研究。以上多層合采氣井動(dòng)態(tài)分析相關(guān)研究多偏重于常規(guī)天然氣藏,然而有別于常規(guī)天然氣,煤層氣主要是以大分子團(tuán)吸附狀態(tài)賦存于煤儲(chǔ)層中,其獨(dú)特的吸附/解吸機(jī)制也決定了煤層氣運(yùn)移產(chǎn)出過程與常規(guī)天然氣截然不同,而且由于不同煤層的壓力系統(tǒng)與物性條件普遍存在差異,在多層合采方式下,煤層氣井的分層產(chǎn)氣動(dòng)態(tài)響應(yīng)特征也將更加復(fù)雜。為此,筆者在考慮井筒儲(chǔ)存和表皮效應(yīng)影響的情況下,建立封閉地層條件下多層合采煤層氣井流動(dòng)產(chǎn)出數(shù)學(xué)模型,探討多層合采方式下煤層氣井的分層產(chǎn)氣量動(dòng)態(tài)響應(yīng)變化規(guī)律,明確影響多層合采煤層氣井分層產(chǎn)氣貢獻(xiàn)的主控因素。
考慮n層均質(zhì)煤層中心一口氣井合采的滲流物理模型如圖1所示。各煤層之間具有良好的隔層,相互之間沒有流體交換,流體混合只發(fā)生在井筒當(dāng)中,并假設(shè):①煤層為均質(zhì)各向同性的連續(xù)雙重孔隙介質(zhì),巖石和流體微可壓縮;②多煤層疊置且各煤層的含氣性、厚度、孔隙度、滲透率、表皮系數(shù)、地層壓力和供氣半徑等參數(shù)存在差異,各層之間具有良好的非滲透隔層,層間流體的竄流僅發(fā)生在井筒內(nèi);③疊置各煤層中飽和甲烷氣體,含游離氣但不含可流動(dòng)水,各煤層中吸附氣解吸后,經(jīng)基質(zhì)擴(kuò)散、裂隙滲流進(jìn)入井筒內(nèi)合并采出,擴(kuò)散服從Fick擴(kuò)散定律,滲流滿足Darcy定律;④單相氣體流動(dòng)為等溫過程。
圖1 多層合采煤層氣井滲流物理模型Fig.1 Multi-layer coal seam gas well seepage physical model
對(duì)于第j煤層的裂隙系統(tǒng),由氣體運(yùn)動(dòng)方程、連續(xù)性方程以及真實(shí)氣體狀態(tài)方程,可以得到描述煤層氣體流動(dòng)的方程為
(1)
式中,j=1,2,…,n;pj為第j煤層的壓力,MPa;kj為第j煤層的裂隙滲透率,10-3μm2;μj為第j煤層的氣體黏度,mPa·s;Zj為第j煤層的氣體偏差因子;φj為第j煤層的裂隙孔隙度;psc為地面標(biāo)準(zhǔn)壓力,MPa;Tsc為地面標(biāo)準(zhǔn)溫度,K;Tj為第j煤層的溫度,K;FG為幾何相關(guān)系數(shù);Vj為第j煤層基質(zhì)內(nèi)氣體平均體積分?jǐn)?shù);t為時(shí)間,d。
對(duì)于第j煤層的基質(zhì)系統(tǒng),由Fick第一擴(kuò)散定律[14]可得:
(2)
其中
式中,τj為第j煤層的吸附時(shí)間,d;VE為第j煤層中基質(zhì)內(nèi)表面氣體體積分?jǐn)?shù),服從Langmuir等溫吸附方程;VLj為第j煤層的Langmuir體積分?jǐn)?shù);pLj為第j煤層的Langmuir壓力,MPa。
如果煤層之間的距離較小,考慮到井筒內(nèi)氣體流動(dòng)壓力損失也相對(duì)較小,可以認(rèn)為各煤層在井底處的流動(dòng)壓力相等,考慮井筒儲(chǔ)集及表皮效應(yīng)的影響,則內(nèi)邊界條件為
(3)
(4)
式中,qsc為地面標(biāo)準(zhǔn)狀況下的煤層氣井產(chǎn)量,m3/d;pw為井底流動(dòng)壓力,MPa;rw為井筒半徑,m;hj為第j煤層的有效厚度,m;Sj為第j煤層的表皮系數(shù),用以表征第j煤層附近井筒壁被污染或被改善的情況;C為井筒儲(chǔ)存系數(shù),m3/MPa。
初始條件為
pj|t=0=pIj.
(5)
外邊界條件為
(6)
第j煤層的分層產(chǎn)氣量為
(7)
式中,pIj為第j煤層的初始?jí)毫?MPa;qjsc為地面標(biāo)準(zhǔn)狀況下第j煤層的分層響應(yīng)產(chǎn)氣量,m3/d;rej為第j煤層的外邊界半徑,m。
煤層氣藏儲(chǔ)層壓力通常較低,可假設(shè)μZ為常數(shù)[15-16]。定義無因次變量如下:
式中,ψjD為第j煤層的無因次儲(chǔ)層壓力;ψwD為無因次井底流動(dòng)壓力;qjD為第j煤層的無因次響應(yīng)產(chǎn)氣量;ctj為第j煤層的綜合壓縮系數(shù),MPa-1;tD為無因次時(shí)間;ωj為第j煤層的彈性儲(chǔ)容比;λj為第j煤層的竄流系數(shù);αj為第j煤層的吸附系數(shù);βj為第j煤層的地層系數(shù)比;reDj為第j煤層的無因次外邊界半徑;VjD為第j煤層的無因次吸附量;VEDj為第j煤層的無因次平衡吸附量;CD為無因次井筒儲(chǔ)集系數(shù)。
采用上述無因次變量對(duì)方程(1)、(2)進(jìn)行無因次化處理,可得:
(8)
(9)
多層合采煤層氣藏?cái)M穩(wěn)態(tài)流動(dòng)問題的初邊值條件及第j煤層的無因次響應(yīng)產(chǎn)氣量[12,14]為
ψjD|tD=0=0;VjD|tD=0=0,
(10)
(11)
(12)
(13)
(14)
由方程(8)~(13)構(gòu)成封閉外邊界條件下n層均質(zhì)煤層中心一口氣井合采的無因次流動(dòng)數(shù)學(xué)模型,對(duì)上述方程進(jìn)行Laplace變換:
(15)
(16)
(17)
(18)
(19)
(20)
由方程(16)整理可得:
(21)
將式(21)代入式(15),整理可得:
(22)
其中
式(22)的無因次壓力通解為
(23)
其中
式中,I0(x)和K0(x)分別為第一類和第二類修正的0階貝塞爾函數(shù);Aj和Bj為待定系數(shù)。
令
F=K1(σjreDj)[I0(σj)-SjσjI1(σjrD)]+
I1(σjreDj)[K0(σj)+SjσjK1(σjrD)],
由邊界條件(18)及(19)可得Aj和Bj分別為
(24)
(25)
式中,I1(x)和K1(x)分別為第一類和第二類修正的1階貝塞爾函數(shù)。
由此可得:
(26)
再將式(26)代入邊界條件式(17),整理可得:
由式(26)和(20)可得到第j煤層的無因次儲(chǔ)層壓力及無因次響應(yīng)產(chǎn)氣量的Laplace空間解為
利用Stehfest方法[16]對(duì)無因次井底壓力以及各煤層的無因次儲(chǔ)層壓力、無因次響應(yīng)產(chǎn)氣量的Laplace空間解進(jìn)行數(shù)值反演,即可求出上述多層合采煤層氣藏流動(dòng)數(shù)學(xué)模型的近似數(shù)值解。
考慮封閉邊界條件下雙層疊置煤層氣藏中心一口合采氣井的情況,其中在地層系數(shù)比相同的條件下,上下煤層的儲(chǔ)層邊界、壓力系統(tǒng)、含氣性與儲(chǔ)層物性條件等方面存在差異。根據(jù)上述數(shù)學(xué)模型的無因次井底壓力與分層無因次響應(yīng)產(chǎn)氣量的數(shù)值解,繪制無因次井底壓力及其導(dǎo)數(shù)、分層無因次響應(yīng)產(chǎn)氣量與無因次時(shí)間的典型曲線圖版,如圖2所示。
從圖2中可以看出,雙層合采煤層氣井的無因次井底壓力及其導(dǎo)數(shù)的雙對(duì)數(shù)曲線可劃分為5個(gè)明顯的特征階段:階段Ⅰ為早期續(xù)流段,該階段主要受純井儲(chǔ)效應(yīng)的影響,無因次井底壓力與其導(dǎo)數(shù)曲線重合,且均為單位斜率直線段,隨后壓力導(dǎo)數(shù)曲線出現(xiàn)“駝峰”;階段Ⅱ?yàn)樵缙谶^渡流段,在該階段無因次井底壓力導(dǎo)數(shù)曲線上,反映裂縫徑向流特征的水平直線段已不復(fù)存在,壓力導(dǎo)數(shù)曲線出現(xiàn)下凹,出現(xiàn)反映各煤層中基質(zhì)向天然裂隙竄流過程的“V”型特征曲線段;階段Ⅲ為系統(tǒng)徑向流段,該階段發(fā)生在邊界響應(yīng)之前,各煤層的介質(zhì)間流體交換達(dá)到某種動(dòng)態(tài)平衡,無因次井底壓力導(dǎo)數(shù)曲線水平直線段對(duì)應(yīng)0.5的直線;階段Ⅳ為邊界差異響應(yīng)流段,該階段外邊界距離較小或物性條件較好的煤層壓力波傳遞到了外邊界后產(chǎn)生邊界效應(yīng),破壞了原有的單煤層流體交換動(dòng)態(tài)平衡關(guān)系,且由于合采各煤層的外邊界條件差異影響,導(dǎo)致無因次壓力導(dǎo)數(shù)曲線上呈現(xiàn)上翹,并出現(xiàn)邊界疊加響應(yīng)的“下凹”曲線段;階段Ⅴ為系統(tǒng)邊界線性流段,該階段各煤層的壓力波傳播均到達(dá)外邊界,顯現(xiàn)為單位斜率的線性流直線段。此外,從分層無因次響應(yīng)產(chǎn)氣量曲線可以看出:由于上下煤層在儲(chǔ)層邊界、壓力系統(tǒng)、含氣性與儲(chǔ)層物性等方面存在較大差異,生產(chǎn)過程中不可避免地出現(xiàn)了層間矛盾或干擾;在雙層合采方式下,上部煤層在一定程度上抑制了下部煤層的產(chǎn)氣釋放能力,因此在地層系數(shù)相等的條件下,上部煤層的分層產(chǎn)氣貢獻(xiàn)明顯高于下部煤層,分層產(chǎn)氣貢獻(xiàn)比例明顯不等于分層的地層系數(shù)比,尤其當(dāng)下部煤層壓降傳遞到外邊界后,其無因次響應(yīng)產(chǎn)氣量快速下降,而對(duì)應(yīng)的上部煤層產(chǎn)氣量則快速上升,最終上下煤層的分層產(chǎn)氣貢獻(xiàn)均維持在恒定水平,與郭平等[17]通過物理模擬實(shí)驗(yàn)與數(shù)值模擬計(jì)算得到的多層合采氣藏各小層產(chǎn)出分配規(guī)律認(rèn)識(shí)基本一致,該階段分層產(chǎn)氣貢獻(xiàn)與各層的地質(zhì)儲(chǔ)量比具有一定的關(guān)聯(lián)性。
圖2 無因次壓力及其導(dǎo)數(shù)和分層產(chǎn)氣量曲線圖版Fig.2 Type curves of dimensionless pressure and its derivative curves & response gas production
對(duì)于含煤盆地內(nèi)疊置的多煤層,由于受煤層的原生及后生變化影響,可使煤層出現(xiàn)尖滅、分叉、變薄和切斷等現(xiàn)象,從而造成疊置各煤層的邊界條件不同。圖3為依賴于下部煤層外邊界半徑的無因次井底壓力及其導(dǎo)數(shù)、分層無因次響應(yīng)產(chǎn)氣量的變化曲線。由圖3(a)可以看出,在無因次井底壓力及其導(dǎo)數(shù)特征曲線上,外邊界半徑reD2主要影響過渡流段以后的各流動(dòng)階段,在其他參數(shù)不變的條件下,reD2越大,下部煤層的邊界效應(yīng)出現(xiàn)時(shí)間越晚,系統(tǒng)徑向流段持續(xù)時(shí)間就越長,而且由于上下煤層之間的邊界差異越小,導(dǎo)致邊界差異響應(yīng)流段在壓力導(dǎo)數(shù)特征曲線上越加不明顯,乃至消失。從圖3(b)可以看出,在下部煤層壓力尚未傳遞至外邊界之前,各煤層的無因次響應(yīng)產(chǎn)氣貢獻(xiàn)基本保持不變。對(duì)于邊界差異響應(yīng)流段之后的分層產(chǎn)氣貢獻(xiàn)比例,reD2的影響顯著,隨著reD2的增大,下部煤層的無因次響應(yīng)產(chǎn)氣量逐漸增加,而上部煤層無因次響應(yīng)產(chǎn)氣量則隨之降低。
圖3 不同reD2下無因次壓力動(dòng)態(tài)變化和無因次響應(yīng)產(chǎn)氣量Fig.3 Dimensionless pressure dynamic and dimensionless gas production at different reD2
地層系數(shù)比為各煤層的地層系數(shù)與合采煤層系統(tǒng)的總地層系數(shù)之比,不同地層系數(shù)比對(duì)無因次井底壓力及響應(yīng)產(chǎn)氣量的影響如圖4所示。對(duì)于雙煤層合采系統(tǒng),在無因次井底壓力及其導(dǎo)數(shù)曲線上,對(duì)于邊界差異響應(yīng)流段出現(xiàn)之前的井底壓力動(dòng)態(tài),地層系數(shù)比的影響很小,而分層無因次響應(yīng)產(chǎn)氣量受地層系數(shù)比的影響較為明顯,早期續(xù)流段過后各煤層的分層產(chǎn)氣貢獻(xiàn)與地層系數(shù)比近似成正比。然而在下部煤層壓降傳遞到邊界之后,地層系數(shù)比對(duì)無因次井底壓力及響應(yīng)產(chǎn)氣量的影響顯著,其中下部煤層的地層系數(shù)β2值越大,在無因次井底壓力及其導(dǎo)數(shù)特征曲線上,邊界差異響應(yīng)流段出現(xiàn)時(shí)間越早且持續(xù)時(shí)間越長。該階段下部煤層的無因次響應(yīng)產(chǎn)氣量逐漸減小并趨于穩(wěn)定,與之同步,上部煤層的無因次響應(yīng)產(chǎn)氣量則逐漸增加,最后也趨于穩(wěn)定。
圖4 不同β2下無因次壓力動(dòng)態(tài)變化和無因次響應(yīng)產(chǎn)氣量Fig.4 Dimensionless pressure dynamic and dimensionless gas production at different β2
圖5為封閉邊界條件下吸附系數(shù)α2對(duì)無因次井底壓力及響應(yīng)產(chǎn)氣量的影響關(guān)系曲線。吸附系數(shù)主要用以表征各煤層中基質(zhì)微孔隙表面吸附氣體能力的程度[16],并決定著解吸氣向天然裂隙系統(tǒng)發(fā)生竄流的時(shí)刻和程度,即在其他參數(shù)一定的條件下,α2越小,其與上部煤層的吸附系數(shù)差異越明顯,在無因次井底壓力及其導(dǎo)數(shù)曲線圖上,早期過渡流段出現(xiàn)時(shí)間越晚且下凹程度越低,尤其當(dāng)α2=0.000 1時(shí),早期過渡流段出現(xiàn)了兩個(gè)反映雙煤層各自雙重孔隙介質(zhì)特征的“V”型曲線段,即各煤層中解吸氣體從基質(zhì)系統(tǒng)向天然裂隙系統(tǒng)竄流過程的同步性變差,且由于α2越小,代表維持多煤層合采系統(tǒng)流體交換平衡關(guān)系的吸附氣資源總量越小,系統(tǒng)徑向流段消失,邊界差異響應(yīng)流段出現(xiàn)時(shí)間越早,而且邊界差異響應(yīng)流段的“下凹”曲線段位置越低;從圖5(b)可以看出,α2對(duì)于早期續(xù)流段后的無因次響應(yīng)產(chǎn)氣量影響明顯,α2越大,則說明吸附氣資源量越大,在地層系數(shù)比一定的條件下,早期分層產(chǎn)氣貢獻(xiàn)比例越接近于地層系數(shù)比,尤其當(dāng)上下煤層的吸附系數(shù)相同時(shí),早期分層產(chǎn)氣貢獻(xiàn)比例即等于地層系數(shù)比,而當(dāng)下部煤層壓降傳遞到邊界之后,α2越大,系統(tǒng)邊界線性流階段的分層無因次響應(yīng)產(chǎn)氣量差異越小。
圖5 不同α2下無因次壓力動(dòng)態(tài)變化和無因次響應(yīng)產(chǎn)氣量Fig.5 Dimensionless pressure dynamic and dimensionless gas production at different α2
彈性儲(chǔ)容比用來描述各煤層中游離氣和吸附氣的彈性儲(chǔ)容能力相對(duì)大小。在其他參數(shù)一定的條件下,改變下部煤層的彈性儲(chǔ)容比ω2,并繪制了受其變化影響的無因次井底壓力及響應(yīng)產(chǎn)氣量的曲線變化關(guān)系,如圖6所示。在無因次井底壓力及其導(dǎo)數(shù)曲線圖上,彈性儲(chǔ)容比ω2的影響主要體現(xiàn)在早期過渡流段。在其他參數(shù)一定的條件下,下部煤層的彈性儲(chǔ)容比ω2越小,表示下部煤層中游離氣占有比例及資源量均越小,在合層開采時(shí)地層壓力下降越快。煤基質(zhì)中吸附氣體越早解吸進(jìn)入裂隙系統(tǒng),故無因次井底壓力導(dǎo)數(shù)曲線上的“V”型曲線段出現(xiàn)越早、下凹越深。從圖6(b)可以看出,下部煤層的彈性儲(chǔ)容比ω2對(duì)各煤層的早期無因次響應(yīng)產(chǎn)氣量具有一定的影響:ω2越大,下部煤層的早期無因次響應(yīng)產(chǎn)氣量越大,而對(duì)應(yīng)的上部煤層早期無因次響應(yīng)產(chǎn)氣量越小;當(dāng)上、下煤層中吸附氣體均大量解吸以后,ω2對(duì)分層產(chǎn)氣量的影響作用逐漸減小。
圖6 不同ω2下無因次壓力動(dòng)態(tài)變化和無因次響應(yīng)產(chǎn)氣量Fig.6 Dimensionless pressure dynamic and dimensionless gas production at different ω2
竄流系數(shù)是煤層多孔介質(zhì)的重要表征參數(shù)之一,用以表征煤基質(zhì)系統(tǒng)中吸附氣體解吸后向裂隙系統(tǒng)擴(kuò)散的難易程度。圖7為在保持其他參數(shù)不變、改變下部煤層的竄流系數(shù)λ2時(shí)的無因次井底壓力及響應(yīng)產(chǎn)氣量變化關(guān)系。從圖7(a)中可以看出,下部煤層的竄流系數(shù)λ2主要影響早期過渡流段與邊界響應(yīng)差異流段,其中λ2越大,表明下部煤層基質(zhì)孔隙中的吸附氣體越難以解吸擴(kuò)散進(jìn)入裂隙系統(tǒng)。因此無因次井底壓力及其導(dǎo)數(shù)曲線上的早期過渡流段出現(xiàn)時(shí)間越晚,且壓力導(dǎo)數(shù)曲線上的“V”型曲線段下凹越淺。尤其當(dāng)λ2=104時(shí),由于上下煤層中吸附氣體解吸擴(kuò)散難易程度差別明顯,造成相對(duì)于上部煤層,下煤層進(jìn)入早期過渡流段的時(shí)間出現(xiàn)明顯的滯后性,因此壓力導(dǎo)數(shù)曲線上早期過渡流段出現(xiàn)了兩個(gè)明顯的“V”型曲線段。如圖7(b)所示,λ2對(duì)初期和晚期的分層無因次響應(yīng)產(chǎn)氣量基本無影響,而對(duì)于早期續(xù)流段過后的分層產(chǎn)氣量影響較大,其中λ2越大,分層產(chǎn)氣貢獻(xiàn)比例越偏離分層地層系數(shù)比,即上下煤層之間的層間矛盾越突出,上部煤層對(duì)下部煤層產(chǎn)氣能力的抑制作用越明顯。值得一提的是,在地層系數(shù)比一定的條件下,當(dāng)λ2=102時(shí),早期上下煤層的分層產(chǎn)氣貢獻(xiàn)比例等于地層系數(shù)比。
圖7 不同λ2下無因次壓力動(dòng)態(tài)變化和無因次響應(yīng)產(chǎn)氣量Fig.7 Dimensionless pressure dynamic and dimensionless gas production at different λ2
綜合以上對(duì)吸附系數(shù)、彈性儲(chǔ)容比、竄流系數(shù)以及地層系數(shù)比的分析,可以得出在其他參數(shù)一定的條件下,多層合采系統(tǒng)中各煤層的吸附系數(shù)、彈性儲(chǔ)容比與竄流系數(shù)是影響分層產(chǎn)氣貢獻(xiàn)比例分配的主控因素。
(1)考慮井筒儲(chǔ)存與表皮效應(yīng)的影響,構(gòu)建了封閉邊界條件下多層合采煤層氣藏滲流產(chǎn)出數(shù)學(xué)模型,繪制了雙層合采煤層氣井的無因次井底壓力動(dòng)態(tài)響應(yīng)和分層產(chǎn)氣量曲線圖版,其中典型曲線圖版具有早期續(xù)流段、早期過渡流段、系統(tǒng)徑向流段、邊界差異響應(yīng)流段以及系統(tǒng)邊界線性流段5個(gè)明顯的特征階段。由于多煤層合采氣井的層間矛盾問題,導(dǎo)致分層產(chǎn)氣貢獻(xiàn)比例大多偏離于其地層系數(shù)比,晚期分層產(chǎn)氣貢獻(xiàn)比例與各層的地質(zhì)儲(chǔ)量比基本一致。
(2)吸附系數(shù)與竄流系數(shù)對(duì)初期和晚期的分層無因次響應(yīng)產(chǎn)氣量基本無影響,彈性儲(chǔ)容比主要影響早期過渡流段,而外邊界半徑、地層系數(shù)比則影響邊界差異響應(yīng)流段與系統(tǒng)邊界線性流段,而且各煤層的吸附系數(shù)、彈性儲(chǔ)容比與竄流系數(shù)的差異是影響多層合采煤層氣井分層產(chǎn)氣貢獻(xiàn)比例分配的主控因素。