基于陣列聲波測(cè)井的海陸過渡相碎屑巖地層裂縫識(shí)別方法

尹 帥， 丁文龍， 趙 威， 孫圓輝， 袁江如， 叢 森

(1.中國地質(zhì)大學(xué)(北京)能源學(xué)院，北京 100083；2.海相儲(chǔ)層演化與油氣富集機(jī)理教育部重點(diǎn)實(shí)驗(yàn)室(中國地質(zhì)大學(xué)(北京))，北京 100083；3.中國石油勘探開發(fā)研究院，北京 100083)

基于陣列聲波測(cè)井的海陸過渡相碎屑巖地層裂縫識(shí)別方法

尹 帥1，2， 丁文龍1，2， 趙 威1，2， 孫圓輝3， 袁江如3， 叢 森1，2

(1.中國地質(zhì)大學(xué)(北京)能源學(xué)院，北京 100083；2.海相儲(chǔ)層演化與油氣富集機(jī)理教育部重點(diǎn)實(shí)驗(yàn)室(中國地質(zhì)大學(xué)(北京))，北京 100083；3.中國石油勘探開發(fā)研究院，北京 100083)

針對(duì)海陸過渡相致密碎屑巖地層裂縫測(cè)井識(shí)別難度大的問題，提出結(jié)合Gassmann方程、自適應(yīng)基質(zhì)礦物及骨架模量提取方法及DEM理論模型，利用陣列聲波測(cè)井資料對(duì)海陸過渡相碎屑巖地層單井裂縫發(fā)育段進(jìn)行定量識(shí)別的方法。研究結(jié)果顯示，提取的碎屑巖地層巖石基質(zhì)礦物體積模量Ko分布在13～58 GPa，巖石體積模量Ks和干巖石骨架體積模量Kd略小于Ko；地層巖石基質(zhì)礦物剪切模量μo分布在5～18 GPa，巖石剪切模量μs略小于μo；利用自適應(yīng)方法可對(duì)地層橫波時(shí)差進(jìn)行有效預(yù)測(cè)，相對(duì)誤差為3.5%；將采用自適應(yīng)方法提取的模量參數(shù)代入DEM理論模型，可實(shí)現(xiàn)對(duì)地層裂縫的預(yù)測(cè)，該方法的識(shí)別效果明顯好于常規(guī)測(cè)井、裂縫參數(shù)法、多因素概率判別法及R/S法等常規(guī)裂縫識(shí)別方法。研究表明，提出的新識(shí)別方法其計(jì)算結(jié)果與實(shí)際值的吻合度高，裂縫識(shí)別結(jié)果可靠。

陣列聲波測(cè)井 海陸過渡相 碎屑巖 裂縫識(shí)別

我國擁有豐富的致密砂巖氣資源，其中鄂爾多斯盆地和四川盆地的技術(shù)可采資源量較大[1]。鄂爾多斯盆地海陸過渡相致密砂巖氣儲(chǔ)層油氣勘探實(shí)踐表明，致密砂巖氣藏雖然分布面積廣，但烴類聚集程度極為不均[2]，造成烴類富集出現(xiàn)差異的主導(dǎo)因素為裂縫發(fā)育程度[3]。在裂縫發(fā)育區(qū)的氣井，采取一定壓裂增產(chǎn)措施之后通常可獲得高產(chǎn)[3-6]。因此，對(duì)致密碎屑巖儲(chǔ)層裂縫發(fā)育層段及發(fā)育區(qū)準(zhǔn)確預(yù)測(cè)，會(huì)直接影響油氣勘探有利區(qū)的劃定及后期儲(chǔ)層改造措施的有效制定[3-4]。

利用測(cè)井方法識(shí)別單井裂縫發(fā)育段可以為壓裂選層提供指導(dǎo)性建議[7-9]，并且該類方法的識(shí)別精度要高于地震方法[10-12]。測(cè)井識(shí)別地層裂縫的方法包括常規(guī)測(cè)井(如聲波、中子、密度及電阻率等)、特殊測(cè)井(如FMI[13-14])及數(shù)值轉(zhuǎn)換方法(如神經(jīng)網(wǎng)絡(luò)、概率判別等[15])等，針對(duì)不同地區(qū)儲(chǔ)層，不同方法的應(yīng)用效果具有一定差異。為此，筆者以沁南樊莊區(qū)塊北部的固縣井區(qū)山西組海陸過渡相碎屑巖地層為例，提出采用基于陣列聲波測(cè)井的裂縫識(shí)別方法，并與其他常規(guī)裂縫識(shí)別方法進(jìn)行了對(duì)比。

1 山西組碎屑巖地層裂縫特征

沁南樊莊區(qū)塊北部的固縣井區(qū)，目前以開采山西組煤層氣為主[16-17]，地層自下而上包括奧陶系、石炭系(本溪組及太原組)、二疊系(山西組、下石盒子組及上石盒子組)、三疊系及第四系。其中，山西組地層埋深300～800 m，巖性以細(xì)粒砂巖、粉砂巖、泥巖互層為主，含煤4～5層。

