瓊東南盆地北部天然氣水合物重點(diǎn)礦體預(yù)測(cè)及定量評(píng)價(jià)

路允乾,孟大江,文鵬飛,張寶金,2

(1. 廣州海洋地質(zhì)調(diào)查局,廣州 510075;2. 南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(廣州),廣州 511458;3. 天然氣水合物勘查開(kāi)發(fā)國(guó)家工程研究中心,廣州 510075)

0 引言

天然氣水合物是一種資源潛力巨大的新型清潔能源,是由天然氣和水在高壓低溫環(huán)境條件下形成的冰態(tài)、結(jié)晶狀籠形化合物[1-5]。近20年來(lái),中國(guó)在南海北部海域開(kāi)展了大量水合物調(diào)查和研究,在珠江口盆地神狐海域以及瓊東南盆地均發(fā)現(xiàn)了水合物賦存的證據(jù)[6-8]。瓊東南盆地是在準(zhǔn)被動(dòng)大陸邊緣構(gòu)造地質(zhì)背景下形成的富油氣盆地,近年來(lái)在瓊東南盆地深水區(qū)的勘探活動(dòng)中,發(fā)現(xiàn)了一系列與水合物賦存有關(guān)的地質(zhì)、地球物理及地球化學(xué)異常,廣州海洋地質(zhì)調(diào)查局于2018年在瓊東南盆地東部首次實(shí)施水合物深部鉆探并獲得水合物實(shí)物樣品[9-12]。前人研究及勘探實(shí)踐證明,瓊東南盆地油氣資源豐富,具備天然氣水合物形成的成藏地質(zhì)條件,水合物資源潛力巨大,且水合物層下方往往伴生大量淺層氣,具有良好的資源前景和勘探潛力[10-18];但因研究區(qū)仍處在較低的勘探開(kāi)發(fā)階段,目前對(duì)天然氣水合物的地球物理識(shí)別主要依靠地震剖面上的似海底反射波(BSR),由于BSR存在的多解性問(wèn)題,尚存在水合物分布不清的問(wèn)題,對(duì)水合物及其氣態(tài)烴儲(chǔ)層物性的研究也較為缺乏。前人在利用體積法預(yù)測(cè)水合物資源量時(shí),往往通過(guò)統(tǒng)計(jì)水合物層的平均厚度和大致分布面積來(lái)獲取水合物層體積,并在這個(gè)體積內(nèi)計(jì)算水合物平均的物性參數(shù),這種方法存在較大誤差。水合物的賦存不需要特定的圈閉,其分布和形態(tài)受溫壓穩(wěn)定帶控制,水合物本身就可以作為下賦氣態(tài)烴的蓋層,因此,精確刻畫(huà)水合物及氣態(tài)烴礦體的展布形態(tài)和厚度并計(jì)算其物性參數(shù),對(duì)于提升水合物資源量預(yù)測(cè)的準(zhǔn)確度十分重要。

該文在水合物成藏地質(zhì)條件分析的基礎(chǔ)上,基于水合物勘探評(píng)價(jià)的實(shí)際,以瓊東南盆地北部重點(diǎn)礦區(qū)三維地震資料為依托,尋找指示水合物儲(chǔ)層的地球物理響應(yīng)特征,利用疊后波阻抗反演對(duì)研究區(qū)水合物及氣態(tài)烴目標(biāo)儲(chǔ)層進(jìn)行精細(xì)刻畫(huà),預(yù)測(cè)水合物及氣態(tài)烴的分布范圍并計(jì)算其厚度和物性參數(shù),預(yù)測(cè)有利成礦區(qū)。該研究有利于認(rèn)識(shí)該區(qū)域水合物的分布特征及富集成藏位置,為瓊東南盆地北部水合物勘探開(kāi)發(fā)提供依據(jù)。

1 區(qū)域地質(zhì)背景

