哈拉湖地區(qū)低頻可控震源天然氣水合物地球物理響應(yīng)特征研究

田玉昆，李桂林，劉暉*，祝有海，馬彥彥，李娟，周惠，王德潤，王林

1 中國地質(zhì)調(diào)查局非常規(guī)油氣地質(zhì)重點(diǎn)實(shí)驗(yàn)室，北京 1000292 中國石油集團(tuán)東方地球物理公司吐哈物探處，新疆哈密 8390093 中國石油集團(tuán)東方地球物理公司國際勘探事業(yè)部，河北涿州 072750

?

哈拉湖地區(qū)低頻可控震源天然氣水合物地球物理響應(yīng)特征研究

田玉昆1，李桂林1，劉暉1*，祝有海1，馬彥彥1，李娟1，周惠1，王德潤2，王林3

1 中國地質(zhì)調(diào)查局非常規(guī)油氣地質(zhì)重點(diǎn)實(shí)驗(yàn)室，北京 1000292 中國石油集團(tuán)東方地球物理公司吐哈物探處，新疆哈密 8390093 中國石油集團(tuán)東方地球物理公司國際勘探事業(yè)部，河北涿州 072750

我國陸域天然氣水合物主要分布于青藏高原和漠河地區(qū).由于永久凍土層的存在，地震勘探很難獲得高品質(zhì)的資料，給天然氣水合物勘探帶來了諸多困難.為解決凍土層對(duì)地震信號(hào)的衰減問題，在哈拉湖地區(qū)采用低頻可控震源進(jìn)行地震資料采集試驗(yàn)，通過提高覆蓋次數(shù)，獲得了較高信噪比的地震資料.在高質(zhì)量地震資料基礎(chǔ)上，進(jìn)行精細(xì)速度分析，獲得了較準(zhǔn)確的疊加速度譜資料；然后以層速度剖面為基礎(chǔ)建立正演模型，開展天然氣水合物地震正演模擬研究；最后利用疊后偏移地震數(shù)據(jù)進(jìn)行地震屬性分析.通過正演模擬和地震屬性綜合研究，總結(jié)了天然氣水合物的地球物理響應(yīng)特征，速度突變和空白反射帶可作為哈拉湖地區(qū)陸域天然氣水合物識(shí)別的敏感因素.

天然氣水合物；低頻可控震源；凍土區(qū)；縱波速度；地震屬性;哈拉湖地區(qū)

1 引言

1965年，蘇聯(lián)在西伯利亞Messoyakha油氣田區(qū)首次發(fā)現(xiàn)天然產(chǎn)出的水合物(張洪濤等，2007).自此以后，各國科學(xué)家對(duì)凍土區(qū)水合物的類型、物理化學(xué)性質(zhì)、地質(zhì)成因、勘探手段、識(shí)別方法、資源評(píng)價(jià)等進(jìn)行了一系列的研究.目前，國外凍土區(qū)水合物研究區(qū)塊主要集中在環(huán)北極地區(qū)的美國阿拉斯加北部斜坡、加拿大Mackenzie三角洲以及俄羅斯Messoyakha油氣田等(郭星旺,2011).

1972年，美國在阿拉斯加北坡普拉德霍見西北部Eileen State-2井中首次從永久凍土層中取到包含水合物的巖心(Collett,1993).同年，加拿大西北地區(qū)的Mackenzie三角洲凍土區(qū)的Mallik L-38井的測井資料也顯示有水合物的存在.此后，各國科學(xué)家對(duì)水合物進(jìn)行了多學(xué)科多領(lǐng)域的研究(Moridis et al.,2004;Dallimore et al.,2005;Pratt et al.,2005;Craven et al.,2012).