山西組地層屬海陸過渡相沉積，經(jīng)歷了多期構(gòu)造運(yùn)動(dòng)[18]，區(qū)內(nèi)北東及北西向褶曲分別與燕山期北西向擠壓及喜山期北東向擠壓相關(guān)。沁南晉試1井地層埋藏史(見圖1)[19]表明，在三疊紀(jì)末期，地層埋深達(dá)到最大，鏡質(zhì)體反射率Ro達(dá)到1.3%，深埋藏強(qiáng)烈壓實(shí)作用是導(dǎo)致地層巖石致密化的重要原因。同時(shí)，研究區(qū)地層發(fā)育疊合不整合面，多期次構(gòu)造運(yùn)動(dòng)、構(gòu)造疊加、強(qiáng)烈壓實(shí)作用及復(fù)雜成巖作用等因素共同作用于山西組地層，導(dǎo)致裂縫發(fā)育特征具有復(fù)雜性和特殊性。

對(duì)研究區(qū)內(nèi)11口取心井山西組碎屑巖地層進(jìn)行裂縫觀察統(tǒng)計(jì)發(fā)現(xiàn)，地層中裂縫類型以構(gòu)造裂縫為主，同時(shí)發(fā)育多種非構(gòu)造裂縫。構(gòu)造裂縫包括張性裂縫、剪切裂縫及擠壓裂縫，這些裂縫往往發(fā)生了一定程度的石英及方解石充填；非構(gòu)造裂縫包括層間頁理縫、溶蝕縫、成巖收縮縫及滑脫縫，其中滑脫縫多見于泥巖[20]。此外，碎屑巖地層中可見大量炭質(zhì)條帶，多呈細(xì)(或粗)條帶狀、斷續(xù)、切割條紋、網(wǎng)狀及斑點(diǎn)狀分布，常具有一定組系，多為近水平或低角度產(chǎn)狀。

研究區(qū)3口井山西組碎屑巖地層FMI成像測(cè)井統(tǒng)計(jì)獲得的地層條件下的裂縫傾角及充填特征結(jié)果如圖2所示。

從圖2可以看出，研究區(qū)目的層中以小于60°的中-低角度縫為主；地層條件下裂縫以低阻未充填及半充填裂縫為主，這些裂縫均為有效裂縫，可以為油氣的儲(chǔ)集和運(yùn)移提供場(chǎng)所或通道[21]；地層中同時(shí)發(fā)育少量高阻全充填縫。

2 基于陣列聲波測(cè)井的裂縫發(fā)育段識(shí)別方法

2.1 利用Gassmann方程計(jì)算相關(guān)參數(shù)

利用Gassmann方程[22]可以對(duì)地層條件下的巖石體積模量和剪切模量進(jìn)行計(jì)算，進(jìn)而根據(jù)巖石體積模量和剪切模量獲得巖石縱波、橫波速度。上述4參數(shù)可表示為：

(1)

μs=μd

(2)

(3)

(4)

式中：Ks和μs分別為地層條件下的巖石體積模量和剪切模量，GPa；Kd和μd分別為干巖石骨架的體積模量和剪切模量，GPa；Ko為巖石基質(zhì)礦物體積模量，GPa；Kf為巖石中所飽和流體的體積模量，GPa；μd為巖石干巖石骨架剪切模量，GPa；vp和vs分別為縱波速度和橫波速度，km/s；φ為巖石孔隙度；ρs為地層條件下的巖石密度，g/cm3。

2.2 利用自適應(yīng)方法提取巖石基質(zhì)礦物及骨架模量

自適應(yīng)方法的思路為：在合理取值區(qū)間內(nèi)對(duì)不同巖性巖石的基質(zhì)礦物及骨架模量進(jìn)行賦值，利用Gassmann方程[22]對(duì)反算的地層巖石縱橫波速度與陣列聲波測(cè)井值進(jìn)行對(duì)比。通過反復(fù)校驗(yàn)，不斷調(diào)整各基質(zhì)礦物及骨架模量參數(shù)值至最優(yōu)，使所反算的縱橫波速度預(yù)測(cè)值與實(shí)測(cè)值吻合度最高，此時(shí)所確定的巖石基質(zhì)礦物及骨架模量即為最優(yōu)解。具體步驟詳述如下。

首先引入一個(gè)參數(shù)σd(干巖石泊松比)[23]，其取值為0～0.45；同時(shí)定義：

(5)

將Biot系數(shù)β引入到Gassman-Biot-Geertsma方程[23]中，可得到有關(guān)Biot系數(shù)的一元二次方程：

(6)

其中，中間變量M可表示為：

(7)

Biot系數(shù)β可根據(jù)陣列聲波測(cè)井?dāng)?shù)據(jù)用式(8)進(jìn)行計(jì)算和解釋[24]。對(duì)于山西組碎屑巖地層，ρma取2.65 g/cm3，Δtmac取182.09 μs/m，Δtmas取288.71 μs/m。β可表示為：

(8)