瓊東南盆地位于東經(jīng)108°50′～111°50′、北緯15°30′～19°00′的海域,在海南島南部陸緣與西沙群島之間,海域面積約6.3×104km2。盆地發(fā)育于南海北部陸緣西側(cè)的華南陸塊、中—西沙陸塊和印支陸塊的結(jié)合部位,屬南海北部準(zhǔn)被動(dòng)大陸邊緣伸展裂谷型盆地[9,11,19]。瓊東南盆地西寬東窄、北寬南窄,在地形上跨越陸架和陸坡區(qū),海水由北西向東南變深,陸架區(qū)水深變化小,一般為90～200 m,往陸坡至海槽區(qū)水體急劇加深,從200 m迅速加深到約2 500 m。圖1所示為瓊東南盆地構(gòu)造單元?jiǎng)澐?、綜合地層柱狀圖及礦區(qū)位置。瓊東南盆地呈北東走向展布,由一系列北東向和北北東向的凹陷和隆起組成[20-21](如圖1a所示)。從區(qū)域構(gòu)造演化歷史看,瓊東南盆地大致經(jīng)歷了裂陷期和裂后坳陷期2個(gè)演化時(shí)期。裂陷期始于始新世,由于印度洋板塊的北向擠壓,新特提斯洋殼俯沖于華南—印支陸塊的陸殼下,南海發(fā)生珠瓊運(yùn)動(dòng)一幕,基底收到拉張作用開(kāi)始張裂,南海北部發(fā)生抬升與剝蝕,并伴隨斷裂和巖漿作用;珠瓊運(yùn)動(dòng)二幕使南海北部再次遭受強(qiáng)烈抬升和剝蝕。早中新世至中中新世末為坳陷早期,南海運(yùn)動(dòng)使盆地北部發(fā)生區(qū)域性抬升和強(qiáng)烈剝蝕,形成中新世與漸新世之間的巨大不整合面。坳陷晚期開(kāi)始于晚中新世,持續(xù)至第四紀(jì),盆地以坳陷作用為主,地殼相對(duì)穩(wěn)定[22-25]。

圖1 瓊東南盆地構(gòu)造單元?jiǎng)澐?、綜合地層柱狀圖及礦區(qū)位置Fig.1 Structural unit division, comprehensive stratigraphic histogram and location of mining area of Qiongdongnan Basin

瓊東南盆地古近系沉積以21 Ma巨大不整合為界,可以劃分出古近系的裂陷沉積和新近系的坳陷沉積2個(gè)沉積構(gòu)造層[21,27](如圖1b所示)。瓊東南盆地古近系斷層較為發(fā)育,新近系斷層少,規(guī)模小,且多表現(xiàn)為層間斷層。盆地中新統(tǒng)以來(lái)斷層可分為北東、北西和北北東向3組斷裂,其中北東向斷層規(guī)模大且斷距大,在區(qū)域上主要分布于中央坳陷帶內(nèi)的凸起與凹陷的分界位置,主要活動(dòng)時(shí)間為早—中中新世,終止于晚中新世時(shí)期,推測(cè)其受控于古近系南海北部廣泛發(fā)育的南北向拉張大環(huán)境,部分為基底繼承性斷裂。相對(duì)而言,北西向斷層規(guī)模小,垂向斷距小,斷層開(kāi)始活動(dòng)時(shí)間晚,集中于早—中中新世,終止于晚中新世,推測(cè)其與盆地西部的紅河—鶯歌海大斷裂新近紀(jì)左旋剪切活動(dòng)有關(guān)。

瓊東南盆地主要烴源巖為陵水組、三亞組和鶯歌海組,油氣勘探表明,臨近研究區(qū)的主要生烴凹陷包括寶島凹陷和松南凹陷,其烴源巖成熟度Ro多數(shù)大于2.0%,以生干氣為主,因此研究區(qū)具有非常好的氣源條件[21]。研究區(qū)上新統(tǒng)—第四系的粗粒沉積物和塊體流沉積為水合物的富集和成藏提供了有利的儲(chǔ)層條件[12,20]。瓊東南盆地具有良好的氣源供給條件,研究區(qū)位于寶島凹陷內(nèi),來(lái)自早期盆地高壓的作用,以及寶島凹陷內(nèi)火山活動(dòng)頻繁,這些因素導(dǎo)致研究區(qū)內(nèi)底部流體在超壓作用下突破地層形成底辟構(gòu)造。底辟作用攜帶深部流體運(yùn)移到淺部地層,底辟構(gòu)造、氣煙囪和斷裂組成了烴類(lèi)向上運(yùn)移的網(wǎng)絡(luò)輸導(dǎo)體系,深部熱解氣源巖與淺部微生物氣源巖構(gòu)成了良好的立體供烴體系,可為天然氣水合物的形成提供充足的氣源供給[29-30]。瓊東南盆地深水區(qū)的天然氣水合物成藏條件優(yōu)良,勘探實(shí)踐表明,瓊東南盆地資源量巨大,是我國(guó)重要的海上油氣富集區(qū)以及天然氣水合物成藏遠(yuǎn)景區(qū)。

2 水合物成藏系統(tǒng)的地球物理特征