國內(nèi)的陸域天然氣水合物研究起步較晚.2002年以來，中國地質(zhì)調(diào)查局先后設(shè)立4個(gè)調(diào)查研究項(xiàng)目，包括“青藏高原多年凍土區(qū)天然氣水合物地球化學(xué)勘查預(yù)研究”、“青藏鐵路沿線天然氣水合物遙感識(shí)別標(biāo)志研究”、“我國陸域永久凍土帶天然氣水合物資源遠(yuǎn)景調(diào)查”和“陸地永久凍土天然氣水合物鉆探技術(shù)研究”(郭星旺,2011)，初步圈定了青藏高原羌塘盆地、祁連山木里地區(qū)和東北漠河盆地等為水合物成礦條件和找礦前景較好的地區(qū).2008年，在祁連山凍土區(qū)成功獲得水合物樣品.

天然氣水合物與常規(guī)油氣一樣，均存儲(chǔ)在具有連通孔隙、裂縫等的巖石或沉積物中，因此水合物儲(chǔ)層必然會(huì)與周圍圍巖存在一定的差異.地震勘探是一種有效的探明地下異常巖層的方法，能夠獲得地下地質(zhì)體的各種信息(Yilmaz,2001;Yuan et al.,2015).但由于極地和高原地區(qū)特殊的自然環(huán)境，很難在陸域水合物遠(yuǎn)景區(qū)開展地震勘探工作.目前，在發(fā)現(xiàn)陸域水合物的地區(qū)中，僅在加拿大的Mackenzie三角洲Mallik地區(qū)和我國的青海木里地區(qū)開展過地震勘探工作(徐明才等,2012).但是由于凍土層的影響，常規(guī)炸藥震源激發(fā)采集的地震數(shù)據(jù)信噪比較低，地震波能量被屏蔽和吸收衰減(Yuan et al.,2016)，嚴(yán)重影響地震勘探的效果.

為了解決凍土區(qū)二維地震勘探資料不理想的問題，查明青海哈拉湖水合物研究區(qū)地層、構(gòu)造特征，探測與天然氣水合物有關(guān)的地震異常響應(yīng)，在青海哈拉湖地區(qū)部署了低頻可控震源二維地震采集工程，地震資料采集觀測系統(tǒng)為3000-10-10-10-3000，低頻可控震源掃描頻率為1.5～96 Hz.研究區(qū)位于哈拉湖坳陷的北部，在大地構(gòu)造上處于南祁連地槽褶皺帶，海拔在4150～4600 m之間，地形整體上較為平坦，適合部署可控震源.哈拉湖坳陷是祁連山西部的一個(gè)山間盆地，位于南祁連盆地中心，屬高山內(nèi)陸湖泊.湖面海拔高度4077 m，最大水深65 m，平均27.4 m.湖泊呈北西西—南東東延伸，長32 km，平均寬度13 km，面積為580 km2.哈拉湖盆地北部有疏勒南山，海拔在5000 m以上，南部為哈拉湖南山，平均海拔不到5000 m，盆地東西兩側(cè)的分水嶺均為低平的丘陵地帶.南北二山由上泥盆統(tǒng)、石炭系和三疊系地層組成，盆地內(nèi)部除零星出露第三紀(jì)地層外，大部分被第四系洪積、冰水沖積物覆蓋.

通過地震采集工作展開的實(shí)驗(yàn)攻關(guān)，獲得了較高品質(zhì)的地震資料，并對(duì)其進(jìn)行了精細(xì)的處理解釋，重點(diǎn)研究了水合物可能分布區(qū)域內(nèi)地震資料的速度、振幅與頻率特點(diǎn).對(duì)整個(gè)研究區(qū)的地震數(shù)據(jù)進(jìn)行了統(tǒng)計(jì)，并提取了頻率和振幅屬性.歸納參數(shù)異常，為預(yù)測天然氣水合物提供參考依據(jù).

2 天然氣水合物地球物理特征分析

明確天然氣水合物的物理性質(zhì)是估算自然界中天然氣水合物含量的必要條件.從組成上來說，天然氣容納于由水分子構(gòu)成的固體晶格中.同為固體，這就意味著天然氣水合物與冰在物理性質(zhì)上必然會(huì)存在著一定的相似性.兩者僅在縱波速度上存在一定差異，冰的縱波速度為3.8 km·s-1，結(jié)構(gòu)I型天然氣水合物為3.3 km·s-1(李新和肖立志,2013).