式中：ρ為巖石密度，g/cm3；ρma為巖石骨架密度，g/cm3；Δtc為巖石縱波時(shí)差，μs/m；Δts為巖石橫波時(shí)差，μs/m；Δtmac為巖石骨架縱波時(shí)差，μs/m；Δtmas為巖石骨架橫波時(shí)差，μs/m。

通過在合理取值區(qū)間內(nèi)不斷調(diào)整參數(shù)σd及Ko，利用式(6)及陣列聲波測(cè)井資料對(duì)干巖石骨架體積模量進(jìn)行求解，將求算結(jié)果代入式(1)，然后利用式(3)對(duì)地層巖石縱波速度進(jìn)行解釋。將地層巖石縱波速度解釋結(jié)果與陣列聲波測(cè)井實(shí)測(cè)值進(jìn)行對(duì)比，兩者差值最小時(shí)所確定的σd及Ko值即認(rèn)為是目標(biāo)函數(shù)的最優(yōu)解。由于μs=μd，可根據(jù)式(5)求取μs，進(jìn)而求得地層巖石基質(zhì)礦物剪切模量，其表達(dá)式為[23]：

(9)

式中：μo為地層巖石基質(zhì)礦物的剪切模量，GPa。

利用自適應(yīng)方法提取巖石基質(zhì)礦物及骨架模量后，將地層巖石縱波、橫波速度預(yù)測(cè)值與陣列聲波測(cè)井值進(jìn)行對(duì)比，當(dāng)兩者很吻合時(shí)，即認(rèn)為所定義各參數(shù)均為最優(yōu)解。

2.3 自適應(yīng)方法提取結(jié)果及驗(yàn)證

研究區(qū)HF-A井(垂直井)鉆穿山西組地層，解釋層段為山西組3號(hào)煤層頂部的海陸過渡相碎屑巖地層(491～514 m井段)，主要為細(xì)粒砂巖、粉砂巖及泥巖互層。含氣性測(cè)井顯示，該地區(qū)山西組地層(主要是致密砂巖段)氣測(cè)異常普遍，好的氣測(cè)顯示往往與裂縫發(fā)育程度較高相關(guān)。解釋井段具有可靠的常規(guī)測(cè)井、陣列聲波測(cè)井及成像測(cè)井資料。

利用上述自適應(yīng)方法提取了山西組碎屑巖地層巖石各模量，利用式(3)和式(4)解釋的巖石縱橫波速度與陣列聲波測(cè)井值間關(guān)系見圖3。從圖3可以看出，地層巖石縱波、橫波速度預(yù)測(cè)值與實(shí)測(cè)值很吻合，表明所提取的地層各模量值可靠。

所提取的地層巖石各模量包括Ko，Kd,Ks，μs及μo(見圖4)。從圖4可以看出，對(duì)于該海陸過渡相碎屑巖地層，Ko>Ks>Kd，其中Ko分布在13～58 GPa；Ks和Kd略小于Ko；μo分布在5～18 GPa，μs略小于μo。

為進(jìn)一步說明所提取碎屑巖地層各模量的可靠性，將所提取的Ko及μo值與通過礦物組分法[25]大致估算的Ko及μo值進(jìn)行對(duì)比。礦物組分法認(rèn)為不同巖性巖石的Ko及μo可以通過對(duì)礦物組分含量及各單礦物組分的Ko及μo值進(jìn)行大致估算。由于Ko和μo的估算方法類似，此處僅列出Ko的估算方程[25]：

(10)

式中：N為巖石中基質(zhì)礦物的種類；Ki為第i種基質(zhì)礦物的體積模量，GPa；fi為第i種礦物的體積分?jǐn)?shù)，%。

從式(10)可以看出，巖石基質(zhì)礦物模量與單種礦物模量值、所含礦物種類及礦物體積分?jǐn)?shù)有關(guān)。研究區(qū)海陸過渡相碎屑巖地層中的致密砂巖以石英為主，此外還含有部分長(zhǎng)石(鉀長(zhǎng)石及斜長(zhǎng)石)、方解石、黃鐵礦、白云石及黏土；而泥巖則以黏土礦物為主，此外還含有石英、長(zhǎng)石、方解石及黃鐵礦，各單礦物的Ko及μo值見表1[24]。

表1 基質(zhì)礦物彈性模量

Table 1 Elastic modulus of matrix minerals GPa

用礦物組分法分析可知，山西組碎屑巖地層的Ko主要分布在20～55 GPa，而μo則主要分布在8～20 GPa。該估值區(qū)間與前述通過自適應(yīng)方法提取的Ko及μo區(qū)間相符，這進(jìn)一步表明采用自適應(yīng)方法提取的山西組地層各巖石力學(xué)參數(shù)可靠。

圖5為幾組地層模量與波阻抗之間的關(guān)系。