2.1 地震反射特征

BSR(Bottom Simulating Reflector)即似海底反射,通常代表水合物分布的底界,是識(shí)別水合物存在最直接的標(biāo)志。在地震剖面上,BSR一般位于海底以下1 000 ms的范圍內(nèi)(雙程旅行時(shí))。BSR是一個(gè)近似于海底的反射面,其上出現(xiàn)的振幅空白帶也是指示天然氣水合物存在的重要標(biāo)志[13,31-32]。研究區(qū)的三維地震資料上可以識(shí)別出明顯的BSR反射特征,以礦體A和礦體B最為典型(位置見(jiàn)圖1c),表現(xiàn)為1套穩(wěn)定的中-強(qiáng)振幅和中-高頻率反射,橫向連續(xù)性好,大致與海底平行,極性與海底相反(如圖2所示)。從BSR的垂向分布深度來(lái)看,BSR處于海底以下200～300 ms(雙程旅行時(shí)),通過(guò)時(shí)深轉(zhuǎn)換計(jì)算,研究區(qū)內(nèi)BSR位于海底以下120～359 m。BSR界面之下的空白反射指示了氣態(tài)烴的存在,氣態(tài)烴通常賦存于天然氣水合物之下,水合物作為其封蓋層存在。

圖2 瓊東南盆地北部似海底反射層(BSR)地震反射特征Fig.2 Seismic reflection characteristics of bottom simulating reflector(BSR) in northern Qiongdongnan Basin

2.2 礦體雕刻

BSR的存在與水合物分布并不是一一對(duì)應(yīng)的關(guān)系,深水海底淺層含砂層也常表現(xiàn)為與海底極性相反的強(qiáng)地震反射特征;若沉積物孔隙中充填少量氣體,也會(huì)產(chǎn)生強(qiáng)烈的地層反射;早先生成的水合物因?yàn)闃?gòu)造作用或溫壓條件發(fā)生改變導(dǎo)致泄漏逸散,水合物藏造成破壞,也會(huì)存有BSR標(biāo)志[8,33];所以?xún)H用BSR來(lái)評(píng)價(jià)水合物的發(fā)育分布,總會(huì)有很多不確定性。通常儲(chǔ)層含天然氣水合物后具有較高的速度特征,含氣態(tài)烴后具有較低的速度特征,二者可以形成較大的波阻抗差異,因此可以通過(guò)波阻抗刻畫(huà)出水合物和氣態(tài)烴礦體。根據(jù)南海北部標(biāo)準(zhǔn)化后的測(cè)井波阻抗統(tǒng)計(jì)直方圖(如圖3所示),水合物和氣態(tài)烴地層波阻抗值存在明顯差異,水合物地層縱波阻抗值大多為3 100～4 500 (m/s)×(g/cm3),氣態(tài)烴地層縱波阻抗值均小于2 900 (m/s)×(g/cm3),表明可以利用波阻抗反演來(lái)識(shí)別天然氣水合物,從而預(yù)測(cè)天然氣水合物的分布范圍。

圖3 波阻抗頻率分布直方圖Fig.3 Distribution histogram of wave impedance frequency

該文采用約束的稀疏脈沖反演方法,在地質(zhì)解釋層位的約束下利用疊前深度偏移速度體構(gòu)建低頻模型,在三維地震體上提取零相位子波,在建立合理的低頻模型并提取高質(zhì)量子波的基礎(chǔ)上,對(duì)反演參數(shù)進(jìn)行測(cè)試調(diào)整,如地震記錄與子波的比例因子、模型可信度以及地震數(shù)據(jù)可信度等,最終優(yōu)選出最符合實(shí)際地質(zhì)情況的反演參數(shù)。通過(guò)對(duì)反演參數(shù)的反復(fù)測(cè)試,發(fā)現(xiàn)反演子波比例因子選擇4×10-6并且模型可信度為30%時(shí)反演結(jié)果最為合理。