自然界中，天然氣水合物所在地層多為未固化或者半固化的沉積層.在測井上具有典型的響應(yīng)特征，例如固態(tài)水合物的電阻率較高(蔣國盛,2002)，聲波時(shí)差低也是識(shí)別水合物的有效標(biāo)志(Murrayet al.,2005;郭星旺和祝有海,2011).相較于巖層，天然氣水合物的密度較低，接近于冰.對(duì)于半固化的沉積層，天然氣水合物的縱波速度相對(duì)較高.因此，對(duì)于由天然氣水合物粘結(jié)碎屑顆粒混合半固化沉積巖而成的天然氣水合物儲(chǔ)層，具有低密度和高縱波速度的特點(diǎn)，這一特點(diǎn)在多數(shù)的天然氣水合物勘探中得到了驗(yàn)證.

天然氣水合物儲(chǔ)層底部以下的地層巖石孔隙中常含有較多的游離天然氣，例如在青藏高原DK-1科學(xué)鉆探試驗(yàn)井中，在鉆遇天然氣水合物樣品后繼續(xù)鉆進(jìn)的過程中遭遇異常高壓氣體.這種上覆地層含有高縱波速度的天然氣水合物，下伏地層含有游離氣的地層結(jié)構(gòu)形成一個(gè)負(fù)波阻抗界面，在海域天然氣水合物的地震剖面上表現(xiàn)為BSR(Bottom Simulating Reflector)(王秀娟等,2010).而陸域天然氣水合物由于凍土層的存在，這種阻抗上的差異表現(xiàn)得不如海域天然氣水合物明顯，因此在陸地永久凍土區(qū)的天然氣水合物勘探中，很少有清晰可辨的BSR出現(xiàn).由于天然氣水合物與沉積物膠結(jié)，會(huì)使得其儲(chǔ)層在聲學(xué)上呈現(xiàn)更好的均一響應(yīng)，在剖面上會(huì)表現(xiàn)為空白或者均勻反射.當(dāng)然，對(duì)于陸域凍土地區(qū)，冰的膠結(jié)也可以產(chǎn)生類似的效果，但是由于凍土層一般分布于天然氣水合物層之上或者與之混合，所以該特征還是可以作為預(yù)測天然氣水合物儲(chǔ)層底界面的一個(gè)標(biāo)志.

盡管阻抗上的差異不是有效識(shí)別陸域天然氣水合物的有效標(biāo)志，但是如果沉積地層中含有天然氣水合物，不論含量多少，其縱波速度均會(huì)呈現(xiàn)明顯增大的現(xiàn)象.因此，在水合物穩(wěn)定帶內(nèi)，只要存在高速異常，就可以認(rèn)為有天然氣水合物分布的可能.由于鉆井資料的匱乏，利用地震數(shù)據(jù)獲取速度場信息便成為一種行之有效的手段，關(guān)鍵是需要在地震資料處理階段開展針對(duì)性的精細(xì)速度分析工作.

綜上，天然氣水合物是一種碳?xì)浠衔?，?chǔ)層下面存在游離氣層，具有高縱波速度和低密度的特點(diǎn)，因此其必然具有一定規(guī)律的地震屬性特征.研究頻率、振幅等地震屬性對(duì)天然氣水合物儲(chǔ)層的識(shí)別具有一定的指導(dǎo)意義.

哈拉湖研究區(qū)缺少鉆井、測井等直接資料，因此首先選取臨區(qū)鉆測資料開展天然氣水合物的特征分析和研究.

南祁連地區(qū)目前只有木里坳陷探測到天然氣水合物.綜合木里坳陷的鉆探結(jié)果和解釋成果，可以得出如下幾點(diǎn)認(rèn)識(shí)：(1) 木里坳陷凍土層厚度在80～120 m；(2) 天然氣水合物儲(chǔ)層分布深度范圍為150～300 m；(3) 地質(zhì)層位為侏羅系，巖性以泥巖為主，含少量粉砂巖、細(xì)砂巖；(4) 儲(chǔ)集空間類型為裂隙型(泥巖)和孔隙型(砂巖).