從圖5可以看出，所提取的地層巖石各模量參數(shù)之間或模量參數(shù)與波阻抗之間均具有非常好的二項(xiàng)式關(guān)系，相關(guān)系數(shù)均非常高。圖5(c)中μs與波阻抗之間的轉(zhuǎn)換關(guān)系可用于研究區(qū)橫波時(shí)差未知井的相應(yīng)預(yù)測(cè),該方法對(duì)地層巖石橫波時(shí)差預(yù)測(cè)的相對(duì)誤差較小(約為3.5%)。分析表明，自適應(yīng)方法不僅提取了可靠的碎屑巖目地層各模量參數(shù)，而且可對(duì)目的層巖石橫波時(shí)差進(jìn)行高精度預(yù)測(cè)。

2.4 基于DEM理論模型的裂縫識(shí)別方法

對(duì)目的層巖石基質(zhì)礦物模量及干巖石骨架模量進(jìn)行準(zhǔn)確提取后，利用DEM(微分等效介質(zhì))理論模型對(duì)碎屑巖地層不同巖性巖石孔隙縱橫比α進(jìn)行反演[26]。前人定義α為巖石中橢球狀孔隙的短軸與長(zhǎng)軸之比，其取值為0～1[27]。根據(jù)不同的α取值，地層巖石孔隙形態(tài)可分為4類[26]：1)α>0.90，球形孔；2)0.10<α<0.90，針形孔；3)0.01<α<0.10，硬幣形裂縫；4)α<0.01，無限型狹長(zhǎng)裂縫，開度往往較小。很顯然，可將α=0.10作為地層巖石中主要發(fā)育孔隙或裂縫的分界：當(dāng)α>0.10時(shí)，表示地層中主要發(fā)育孔隙；而當(dāng)α<0.10時(shí)，表示地層中發(fā)育裂縫，為裂縫性地層。

α的反演依據(jù)為，孔裂隙的存在會(huì)在巖石彈性波速上產(chǎn)生顯著響應(yīng)，前人對(duì)此多有研究[28-29]。J.G.Berryman等人[30]通過大量試驗(yàn)，最早給出了基于DEM理論的地層巖石體積及剪切模量的耦合常微分表達(dá)方程：

(11)

(12)

式中：y為巖石孔隙度；Ki為孔隙內(nèi)組分的體積模量，GPa；μi為孔隙內(nèi)組分的剪切模量，GPa；K*和μ*分別為考慮基質(zhì)礦物及孔隙組分的巖石等效體積及剪切模量，GPa；P*i和Q*i分別為K*和μ*的激化因子。

通過對(duì)式(11)、式(12)進(jìn)行解耦，可得到干巖石骨架體積模量與剪切模量之比[31]：

(13)

式(13)中，a和b均為α的函數(shù)，分別為下式的梯度和截距[30-31]：

(14)

采用線性方法對(duì)式(14)中P*i-Q*i進(jìn)行一階求導(dǎo)，可確定b和a值[25,27]：

b=(P*i-Q*i)′

(15)

(16)

式(15)、式(16)中，P*i-Q*i和(P*i-Q*i)′均為α的函數(shù)，可通過A,B，R，θ和g等幾個(gè)中間過程參數(shù)進(jìn)行表達(dá)，具體表達(dá)式見文獻(xiàn)[26]。首先利用α對(duì)a和b進(jìn)行定量表征，然后根據(jù)式(13)對(duì)Kd/μd進(jìn)行測(cè)井解釋。將所解釋的Kd/μd值與陣列聲波測(cè)井值進(jìn)行對(duì)比，確定滿足最小誤差前提下的α值，從而完成地層巖石孔隙縱橫比反演。

3 裂縫識(shí)別結(jié)果及對(duì)比

3.1 基于陣列聲波測(cè)井的裂縫發(fā)育段識(shí)別

利用式(13)所反演的HF-A井山西組海陸過渡相碎屑巖地層各類巖石的Kd/μd值及與前述利用自適應(yīng)方法提取的Kd/μd的對(duì)比結(jié)果見圖6。

從圖6可以看出，目的層致密砂巖的Kd/μd相對(duì)高一些，往往大于3.0；而泥巖Kd/μd相對(duì)低一些，但分布范圍比較廣，為1.0～3.3。2種方法所得Kd/μd很吻合，分布在中值線上，表明反演的巖石孔隙縱橫比α結(jié)果可靠。

3.2 其他裂縫識(shí)別方法

為了考察基于陣列聲波測(cè)井的裂縫識(shí)別方法的效果，將其與測(cè)井裂縫參數(shù)計(jì)算方法、R/S變尺度方法[32]及多因素裂縫判別方法[33]對(duì)HF-A井目的層的裂縫識(shí)別結(jié)果進(jìn)行了對(duì)比。測(cè)井裂縫參數(shù)主要計(jì)算了裂縫孔隙度、裂縫開度及裂縫滲透率等3個(gè)參數(shù)。