波阻抗剖面能更清晰地凸顯含水合物地層的高波阻抗特征和下部氣態(tài)烴的低阻抗特征,水合物和氣態(tài)烴的頂、底界面清楚,可以直觀地看出礦體的空間展布形態(tài)。圖4所示為瓊東南盆地北部礦體波阻抗反演剖面,紅色高阻抗指示了高飽和度水合物礦體分布,黃綠色低阻抗指示了相對(duì)低飽和度水合物礦體分布,白色和黑色虛線(xiàn)分別代表了水合物礦體的頂?shù)捉缑妗姆囱萁Y(jié)果看,以BSR(黑色虛線(xiàn))為界,水合物高阻抗特征和氣態(tài)烴低阻抗特征明顯,天然氣水合物儲(chǔ)層波阻抗值超過(guò)3 100 (m/s)×(g/cm3),氣態(tài)烴波阻抗值均小于2 900 (m/s)×(g/cm3)。礦體A水合物波阻抗值為3 100～4 000 (m/s)×(g/cm3),水合物大多聚集在1.6 s以上,呈條帶狀分布(如圖4a所示);礦體B水合物層與地層平行分布,波阻抗值為3 800～5 500 (m/s)×(g/cm3),水合物聚集在1.9～2.1 s(如圖4b所示)。氣態(tài)烴表現(xiàn)為藍(lán)色低縱波阻抗特征,波阻抗值均小于2 900 (m/s)×(g/cm3),賦存于天然氣水合物礦體底部,礦體A和礦體B底部均有較厚的氣態(tài)烴層發(fā)育。

圖4 瓊東南盆地北部礦體波阻抗反演剖面Fig.4 Wave impedance inversion profile of ore body in northern Qiongdongnan Basin

3 水合物及氣態(tài)烴儲(chǔ)層物性預(yù)測(cè)

3.1 儲(chǔ)層物性預(yù)測(cè)分析

海洋天然氣水合物儲(chǔ)層特性直接影響天然氣水合物富集程度,儲(chǔ)層飽和度是影響天然氣水合物資源量的關(guān)鍵參數(shù)。圖5所示為水合物及氣態(tài)烴飽和度與波阻抗巖石物理量板。儲(chǔ)層含天然氣水合物后波阻抗值變大,且波阻抗值隨著儲(chǔ)層水合物飽和度的增大而增大,二者具有良好的線(xiàn)性關(guān)系[34]。根據(jù)南海北部天然氣水合物鉆井實(shí)測(cè)資料,繪制了天然氣水合物飽和度與波阻抗交匯圖(如圖5a所示),確立了水合物飽和度與波阻抗的量化關(guān)系為:

圖5 水合物及氣態(tài)烴飽和度與波阻抗巖石物理量版Fig.5 Petrophysical template of hydrate and gaseous hydrocarbon saturation and wave impedance

Sathy=0.000 339×AI-0.976

其中:Sathy為水合物飽和度,AI為波阻抗,二者相關(guān)系數(shù)約為0.74。

氣態(tài)烴飽和度與波阻抗的關(guān)系受地層孔隙度影響,因此需要建立不同孔隙度下的氣態(tài)烴波阻抗與飽和度的擬合關(guān)系。研究區(qū)礦物成分參考了鉆井的取心分析資料取平均獲得,其中石英含量為43.2%,斜長(zhǎng)石含量為6.8%,方解石含量為18.0%,黏土含量為32.0%,由此建立了氣態(tài)烴飽和度和波阻抗的巖石物理量版(如圖5b所示)。研究區(qū)氣態(tài)烴平均孔隙度為40%～50%,該文采用40%的孔隙度曲線(xiàn)預(yù)測(cè)氣態(tài)烴礦體飽和度。

3.2 水合物及氣態(tài)烴的厚度和飽和度預(yù)測(cè)

研究區(qū)三維地震資料及波阻抗反演結(jié)果顯示,礦體A和礦體B是水合物及氣態(tài)烴分布的有利區(qū)。根據(jù)礦體雕刻結(jié)果,人工拾取這2個(gè)礦體范圍內(nèi)的海底、BSR界面以及氣態(tài)烴底界并提取礦體的平面分布特征。篩選出海底與BSR界面之間大于3 100 (m/s)×(g/cm3)的波阻抗值,并統(tǒng)計(jì)它們的時(shí)間域厚度;利用相同的方法,在BSR界面與氣態(tài)烴底界之間篩選出小于2 900 (m/s)×(g/cm3)的波阻抗值,得到氣態(tài)烴的時(shí)間域厚度。圖6所示為礦體A和礦體B物性參數(shù)平面切片。

圖6 礦體A和礦體B物性參數(shù)平面切片(暖色代表高數(shù)值,冷色代表低數(shù)值)Fig.6 Plane section of physical property parameters of ore body A and ore body B(Warm colors represent high values and cool colors represent low values)

時(shí)間域厚度與疊前深度偏移速度相乘可得到礦體厚度(如圖6a、圖6c、圖6e和圖6g所示)。礦體A和B水合物與氣態(tài)烴飽和度則根據(jù)圖5巖石物理量版求取(如圖6b、圖6d、圖6f和圖6h所示)。