木里坳陷DK8、DK9井發(fā)現(xiàn)水合物，如圖1所示.DK9井發(fā)現(xiàn)水合物層段為188～200 m，厚度12 m；井位置處插入的聲波測井曲線顯示水合物層段具有相對(duì)高速的特征，其下伏地層顯示明顯低速；在過兩口井的地震剖面上，水合物層段地震反射為弱反射特征；此外，淺層有一套較穩(wěn)定的水平狀低頻反射(圖1文本框所示)，推測其為凍土層的地震反映.

3 速度異常分析

在對(duì)哈拉湖研究區(qū)地震資料處理過程中發(fā)現(xiàn)了地下存在高速異常的情況.圖2展示了SN2測線上CMP號(hào)分別為2950和3000的兩個(gè)地震道的速度譜，可以看到明顯速度反轉(zhuǎn)的現(xiàn)象.圖3展示了同一條測線上與之相鄰的CMP號(hào)分別為3050和3100的兩個(gè)地震道的速度譜，可以看到并未出現(xiàn)速度反轉(zhuǎn)異常.對(duì)比圖2和圖3，表明工區(qū)局部區(qū)域在淺層存在高速異常導(dǎo)致速度反轉(zhuǎn).圖4為四道速度譜值的疊合顯示，更加直觀地反映了上述速度反轉(zhuǎn)異常.

經(jīng)過進(jìn)一步逐道分析統(tǒng)計(jì)，繪制了速度反轉(zhuǎn)異常區(qū)域的平面分布范圍(如圖5中粉紅色線段標(biāo)示).其中，EW1測線，未見明顯反轉(zhuǎn)；EW2測線，出現(xiàn)速度反轉(zhuǎn)的CMP范圍為980～4000；SN1測線，出現(xiàn)速度反轉(zhuǎn)的CMP范圍為900～2100；SN2測線，出現(xiàn)速度反轉(zhuǎn)的CMP范圍為1800～3000.總體看，速度反轉(zhuǎn)異常區(qū)域基本分布在圖5所示構(gòu)造圖的南部凹陷及其周邊，圖中用黑色虛線勾繪了其平面范圍.

圖1 DK8、DK9鉆井及過井地震剖面反射特征分析Fig.1 DK8 and DK9drilling and reflection characteristics analysis of seismic profile across wells

圖2 SN2線CMP2950、CMP3000道集和速度譜Fig.2 Trace gathers of CMP2950 and CMP3000 on Line SN2 and the velocity spectrum

圖3 SN2線CMP3050、CMP3100道集和速度譜Fig.3 Trace gathers of CMP3050 and CMP3100 on Line SN2and the velocity spectrum

圖4 SN2線CMP2950至CMP3100速度譜疊合對(duì)比Fig.4 Comparison of velocity spectrum of CMP2950 to CMP3100on Line SN2

將EW2、SN1、SN2三條測線速度分析得到的速度譜經(jīng)DIX公式轉(zhuǎn)化為層速度剖面(圖6)，三條速度剖面的淺部(圖中500ms附近)均分布有一個(gè)相對(duì)高速層，即速度反轉(zhuǎn)區(qū)域.從地表起算深度在200 ms，按2500 m·s-1層速度計(jì)算，埋深應(yīng)在250 m左右，與水合物的形成深度一致.

4 模型正演研究

地質(zhì)建模正演模擬是地質(zhì)體識(shí)別的基礎(chǔ)分析方法.地層中的水合物儲(chǔ)層屬于地質(zhì)異常體，它的地震反射特征以及不同厚度的水合物藏地震反射特征的差異，都可以通過模型正演分析來總結(jié)規(guī)律、尋找答案.

為了更加接近真實(shí)情況，選取速度反轉(zhuǎn)異常最為明顯的EW2線作為實(shí)驗(yàn)線，將其速度譜轉(zhuǎn)化為層速度剖面，參考該層速度剖面建立水合物儲(chǔ)層地質(zhì)模型，如圖7所示.模型中部的不規(guī)則團(tuán)塊模擬的是厚層水合物，層速度2800 m·s-1.合成記錄如圖8.由記錄可見，厚層水合物具有明顯的弱振幅空白反射或均勻反射特征.