對(duì)于R/S變尺度方法，其識(shí)別裂縫的依據(jù)為裂縫的存在能增加測(cè)井?dāng)?shù)據(jù)的復(fù)雜性[32]。通過研究單測(cè)井參數(shù)序列的過程序列全層段極差R(N)與標(biāo)準(zhǔn)差S(N)的比值與測(cè)井點(diǎn)數(shù)(N)之間的變化關(guān)系，對(duì)天然裂縫發(fā)育段進(jìn)行預(yù)測(cè)。研究發(fā)現(xiàn)，對(duì)于該地區(qū)海陸過渡相碎屑巖地層而言，測(cè)井序列(筆者采用淺側(cè)向電阻率)的R/S與N在雙對(duì)數(shù)坐標(biāo)中時(shí)對(duì)裂縫的識(shí)別能力較差(見圖7)。而測(cè)井序列的R/S與N的線性交匯圖則能較好地識(shí)別出裂縫發(fā)育段(見圖8)，從圖8可以識(shí)別出目的層中的5個(gè)裂縫發(fā)育段，表現(xiàn)為斜率突然發(fā)生降低。雖然圖7中R/S與N在雙對(duì)數(shù)坐標(biāo)中時(shí)不能有效識(shí)別出地層裂縫發(fā)育段，但對(duì)其整體進(jìn)行線性擬合，可確定單井目的層整體斜率，該斜率值可用于不同單井裂縫發(fā)育程度的對(duì)比，其值越低，表示單井裂縫發(fā)育程度相對(duì)更高一些。

對(duì)于多因素裂縫概率判別[33]，通過對(duì)目的層裂縫(及不含裂縫)測(cè)井參數(shù)進(jìn)行提取(50組，裂縫及非裂縫參數(shù)各25組)，根據(jù)貝葉斯判別分析準(zhǔn)則[33]逐步建立裂縫判別方程Y?？紤]到實(shí)際各單井所含測(cè)井資料，最終選用縱波時(shí)差、井徑、體積密度、深側(cè)向電阻率及淺側(cè)向電阻率對(duì)目的層進(jìn)行多因素裂縫概率判別。當(dāng)判別函數(shù)Y大于1.5時(shí)，為裂縫發(fā)育地層段；當(dāng)Y小于1.5時(shí)，地層裂縫不發(fā)育。50組樣本參數(shù)的判別準(zhǔn)確率約為80%。

3.3 裂縫識(shí)別方法及結(jié)果對(duì)比

綜合常規(guī)測(cè)井、裂縫參數(shù)、多因素概率判別、R/S法裂縫識(shí)別及基于陣列聲波測(cè)井裂縫識(shí)別方法的HF-A井山西組目的層裂縫測(cè)井綜合識(shí)別成果見圖9。該井段具有經(jīng)過巖性歸位校正后的FMI成像測(cè)井，其中裂縫綜合識(shí)別結(jié)果為FMI成像測(cè)井中所識(shí)別出的裂縫發(fā)育段。

從圖9可以看出，該井段山西組地層發(fā)育5個(gè)裂縫發(fā)育段。這些裂縫發(fā)育段特征在部分常規(guī)測(cè)井(如井徑、密度)及裂縫參數(shù)(如裂縫滲透率)上的裂縫響應(yīng)并不是特別顯著，這與該地區(qū)巖石致密化所導(dǎo)致的裂縫特征較為復(fù)雜有關(guān)，導(dǎo)致基于常規(guī)測(cè)井?dāng)?shù)據(jù)的多因素裂縫概率判別方法的裂縫識(shí)別效果一般。對(duì)比裂縫綜合識(shí)別結(jié)果，多因素裂縫概率判別方法的吻合度僅為67%。對(duì)于R/S法，該方法能有效識(shí)別出碎屑巖地層中的5個(gè)裂縫發(fā)育段。但在單裂縫發(fā)育段的縱向識(shí)別精度上具有一定出入，該方法所解釋的裂縫發(fā)育段約為80%，相比前述多因素判別法有顯著提高。

基于陣列聲波測(cè)井的裂縫識(shí)別方法能從聲學(xué)角度對(duì)地層巖石孔縫形態(tài)進(jìn)行反演，圖9中α<0.1所指示的地層埋深段即為裂縫發(fā)育段。該方法不僅考慮了地層中常規(guī)縱波的變化，而且考慮了橫波的變化，這2個(gè)參數(shù)一般來說均能對(duì)地層中裂縫的有無及存在狀態(tài)產(chǎn)生顯著響應(yīng)。因此，從圖9的解釋結(jié)果來看，裂縫識(shí)別效果最好。除了497.1～498.0 m的裂縫發(fā)育段未能識(shí)別出來之外，其他4個(gè)裂縫發(fā)育段均能有效識(shí)別，且識(shí)別精度較高，吻合度約為92%。該方法的裂縫識(shí)別精度依賴于橫波時(shí)差的可靠性，對(duì)具有可靠陣列聲波測(cè)井的單井地層識(shí)別效果較好；當(dāng)沒有該項(xiàng)測(cè)井時(shí)，其裂縫識(shí)別精度取決于地層巖石橫波時(shí)差預(yù)測(cè)效果[34-35]。

上述研究結(jié)果表明，基于陣列聲波測(cè)井的裂縫識(shí)別方法適用于研究區(qū)海陸過渡相碎屑巖地層,相比其他常規(guī)識(shí)別方法具有明顯優(yōu)勢(shì)，因而具有一定應(yīng)用前景。