礦體A和B均位于海底隆起的構(gòu)造脊上,證明海底低凸起構(gòu)造是水合物和氣態(tài)烴賦存的有利區(qū)。礦體A水合物平均厚度為15～50 m,厚度中心位于礦體西部,最厚處超過(guò)50 m(如圖6a所示);水合物平均飽和度為0.1～0.4,高飽和度水合物主要集中在礦體西南部(如圖6b所示)。礦體A氣態(tài)烴平均厚度為30～100 m,厚度中心位于礦體中南部,厚度超過(guò)100 m(如圖6c所示);氣態(tài)烴平均飽和度約為0.1～0.3,高飽和度氣態(tài)烴同樣位于礦體中南部(如圖6d所示)。礦體B水合物平均厚度為15～50 m,由北向南厚度逐漸增大,厚度中心位于礦體南部(如圖6e所示);水合物平均飽和度為0.2～0.4,飽和度同樣具有南高北低的特征(如圖6f所示)。礦體B氣態(tài)烴平均厚度為10～80 m,北部氣態(tài)烴厚度大于南部;氣態(tài)烴平均飽和度較低,分布較為分散。綜上所述,礦體A預(yù)測(cè)的高飽和度水合物及氣態(tài)烴主要集中于礦體的中南部,有利探測(cè)區(qū)為A-1井和A-2井;礦體B預(yù)測(cè)的水合物飽和度較高,但氣態(tài)烴飽和度較低,有利探測(cè)區(qū)為B-1井和B-2井。

廣州海洋地質(zhì)調(diào)查局在瓊東南盆地北部A-1井附近進(jìn)行了鉆探,礦體A的波阻抗反演與測(cè)井波阻抗疊合如圖7所示,圖中黑色曲線(xiàn)為測(cè)井波阻抗曲線(xiàn)?？梢钥闯龇囱萁Y(jié)果與測(cè)井資料吻合良好,天然氣水合物表現(xiàn)為高波阻抗特征,反演剖面清晰地刻畫(huà)出了水合物的展布形態(tài),水合物賦存在海底以下約200 ms的溫壓穩(wěn)定帶,大致和海底平行分布。同時(shí),反演剖面也刻畫(huà)出了氣態(tài)烴的響應(yīng)特征,氣態(tài)烴賦存在水合物層之下,水合物可作為其封蓋層。A-1井鉆遇的水合物層下方存在明顯的低阻抗響應(yīng),而從測(cè)井解釋上可以看出低阻抗為氣態(tài)烴特征。鉆井結(jié)果顯示A-1井鉆遇水合物厚度為16.5 m,平均飽和度為14.3%;鉆遇氣態(tài)烴厚度為51 m,平均飽和度為8.8%,與該文的預(yù)測(cè)結(jié)果吻合。

圖7 瓊東南盆地北部礦體A波阻抗反演與測(cè)井波阻抗疊合Fig.7 Wave impedance inversion and logging wave impedance superposition of ore body A in northern Qiongdongnan Basin

4 結(jié)論

1)瓊東南盆地北部成藏條件優(yōu)越,研究區(qū)重點(diǎn)礦體的水合物BSR反射特征明顯,水合物層之下有氣態(tài)烴發(fā)育;水合物地層縱波阻抗通常情況下大于3 100 (m/s)×(g/cm3),氣態(tài)烴地層縱波阻抗通常小于2 900 (m/s)×(g/cm3),根據(jù)二者的波阻抗差異可以精確刻畫(huà)出礦體的分布范圍和展布形態(tài)。

2)通過(guò)約束稀疏脈沖反演得到的反演結(jié)果可靠,波阻抗異常明顯,地質(zhì)現(xiàn)象清晰,與鉆井解釋結(jié)果吻合良好,能有效地反映研究區(qū)天然氣水合物礦體的空間分布范圍,為后續(xù)水合物儲(chǔ)層特征分析及其他綜合地球物理分析提供依據(jù)。

3)研究區(qū)三維地震資料及波阻抗反演結(jié)果顯示,礦體A和B是水合物及氣態(tài)烴分布的有利礦區(qū)。礦體A預(yù)測(cè)的高飽和度水合物及氣態(tài)烴主要集中于礦體的中南部;礦體B預(yù)測(cè)的水合物飽和度較礦體A高,但氣態(tài)烴飽和度較低。海底隆起帶是水合物及氣態(tài)烴勘探的有利地區(qū)。