地震反射的強(qiáng)弱特征是波阻抗差異的體現(xiàn).地質(zhì)模型中頂部松散沉積物層(藍(lán)色地層)與下伏高速凍土層(青色地層)就是一個(gè)強(qiáng)波阻抗界面，在合成記錄上表現(xiàn)為一個(gè)正極性強(qiáng)振幅同相軸(圖8藍(lán)色箭頭所示).青色地層底界面與下伏地層又構(gòu)成一個(gè)稍弱的同相軸，兩個(gè)同相軸之間是凍土與水合物形成的空白帶.水合物儲(chǔ)層(黃色)雖然是相對(duì)高速層，但因其密度屬中低密度，導(dǎo)致其與圍巖的波阻抗差異較小，因而在合成記錄上表現(xiàn)為弱振幅甚至是空白反射.水合物越厚，剖面上“空白”區(qū)域越大、越明顯.由此推論，實(shí)際地震剖面上的弱振幅、均勻或空白反射可能具有與此相似的成因，是水合物儲(chǔ)層的一種直觀反映.

圖5 速度反轉(zhuǎn)區(qū)域平面位置及平面范圍預(yù)測Fig.5 Plane position and plane range prediction of velocity inversion region

圖7 EW2線層速度剖面(a)及根據(jù)該剖面建立的水合物儲(chǔ)層地質(zhì)模型 (單位：m·s-1)(b)Fig.7 Velocity profile of Line EW2 and gas hydrate reservoir geological model based on this profile

圖8 EW2線模型正演記錄Fig.8 Seismic forward records of EW2 line model

圖9 EW2線偏移地震剖面上速度異常區(qū)的地震相特征Fig.9 Seismic facies characteristics of velocity anomaly area in EW2 line migration seismic profile

圖10 EW2線振幅譜梯度屬性剖面Fig.10 EW2 line amplitude spectrum gradient attribute profile

圖11 速度異常區(qū)底界沿層25Hz單頻振幅屬性Fig.11 25 Hz single frequency amplitude attribute of the bottom boundary along the velocity anomaly area

從EW2測線的地震剖面(圖9)看，速度異常區(qū)在地震剖面上表現(xiàn)為一致的地震相特征，即表層為近水平低頻反射，其下為雜亂弱反射(圖中紅色標(biāo)示位置)，再向下又顯示為清晰的反射.該均勻弱反射區(qū)即為本次預(yù)測的天然氣水合物可能的賦存層位.

5 地震屬性分析

地震屬性分析技術(shù)能夠提取隱藏在地震資料中的有用信息，廣泛運(yùn)用于油藏特征描述等方面，能夠用來定性地討論地層巖性、物性、流體成分等，在一定程度上提高了對(duì)儲(chǔ)層有利區(qū)預(yù)測的準(zhǔn)確度.

目前，廣泛應(yīng)用于油氣檢測的地震屬性有上百種，但是尚無應(yīng)用地震屬性進(jìn)行陸域天然氣水合物的預(yù)測研究.主要原因有二：第一，針對(duì)陸域水合物的地震勘探工作非常少，資料匱乏；第二，由于凍土區(qū)的影響，獲得的地震資料大多品質(zhì)較低，即使開展屬性分析也無法得到有價(jià)值的信息.本次哈拉湖地區(qū)低頻可控震源獲得地震資料信噪比較高，具備開展地震屬性分析預(yù)測天然氣水合物儲(chǔ)層的研究條件.

現(xiàn)有關(guān)于天然氣水合物地震屬性的分析，都針對(duì)海域天然氣水合物(Andreassen et al.,1995;牛濱華等,2006;Petersen et al.,2007;徐華寧等,2010).通過海域天然氣水合物的地球物理特征研究可知，天然氣水合物出現(xiàn)的層位大都存在地震反射振幅與頻率異常.因此，本文主要利用振幅和頻率屬性進(jìn)行預(yù)測研究，并與前述速度異常區(qū)進(jìn)行對(duì)比分析.