4 結(jié) 論

1) 提出了基于陣列聲波測(cè)井資料識(shí)別海陸過渡相碎屑巖地層裂縫的方法，通過聯(lián)合Gassmann方程、自適應(yīng)基質(zhì)礦物及骨架模量提取方法及DEM理論模型，可以實(shí)現(xiàn)對(duì)單井裂縫發(fā)育段的有效識(shí)別。

2) 所提取的地層巖石基質(zhì)礦物體積模量Ko分布在13～58 GPa，巖石體積模量Ks和干巖石骨架體積模量Kd略小于Ko；地層巖石基質(zhì)礦物剪切模量μo分布在5～18 GPa，巖石剪切模量μs略小于μo；利用自適應(yīng)方法同時(shí)可以對(duì)地層橫波時(shí)差進(jìn)行有效預(yù)測(cè)，相對(duì)誤差僅為3.5%。

3) 基于陣列聲波測(cè)井的裂縫識(shí)別方法其效果明顯好于常規(guī)測(cè)井、裂縫參數(shù)法、多因素概率判別法和R/S法等常規(guī)裂縫識(shí)別方法，表明采用該方法對(duì)海陸過渡相碎屑巖地層裂縫進(jìn)行識(shí)別，結(jié)果更為可靠。

References

[1] 李建忠，郭彬程，鄭民，等.中國致密砂巖氣主要類型、地質(zhì)特征與資源潛力[J].天然氣地球科學(xué)，2012，23(4)：607-615. Li Jianzhong,Guo Bincheng,Zheng Min,et al.Main types,geological features and resource potential of tight sandstone gas in China[J].Natural Gas Geoscience,2012,23(4):607-615.

[2] 童曉光，郭彬程，李建忠，等.中美致密砂巖氣成藏分布異同點(diǎn)比較研究與意義[J].中國工程科學(xué)，2012，14(6)：9-15，30. Tong Xiaoguang,Guo Bincheng,Li Jianzhong,et al.Comparison study on accumulation & distribution of tight sandstone gas between China and the United States and its significance[J].Engineering Science,2012,14(6):9-15,30.

[3] 丁文龍，尹帥，王興華，等.致密砂巖氣儲(chǔ)層裂縫評(píng)價(jià)方法與表征[J].地學(xué)前緣，2015，22(4)：173-187. Ding Wenlong,Yin Shuai,Wang Xinghua,et al.Assessment method and characterization of tight sandstone gas reservoir fractures[J].Earth Science Frontiers,2015,22(4):173-187.

[4] Olson J E,Laubach S E,Lander R H.Natural fracture characterization in tight gas sandstones:integrating mechanics and diagenesis[J].AAPG Bulletin,2009,93(11):1535-1549.

[5] Nelson R A.Geological analysis of naturally fractured reservoirs[M].Houston :Gulf Publishing Company,1985:8-26.

[6] 曾義金.頁巖氣開發(fā)的地質(zhì)與工程一體化技術(shù)[J].石油鉆探技術(shù)，2014，42(1)：1-6. Zeng Yijin.Integration technology of geology & engineering for shale gas development[J].Petroleum Drilling Techniques,2014,42(1):1-6.

[7] Zeng Lianbo,Li Xiangyang.Fractures in sandstone reservoirs with ultra-low permeability:a case study of the Upper Triassic Yanchang Formation in the Ordos Basin,China[J].AAPG Bulletin,2009,93(4):461-477.

[8] 趙立強(qiáng)，劉飛，王佩珊，等.復(fù)雜水力裂縫網(wǎng)絡(luò)延伸規(guī)律研究進(jìn)展[J].石油與天然氣地質(zhì)，2014，35(4)：562-569. Zhao Liqiang,Liu Fei,Wang Peishan,et al.A review of creation and propagation of complex hydraulic fracture network[J].Oil & Gas Geology,2014,35(4):562-569.

[9] 尹帥，單鈺銘，周文，等.破裂準(zhǔn)則方程在深層巖石力學(xué)強(qiáng)度響應(yīng)中的應(yīng)用[J].地球物理學(xué)進(jìn)展，2014，29(6)：2942-2949. Yin Shuai,Shan Yuming,Zhou Wen,et al.Application of failure criterion equation in deep rock’s mechanical strength response[J].Progress in Geophysics,2014,29(6):2942-2949.

[10] 孫煒，李玉鳳，付建偉，等.測(cè)井及地震裂縫識(shí)別研究進(jìn)展[J].地球物理學(xué)進(jìn)展，2014，29(3)：1231-1242. Sun Wei,Li Yufeng,Fu Jianwei,et al.Review of fracture identification with well logs and seismic data[J].Progress in Geophysics,2014,29(3):1231-1242.

[11] Murray G H.Quantitative fracture study-Sanish Pool McKenzie country,North Dakota[J].AAPG Bulletin,1968,52(1):57-65.

