許君玉, 許新
(1.中國石油化工股份有限公司石油勘探開發(fā)研究院, 北京 100083; 2.北京市海淀區(qū)職業(yè)介紹服務中心, 北京 100195)
紅河油田位于鄂爾多斯盆地天環(huán)坳陷南段,其東部臨伊陜斜坡,西部接西緣沖斷帶,南部為渭北隆起。由于多期地質(zhì)構(gòu)造運動,區(qū)內(nèi)形成了多條斷裂和裂縫。裂縫是控制低滲透油氣藏開發(fā)產(chǎn)能的關(guān)鍵性因素[1],裂縫的存在不僅改善了地層對流體的儲集性能,而且作為滲流通道大大提高了低孔隙度低滲透率儲集層的生產(chǎn)能力[2]。因此,對紅河油田低孔隙度低滲透率儲集層裂縫的研究意義重大。
通過裂縫觀察和對研究區(qū)21口取心井182條裂縫統(tǒng)計,裂縫巖性以砂巖為主,其含量大于90%;裂縫產(chǎn)狀主要為高角度裂縫和垂直裂縫,占80%以上;裂縫開度范圍0.1~0.3 mm;裂縫充填情況為未充填60%,其余為半充填和全充填各占一半,充填物主要是方解石和碳質(zhì)泥屑。裂縫平均密度為0.25條/m,發(fā)育井段的垂直裂縫達到3~6條/m,水平裂縫密度相對較低,約為3條/m。
裂縫在測井曲線上有2種響應特征。①低電阻率特征,即雙感應、八側(cè)向和陣列感應測井曲線均顯示低電阻率值,八側(cè)向與雙感應電阻率值正異常;雙井徑微擴徑,聲波、中子值增大,密度值減小。②高電阻率特征,即雙感應正異常顯示高電阻率值,八側(cè)向電阻率值在層的上部與雙感應電阻率值正差異,在層的下部與雙感應電阻率值負差異;聲波、密度值減小,中子值增大。
裂縫與基質(zhì)的地質(zhì)和地球物理特征差異性,地層中裂縫的存在會引起不同的測井響應,根據(jù)測井響應的不同變化,就可以識別和研究裂縫[3-18]。研究人員利用測井資料對裂縫進行了識別[19-21],但效果欠佳。在對紅河油田裂縫識別的過程中,發(fā)現(xiàn)不同井的裂縫在測井曲線上的響應特征不同,因此利用測井曲線的響應特征識別裂縫比較困難。在對裂縫研究的基礎上,利用常規(guī)測井曲線作為輸入曲線的變尺度分形方法對紅河油田裂縫進行了識別,與成像測井資料對比,效果較好。
對裂縫識別的探索發(fā)現(xiàn),在利用測井資料對裂縫進行識別的方法中,分形幾何學的應用在研究儲層的孔隙結(jié)構(gòu)、儲層的非均質(zhì)性以及裂縫預測等方面取得了豐富的成果[22-23]。R/S分形(變尺度分形)技術(shù)是應用廣泛且成熟的分形統(tǒng)計方法之一[24-26],方法中R為極差,它是最大累積離差與最小累積離差之差,代表了時間序列的復雜程度;S是標準差,它是變差的平方根,代表了時間序列的平均趨勢,R/S之比代表無因次的時間序列的相對波動程度。對于一條對稱的測井曲線z(i)(i=1,2,…,n),R和S的定義為
(1)
(2)
式中,n為逐點分析層段測井采樣點數(shù);u為由端點開始在0~n之間依次增加的采樣點數(shù);i和j為表示采樣點個數(shù)的變量。
圖1 HH1057-3井2 229~2 248 m井段不同測井曲線變尺度分形圖
根據(jù)以上方法建立不同測井曲線的R(n)/S(n)與n的對應關(guān)系,與成像測井資料對比,可找出識別孔洞、裂縫最敏感的測井曲線。
圖2 HH1057-3井裂縫在成像測井和R(n)/S(n)曲線上的響應特征
用變尺度分形方法對紅河油田研究區(qū)進行了裂縫識別。通過對不同測井曲線作變尺度R(n)/S(n)分形圖與成像測井資料對比,發(fā)現(xiàn)不同測井曲線對裂縫識別的敏感程度不同。圖1是對HH1057-3井在2 229~2 248 m井段用不同測井曲線的作的R(n)/S(n)的變尺度分形圖。為避免巖性影響,在選取識別的裂縫井段時參照自然伽馬、自然電位和井徑測井曲線,選取巖性變化不大的井段作為裂縫識別井段。
從圖1所作的各條測井曲線的變尺度分形圖上看出,圖1(a)、圖1(b)和圖1(c)的自然電位、自然伽馬和井徑這3條巖性測井曲線的變尺度分形圖形態(tài)各不相同,只有自然電位曲線的變尺度分形圖上2 241.125~2 243.25 m的井段為平直段,所以圖1(b)和圖1(c)的自然伽馬和井徑曲線對裂縫識別是不敏感的。圖1(d)至圖1(i)為電阻率曲線中感應、深感應、八側(cè)向、淺側(cè)向、深側(cè)向和陣列感應測井曲線的變尺度分形圖,這6條測井曲線的變尺度分形圖形態(tài)相似,只在淺側(cè)向曲線的變尺度分形圖上2 246.25~2 246.75 m井段有一平直段,其余5條電阻率測井曲線變尺度分形圖上均無平直井段,說明電阻率測井曲線中除淺側(cè)向曲線外的其余曲線對裂縫識別不敏感。圖1(j)、圖1(k)和圖1(l)的聲波、中子孔隙度和密度3條孔隙度曲線的變尺度分形圖中,前二者曲線形態(tài)相似,聲波時差曲線的變尺度分形圖上的平直段深度為2 236.125~2 240.25 m;中子孔隙度曲線的變尺度分形圖上的平直段有2段,上段為2 230.25~2 230.625 m,下段為2 234.375~2 240.25 m;密度測井曲線變尺度分形圖上沒有平直段,說明密度測井曲線對裂縫識別不敏感。
圖2左側(cè)為HH1057-3井成像測井圖識別的裂縫井段,識別出4段裂縫井段,第1段為2 230.30~2 230.80 m,第2段為2 233.38~2 236.20 m,第3段為2 237.00~2 239.10 m,第4段為2 239.80~2 242.12 m。圖3右側(cè)為自然電位、淺側(cè)向、聲波和中子孔隙度4條測井曲線所識別的裂縫井段與其左側(cè)的成像測井圖上識別的裂縫井段的對比。圖2右側(cè)自然電位曲線變尺度分形圖上的直線段深度為2 241.125~2 243.25 m,頂部深度比成像測井圖上的裂縫井段第4段底部深度淺0.995 m,與成像測井圖上的裂縫井段重復率為12.86%,所以認為自然電位曲線變尺度分形法對裂縫識別不敏感。淺側(cè)向曲線變尺度分形圖上的直線段深度為2 246.25~2 246.75 m,在成像測井圖上識別的所有裂縫井段的下方,與成像測井圖上識別的裂縫井段沒有重復,所以淺側(cè)向曲線的變尺度分形法對裂縫識別不敏感。聲波測井曲線變尺度分形圖上的直線段深度2 236.125~2 240.25 m,其頂部深度比成像測井圖上裂縫井段的第2段(2 233.38~2 236.20 m)底部深度淺0.075 m,底部深度比成像測井圖上裂縫井段的第3段(2 237.00~2 239.10 m)底部深1.15 m,聲波測井曲線變尺度分形法識別的裂縫井段與成像測井圖上的裂縫井段重復率為43.5%。中子孔隙度測井曲線變尺度分形法識別出2段裂縫,上段(2 230.25~2 230.625 m)頂部深度比成像測井圖上的裂縫井段第1段(2 230.30~2 230.80 m)頂部深度淺0.05 m,底部深度比成像測井圖上裂縫井段第1段的底部深度淺0.175 m;下段(2 234.375~2 240.25 m)頂部深度比成像測井圖上的裂縫井段第2段(2 233.38~2 236.20 m)頂部深度深0.995 m,底部深度比成像測井圖上的裂縫井段第4段(2 239.80~2 242.12 m)底部深度淺1.87 m。中子
孔隙度測井曲線變尺度分形法識別的裂縫井段與成像測井圖上的裂縫井段重復率為85.0%。所以,變尺度分形法識別裂縫中,中子孔隙度測井曲線是識別裂縫最敏感的曲線。
用變尺度分形方法與常規(guī)測井曲線對裂縫的響應特征相結(jié)合對紅河油田研究區(qū)內(nèi)31口井進行了裂縫識別,裂縫識別率為80%以上。表1為有裂縫的巖心和成像測井裂縫井段與變尺度分形方法識別的裂縫井段頂、底深度對比。
表1 紅河油田取心井和成像測井裂縫井段與變尺度分形法識別裂縫井段對比表
表1(續(xù))
井名取心或成像測井井段/m變尺度分形法裂縫識別井段/m識別裂縫距離取心或成像裂縫頂部深度/m識別裂縫距離取心或成像裂縫底部深度/m裂縫獲取方式HH262128.20~2128.802128.625~2129.250.4250.45成像HH1057-32191.80~2193.202191.625~2192.375-0.175-0.825成像HH1057-32193.20~2194.802193.875~2194.8750.6750.075成像HH1057-32279.60~2282.002278.875~2280.875-0.725-1.125成像HH1057-32297.70~2298.602298.75~2299.1251.050.525成像
用變尺度分形方法對紅河油田研究區(qū)內(nèi)試油資料中獲得工業(yè)油流的井日產(chǎn)液量對比驗證,按已有試油資料中日產(chǎn)液量由大到小排序,將識別的裂縫井列表于表2。由表2可知,日產(chǎn)液量大的井多數(shù)存在裂縫,說明了裂縫的存在對儲層產(chǎn)量起著關(guān)鍵作用;另一方面,在識別的14口高產(chǎn)井中有11口井識別出了裂縫,說明了變尺度分形方法對裂縫識別的有效性。
表2 紅河油田儲層產(chǎn)液量與變尺度分形法識別裂縫
(1) 紅河油田發(fā)育高角度裂縫和垂直裂縫。不同測井曲線對裂縫的測井響應特征不同,用常規(guī)測井曲線很難識別裂縫的存在。
(2) 裂縫在變尺度分形曲線上的響應特征表現(xiàn)為曲線偏離原來的變化趨勢,成為近似一直線段。變尺度分形法裂縫識別中,不同測井曲線對裂縫響應的敏感程度不同,中子孔隙度測井曲線是識別裂縫最敏感的測井曲線,這種方法裂縫識別率較高。
(3) 變尺度分形方法識別裂縫操作簡單、成本費用低,是識別裂縫較經(jīng)濟的方法。
(4) 利用變尺度分形方法識別的裂縫與巖心或成像測井識別的裂縫頂、底深度有一定的誤差。因此,綜合利用各種方法識別裂縫,才能最大限度地提高裂縫識別率。
參考文獻:
[1] 徐朝暉, 徐懷民, 林軍, 等. 常規(guī)測井資料識別砂礫巖儲集層裂縫技術(shù) [J]. 科技導報, 2008, 26(7): 34-37.
[2] 周燦燦, 楊春頂. 砂巖裂縫的成因及其常規(guī)測井資料綜合識別技術(shù)研究 [J]. 石油地球物理勘探, 2003, 38(4): 425-430.
[3] 雍世和. 測井資料綜合解釋與數(shù)字處理 [M]. 北京: 石油工業(yè)出版社, 1999: 146-168.
[4] 葉濤, 蔣有錄, 劉華, 等. 渤南洼陷裂縫成因類型及測井識別方法 [J]. 西安石油大學學報: 自然科學版, 2013, 28(1): 52-60.
[5] 王慶如, 李敬功. 碳酸鹽巖氣藏儲量參數(shù)測井評價方法 [J]. 巖性油氣藏, 2013, 25(6): 98-111.
[6] 王永剛. 濟陽坳陷太古界變質(zhì)巖儲層裂縫識別與定量解釋 [J]. 測井技術(shù), 2012, 36(6): 590-595.
[7] 李健偉, 王磊, 陳利雯, 等. 碳酸鹽巖儲層裂縫識別技術(shù)研究 [J]. 工程地球物理學報, 2013, 10(6): 885-890.
[8] 姚勝林, 馬玉甲. 碳酸鹽巖儲層裂縫識別技術(shù)研究 [J]. 工程地球物理學報, 2013, 10(6): 885-890.
[9] 蒲靜, 秦啟榮. 油氣儲層裂縫預測方法綜述 [J]. 特種油氣藏, 2008, 15(3): 9-13.
[10] 王珂, 張榮虎, 戴俊生, 等. 低滲透儲層裂縫研究進展 [J]. 地球科學與環(huán)境學報, 2015, 37(2): 44-58.
[11] 吳瓊, 林冬萍, 于春燕, 等. 新立油田低滲透油層裂縫測井識別方法 [J]. 大慶石油地質(zhì)與開發(fā), 2007, 26(2): 112-115.
[12] 鄧少貴, 仝兆岐, 范宜仁, 等. 致密砂巖儲集層裂縫的雙側(cè)向測井響應快速計算方法 [J]. 中國石油大學學報: 自然科學版, 2005, 29(3): 31-34.
[13] 柯式鎮(zhèn), 孫貴霞. 井壁電成像測井資料定量評價裂縫的研究 [J]. 測井技術(shù), 2002, 26(2): 101-103.
[14] 肖承文, 朱筱敏, 李進福, 等. 高壓低滲致密裂縫性砂巖測井評價技術(shù) [J]. 新疆石油地質(zhì), 2007, 28(6): 761-763.
[15] 李軍, 張超謨, 肖承文, 等. 庫車地區(qū)砂巖裂縫測井定量評價方法及應用 [J]. 天然氣工業(yè), 2008, 28(10): 25-27.
[16] 鄭軍, 劉鴻博, 周文, 等. 阿曼五區(qū)塊Daleel油田儲層裂縫識別方法研究 [J]. 測井技術(shù), 2010, 34(3): 251-256.
[17] 王瑞雪, 談順佳, 張曉峰. 常規(guī)測井方法在川西南地區(qū)儲層裂縫識別中的應用 [J]. 石油地質(zhì)與工程, 2014, 28(1): 41-43.
[18] 朱雷, 章成廣, 居大海, 等. 利用聲波全波列測井資料識別致密砂巖儲層裂縫 [J]. 長江大學學報: 自科版, 2015, 12(8): 23-26.
[19] 王瑞雪, 張曉峰, 談順佳, 等. 基于成像測井資料多種濾波方法在裂縫識別中的應用 [J]. 測井技術(shù), 2015, 39(2): 155-159.
[20] 趙杰, 李霞, 楊雪冰. 古龍地區(qū)青山口組泥巖裂縫有效儲層識別方法 [J]. 大慶石油地質(zhì)與開發(fā), 2013, 32(4): 147-150.
[21] 楊雪, 王瑩, 李萬才. 應用常規(guī)測井資料進行裂縫識別 [J]. 吉林地質(zhì), 2013, 32(1): 102-105.
[22] 紀發(fā)華, 張一偉. 分形幾何學在儲層非均質(zhì)性描述中的應用 [J]. 中國石油大學學報: 自然科學版, 1994, 18(5): 161-168.
[23] 王廣才, 章成廣, 居大海, 等. 分數(shù)維在識別深層致密砂巖裂縫中的應用 [J]. 長江大學學報: 自科版, 2015, 12(8): 33-36.
[24] 師光輝, 李永權(quán), 王魯. 結(jié)合變尺度分析方法和常規(guī)測井資料識別裂縫的發(fā)育 [J]. 石油儀器, 2013, 27(1): 49-51.
[25] 劉麗麗, 趙中平, 李亮, 等. 變尺度分形技術(shù)在裂縫預測和儲層評價中的應用 [J]. 石油與天然氣地質(zhì), 2008, 29(1): 31-37.
[26] 毛寧波.R/S分析估算儲層參數(shù)分維研究新進展 [J]. 石油物探, 2000, 39(3): 108-126.