提取EW2線的振幅譜梯度屬性(圖10)，與圖9進(jìn)行對(duì)比，振幅梯度變化較快的位置(白色線框內(nèi)區(qū)域)與圖9中所示的均勻弱反射(白色線框內(nèi)紅色區(qū)域)有很好的對(duì)應(yīng)關(guān)系，與正演研究高速異常的位置也能夠很好地吻合.該區(qū)域內(nèi)紅色表示振幅快速衰減，在地震剖面上表現(xiàn)為空白帶.

將圖5的速度反轉(zhuǎn)位置作為預(yù)測的水合物儲(chǔ)層底界，并沿該層位提取25 Hz單頻振幅屬性，如圖11所示，速度反轉(zhuǎn)區(qū)域(圖中紫色線段所示)與振幅低值區(qū)域吻合良好，這也印證了水合物層位顯示空白帶及弱振幅反射地震相的推斷.

6 結(jié)論

針對(duì)陸域天然氣水合物所開展的地震研究很少，由于永久凍土層的存在，能夠獲得高品質(zhì)地震資料非常困難.本次研究對(duì)哈拉湖地區(qū)低頻可控震源獲得的高品質(zhì)地震資料進(jìn)行了精細(xì)處理解釋，以鄰區(qū)發(fā)現(xiàn)天然氣水合物的測井資料為指導(dǎo)，分析該地區(qū)天然氣水合物賦存區(qū)的地球物理特征，獲得如下結(jié)論：

(1) 開展天然氣水合物的地震特征分析對(duì)地震資料的質(zhì)量提出了更高要求，1.5 Hz低頻可控震源寬頻激發(fā)和高覆蓋觀測系統(tǒng)的采集方法在哈拉湖地區(qū)見到了明顯效果，得到了較高品質(zhì)的地震資料，特別是針對(duì)淺層天然氣水合物勘探，成像效果明顯.

(2) 速度反轉(zhuǎn)和空白反射特征可能為凍土區(qū)水合物在地震剖面上的反映，二者同時(shí)出現(xiàn)的異常區(qū)可作為天然氣水合物勘探的重要有利目標(biāo)區(qū).

(3) 哈拉湖坳陷淺部地層具備天然氣水合物成藏條件，且在該區(qū)南部凹陷淺部地層發(fā)現(xiàn)了速度反轉(zhuǎn)異常這一重要特征，推測這是厚層水合物賦存的反映，地震相和振幅屬性分析均印證了這一推測.

Andreassen K,Hart P E,Grantz A.1995.Seismic studies of a bottom simulating reflection related to gas hydrate beneath the continental margin of the Beaufort Sea.Journal of Geophysical Research Solid Earth,100(B7):12659-12673.

Collett T S.1993.Natural gas hydrates of the Prudhoe Bay and Kuparuk River area,North Slope,Alaska.AAPG Bulletin,77(5):793-812.

Craven J A,Roberts B J,Bellefleur G,et al.2012.Recent magnetotelluric measurements at the Mallik gas hydrate production research well site,Northwest Territories.Proteomics,12(12):2045-2059.

Dallimore S R,Collett T S,Taylor A E,et al.2005.Scientific results from the Mallik 2002 gas hydrate production research well program,Mackenzie Delta,Northwest Territories,Canada.Bulletin of the Geological Survey of Canada,(585),doi:10.4095/220702.

Guo X W.2011.Well logging response characteristics and evaluation of gas hydrate in Qilian mountain permafrost.[Ms thesis] (in Chinese).Beijing:Chinese Academy of Geological Sciences.

Guo X W,Zhu Y H.2011.Well logging characteristics and evaluation of hydrates in Qilian Mountain permafrost.Geological Bulletin of China (in Chinese),30(12):1868-1873,doi:1671-2552(2011)12-1868-06.

Jiang G S.2002.Exploration and development of natural gas hydrate (in Chinese).Wuhan:China University of Geosciences Press.

Li X,Xiao L Z.2013.Natural Gas Hydrates Geophysical Characteristics and Well Logging Evaluation (in Chinese).Beijing:Petroleum Industry Press.