[12] 尹帥，單鈺銘，周文，等.Hoek-Brown準(zhǔn)則在致密砂巖彈性參數(shù)測(cè)井解釋中的應(yīng)用[J].石油鉆探技術(shù)，2015，43(1)：88-95. Yin Shuai,Shan Yuming,Zhou Wen,et al.Application of Hoek-Brown criterion for tight sandstone elastic parameters in log interpretation[J].Petroleum Drilling Techniques,2015,43(1):88-95.

[13] 丁文龍，漆立新，呂海濤，等.利用FMI資料分析塔河油田南部中-下奧陶統(tǒng)儲(chǔ)層構(gòu)造應(yīng)力場(chǎng)[J].現(xiàn)代地質(zhì)，2009， 23(5)：852-859. Ding Wenlong,Qi Lixin,Lv Haitao,et al.Analysis of the lower-middle ordovician reservoir tectonic stress field using FMI data in the south of Tahe Oilfield[J].Geoscience,2009,23(5):852-859.

[14] Ding Wenlong,Fan Tailiang,Yu Bingsong,et al.Ordovician carbonate reservoir fracture characteristics and fracture distribution forecasting in the Tazhong Area of Tarim Basin,Northwest China[J].Journal of Petroleum Science and Engineering,2012,86/87(5):62-70.

[15] Jahanbakhshi R,Keshavarzi R.Intelligent prediction of wellbore stability in oil and gas wells:an artificial neural network approach[J].American Rock Mechanics Association,2012,46:1-9.

[16] 秦勇，梁建設(shè)，申建，等.沁水盆地南部致密砂巖和頁巖的氣測(cè)顯示與氣藏類型[J].煤炭學(xué)報(bào)，2014，39(8)：1559-1565. Qin Yong,Liang Jianshe,Shen Jian,et al.Gas logging shows and gas reservoir types in tight sandstones and shales from Southern Qinshui Basin[J].Journal of China Coal Society,2014,39(8):1559-1565.

[17] 丁文龍，梅永貴，尹帥，等.沁水盆地煤系地層孔-裂隙特征測(cè)井反演[J].煤炭科學(xué)技術(shù)，2015，43(2)：53-57. Ding Wenlong,Mei Yonggui,Yin Shuai,et al.Logging inversion on pore-crack features of coal measure strata in Qinshui Basin[J].Coal Science and Technology,2015,43(2):53-57.

[18] 李月，林玉祥，于騰飛.沁水盆地構(gòu)造演化及其對(duì)游離氣藏的控制作用[J].桂林理工大學(xué)學(xué)報(bào)，2011，31(4)：481-487. Li Yue,Lin Yuxiang,Yu Tengfei.Tectonic evolution of Qinshui Basin and free gas reservoir control[J].Journal of Guilin University of Technology,2011,31(4):481-487.

[19] Su X B,Lin X Y,Zhao M J,et al.The upper Paleozoic coalbed methane system in the Qinshui Basin,China[J].AAPG Bulletin,2005,89(1):81-100.

[20] 丁文龍，許長(zhǎng)春，久凱，等.泥頁巖裂縫研究進(jìn)展[J].地球科學(xué)進(jìn)展，2011，26(2)：135-144. Ding Wenlong,Xu Changchun,Jiu Kai,et al.The research progress of shale fractures[J].Advances in Earth Science,2011,26(2):135-144.

[21] Bogdonov A A.The intensity of cleavage as related to the thickness of beds[J].Soviet Geology，1947，16:102-104.

[22] Robert G K,Xu Shiyu.An approximation for the Xu-White Velocity[J].Geophysics,2002,67(5):1406-1414.

[23] 林凱，賀振華，熊曉軍，等.基于基質(zhì)礦物模量自適應(yīng)提取橫波速度反演方法[J].石油地球物理勘探，2013，48(2)：262-267. Lin Kai,He Zhenhua,Xiong Xiaojun,et al.S-wave velocity inversion based on adaptive extraction of matrix mineral modulus[J].Oil Geophysical Prospecting,2013,48(2):262-267.

[24] 雍世和，張超謨.測(cè)井?dāng)?shù)據(jù)處理與綜合解釋[M].東營(yíng)：石油大學(xué)出版社，1996：235-240. Yong Shihe,Zhang Chaomo.Logging data processing and comprehensive interpretation[M].Dongying:Petroleum University Press,1996:235-240.

[25] He Xilei,He Zhenhua,Wang Ruiliang,et al.Calculations of rock matrix modulus based on a linear regression relation[J].Applied Geophysics,2011,8(3):155-162.

[26] Li Hongbing,Zhang Jiajia.Elastic moduli of dry rocks containing spheroidal pores based on differential effective medium theory[J].Journal of Applied Geophysics,2011,75(4):671-678.

[27] 董寧，霍志周，孫贊東，等.泥頁巖巖石物理建模研究[J].地球物理學(xué)報(bào)，2014，57(6)：1990-1998. Dong Ning,Huo Zhizhou,Sun Zandong,et al.An investigation of a new rock physics model for shale[J].Chinese Journal of Geophysics,2014,57(6):1990-1998.