Moridis G J,Collett T S,Dallimore S R,et al.2004.Numerical studies of gas production from sebveral CH4hydrate zones at the Mallik site,Mackenzie Delta,Canada.Journal of Petroleum Science and Engineering,43(3-4):219-238.

Murray D,Kleinberg R,Sinha B,et al.2005.Formation evaluation of gas hydrate reservoirs.∥SPWLA 46th Annual Logging Symposium.New Orleans,Louisiana:Society of Petrophysicists and Well-Log Analysts.

Niu B H,Wen P F,Wen N,et al.2006.Estimation of hydrate concentration based on AVO modeling of BSR.Chinese J.Geophys.(in Chinese),49 (1):143-152,doi:10.3321/j.issn:0001-5733.2006.01.0200001-5733(2006)01-0143-10.

Petersen C J,Papenberg C,Klaeschen D.2007.Local seismic quantification of gas hydrates and BSR characterization from multi-frequency OBS data at northern Hydrate Ridge.Earth and Planetary Science Letters,255(3-4):414-431.

Pratt R G,Hou F,Bauer K,et al.2005.Waveform tomography images of velocity and inelastic attenuation from the Mallik 2002 cross-hole seismic surveys.∥Dallimore S R,Collett T S eds.Scientifique Results from Mallik 2002 Gas Hydrate Production Research Well Program,Mackenzie Delta,Northwest Territories,Canada.GSC Bulletin 585.Bulletin-Geological Survey of Canada,2005,585:122.

Wang X J,Wu S G,Wang D W,et al.2010.The role of polygonal faults in fluid migration and gas hydrate reservoir forming in Southeast Hainan Basin.Oil Geophysical Prospecting (in Chinese),45(1):122-128,doi:10.13810/j.cnki.issn.1000-7210.2010.01.024.

Xu H N,Yang S X,Zheng X D,et al.2010.Seismic identification of gas hydrate and its distribution in Shenhu Area,South China Sea.Chinese J.Geophys.(in Chinese),53(07):1691-1698,doi:10.3969/j.issn.0001-5733.2010.07.020.

Xu M C,Liu J X,Chai M T,et al.2012.An experimental study of natural gas hydrates in the Muli region,Qinghai Province by the seismic reflection method.Geology and Exploration (in Chinese),48(6):1180-1187,doi:0495-5331(2012)06-1180-8.

Yilmaz O.,2001.Seismic Data Analysis.∥Processing,Inversion,and Interpretation of Seismic Data Vol.1 &2.Society of Exploration Geophysicists.

Yuan S Y,Wang S X,Luo C M,et al.2015.Simultaneous multitrace impedance inversion with transform-domain sparsity promotion.Geophysics,80(2):R71-R80.

Yuan S Y,Wang S X,Tian N,et al.2016.Stable inversion-based multitrace deabsorption method for spatial continuity preservation and weak signal compensation.Geophysics,81(3):V199-V212.

Zhang H T,Zhang H Q,Zhu Y H.2007.Gas hydrate investigation and research in China:Present status and progress.Geology in China (in Chinese),34(6):953-961,doi:1000-3657(2007)06-0953-09.

附中文參考文獻(xiàn)

郭星旺.2011.祁連山凍土區(qū)天然氣水合物測井響應(yīng)特征及評(píng)價(jià)[碩士論文].北京:中國地質(zhì)科學(xué)院碩士學(xué)位論文.

郭星旺,祝有海.2011.祁連山凍土區(qū)DK-1鉆孔天然氣水合物測井響應(yīng)特征和評(píng)價(jià).地質(zhì)通報(bào),30(12):1868-1873.

蔣國盛.2002.天然氣水合物的勘探與開發(fā).武漢:中國地質(zhì)大學(xué)出版社.

李新,肖立志.2013.天然氣水合物的地球物理特征與測井評(píng)價(jià).北京:石油工業(yè)出版社.

牛濱華,文鵬飛,溫寧等.2006.基于BSR的AVO正演估算水合物含量方法的研究.地球物理學(xué)報(bào),49(1):143-152,doi:10.3321/j.issn:0001-5733.2006.01.020.