[28] 李宏兵，張佳佳，姚逢昌.巖石的等效孔隙縱橫比反演及其應(yīng)用[J].地球物理學(xué)報(bào)，2013，56(2)：608-615. Li Hongbing,Zhang Jiajia,Yao Fengchang.Inversion of effective pore aspect ratios for porous rocks and its applications[J].Chinese Journal of Geophysics,2013,56(2):608-615.

[29] Hudson J A.Wave speeds and attenuation of elastic waves in material containing cracks[J].Geophysical Journal Royal Astronomical Society,1981,64(1):133-149.

[30] Berryman J G,Berge P A.Critique of two explicit schemes for estimating elastic properties of multiphase composites[J].Mechanics of Materials,1996,22(2):149-164.

[31] Li Hongbing,Zhang Jiajia.Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory[J].Geophysical Prospecting,2012,60(2):281-292.

[32] Pang J,North C P.Fractals and their applicability in geological wireline log analysis[J].Journal of Petroleum Geology,1996，19(3):339-349.

[33] 郭科，龔灝.多元統(tǒng)計(jì)方法及其應(yīng)用[M].成都：電子科技大學(xué)出版社，2003：65-72. Guo Ke,Gong Hao.Multivariate statistical method and application[M].Chengdu:University of Electronic Science & Technology of China Press,2003:65-72.

[34] 路保平，鮑洪志.巖石力學(xué)參數(shù)求取方法進(jìn)展[J].石油鉆探技術(shù)，2005，33(5)：44-47. Lu Baoping,Bao Hongzhi.Advances in calculation methods for rock mechanics parameters[J].Petroleum Drilling Techniques,2005,33(5):44-47.

[35] 陸黃生.測(cè)井技術(shù)在石油工程中的應(yīng)用分析與發(fā)展思考[J].石油鉆探技術(shù)，2012，40(6)：1-7. Lu Huangsheng.Application and development analysis of well logging information in petroleum engineering[J].Petroleum Drilling Techniques,2012,40(6):1-7.

[編輯 令文學(xué)]

Fracture Identification Method for Marine-Continental Transitional Clastic Rocks Based on the Array Acoustic Logging

Yin Shuai1,2, Ding Wenlong1,2, Zhao Wei1,2, Sun Yuanhui3, Yuan Jiangru3, Cong Sen1,2

(1.SchoolofEnergyResources,ChinaUniversityofGeosciences(Beijing),Beijing,100083,China; 2.KeyLaboratoryforMarineReservoirEvolutionandHydrocarbonAbundanceMechanism,MinistryofEducation,ChinaUniversityofGeosciences(Beijing),Beijing,100083,China; 3.PetroChinaResearchInstituteofPetroleumExploration&Development,Beijing,100083,China)

It is much more difficult to identify fractures in marine-continental transitional facies tight clastic rocks from well logging,so a study was proposed to carry out quantitative identification on single-well fracture development sections in marine-continental transitional facies clastic formations by using array sonic logging data, combined with the Gassmann equation, adaptive matrix mineral and skeleton modulus extraction method, and DEM theoretical model. It is shown that the extractedKo(matrix mineral bulk modulus of clastics)ranges from 13 to 58 GPa, which is slightly higher thanKs(bulk modulus)andKd(dry rock skeleton bulk modulus), and μo(matrix mineral shear modulus) ranges from 5 to 18 GPa, which is slightly higher than μs(shear modulus). Shear wave time difference can be predicted effectively with relative error 3.5% from adaptive method. Fractures can be predicted by introducing modulus parameters extracted by the adaptive method into the DEM theoretical model, and its coincidence rate is apparently superior to that of the conventional logging method, the fracture parameter method,the multi-parameters probability discriminance method and theR/Smethod. In conclusion, the newly proposed fracture identification method is reliable, with its calculation results highly coincident with practical value.

array sonic logging；marine-continental transitional facies；clastic rocks；fracture identification

2015-04-18；改回日期：2015-08-16。

尹帥(1989—)，男，山東新泰人，2010年畢業(yè)于山東科技大學(xué)地質(zhì)工程專業(yè)，在讀博士研究生，主要從事石油構(gòu)造分析與控油作用、非常規(guī)油氣構(gòu)造和裂縫及其與含氣量的關(guān)系等方面的研究。

丁文龍，dingwenlong2006@126.com。

國家自然科學(xué)基金項(xiàng)目“渤海灣盆地濟(jì)陽坳陷古近系陸相富有機(jī)質(zhì)頁巖裂縫研究”(編號(hào)：41372139)、“中國南方下古生界海相富有機(jī)質(zhì)頁巖裂縫發(fā)育程度與主控因素定量關(guān)系研究”(編號(hào)：41072098)及國家科技重大專項(xiàng)課題“煤系地層游離氣富集規(guī)律分布與研究”(編號(hào)：2011ZX05033-004)聯(lián)合資助。

?測(cè)井錄井?

10.11911/syztjs.201505013

P631.8+14

A

1001-0890(2015)05-0075-08

聯(lián)系方式：13375633766，speedysys@163.com。