王秀娟,吳時(shí)國,王大偉等.2010.瓊東南盆地多邊形斷層在流體運(yùn)移和天然氣水合物成藏中的作用.石油地球物理勘探,45(1):122-128,doi:10.13810/j.cnki.issn.1000-7210.2010.01.024.徐華寧,楊勝雄,鄭曉東等.2010.南中國海神狐海域天然氣水合物地震識(shí)別及分布特征.地球物理學(xué)報(bào),53(07):1691-1698,doi:10.3969/j.issn.0001-5733.2010.07.020.

徐明才,劉建勛,柴銘濤等.2012.青海木里地區(qū)天然氣水合物反射地震試驗(yàn)研究.地質(zhì)與勘探,48(06):1180-1187,doi:0495-5331(2012)06-1180-8.

張洪濤,張海啟,祝有海.2007.中國天然氣水合物調(diào)查研究現(xiàn)狀及其進(jìn)展.中國地質(zhì),34(6):953-961.

(本文編輯 汪海英)

Geophysical response characteristics of gas hydrate in Har Hu region based on low-frequency vibroseis data

TIAN Yu-Kun1,LI Gui-Lin1,LIU Hui1*,ZHU You-Hai1,MA Yan-Yan1,LI Juan1,ZHOU Hui1,WANG De-Run2,WANG Lin3

1 The Key Laboratory of Unconventional Petroleum Geology,CGS,Beijing 100029,China2 Tuha Division,BGP Inc.,CNPC,Hami Xinjiang 839009,China3 BGP International，CNPC,Zhuozhou Hebei 072750,China

In China,terrestrial gas hydrate is mainly distributed in Qinghai-Tibet Plateau and Mohe County regions.Because of the existence of permafrost,it is difficult to obtain high quality data by seismic exploration,which brings many difficulties to the gas hydrate exploration.In order to solve the problem of attenuation of seismic signal by frozen soil,we used the low-frequency vibroseis as the means to do seismic data acquisition test,and obtained higher signal to noise ratio seismic data by increasing the folds.Based on the high quality seismic data,fine velocity analysis is carried out to obtain the accurate stack velocity spectrum data.And then forward modeling is established on the basis of interval velocity,seismic forward research of gas hydrate is developed.Finally,seismic attribute analysis is carried out on the post stack migration profile.The geophysical response characteristics of natural gas hydrate are summarized through the comprehensive study of the forward and the post stack seismic attributes.The mutational speed and the blank reflection zone can be used as the sensitive factors to identify the terrestrial gas hydrate in Har Hu region.

Gas hydrate；Low frequency vibroseis；Frozen soil region；P-wave velocity；Seismic attribute;Har Hu region

田玉昆，李桂林，劉暉等.2016.哈拉湖地區(qū)低頻可控震源天然氣水合物地球物理響應(yīng)特征研究.地球物理學(xué)報(bào)，59(11):4287-4296,

10.6038/cjg20161129.

Tian Y K,Li G L,Liu H,et al.2016.Geophysical response characteristics of gas hydrate in Har Hu region based on low-frequency vibroseis data.Chinese J.Geophys.(in Chinese),59(11):4287-4296,doi:10.6038/cjg20161129.

地質(zhì)礦產(chǎn)調(diào)查評(píng)價(jià)專項(xiàng)(DD20160181，12120115006601)，天然氣水合物資源勘查與試采工程(GZH201400302)和國家自然科學(xué)基金(41202076)聯(lián)合資助．

田玉昆，男，1984年生，博士，工程師，主要研究方向?yàn)榈厍蛭锢碣Y料處理反演和儲(chǔ)層預(yù)測．E-mail：tyksango@gmail.com

*通訊作者 劉暉，男，1981年生，副研究員，2011年獲得中國地質(zhì)大學(xué)(北京)博士學(xué)位，研究方向?yàn)樘烊粴馑衔锱c油氣地質(zhì)．E-mail：lhpplu@aliyun.com

10.6038/cjg20161129

P631

2016-08-02，2016-09-13收修定稿