貝克盆地衛(wèi)星重力場特征與構造單元劃分

張 科，方 勇，陳全紅，襲著綱，程岳宏

0 引言

貝克盆地位于印尼東部巴布亞省貝克島和Yapen 島之間海域，南緯 0．5°～1．5°，東經(jīng) 135°～137°之間，盆地面積約 2 ×105km2［1］。該區(qū)遠離印尼主要工業(yè)城市，且水深較大，勘探成本較高，因此，多年來該區(qū)尚未開展大規(guī)模油氣勘探工作，勘探程度極低，只有稀疏的二維地震測線，尚未發(fā)現(xiàn)油氣。進入21世紀之后，隨著印尼東部油氣勘探的不斷開展，2007年TGS公司在該盆地陸續(xù)采集了多用戶二維地震資料，船測重力資料，并在貝克島及周緣島嶼開展了地質踏勘工作，據(jù)初步研究表明，該盆地油氣苗比較活躍，且發(fā)育較大規(guī)模的坳陷，具備較大的勘探潛力。

由于資料短缺所限，前人在該區(qū)研究僅局限于島嶼露頭的踏勘及稀疏二維地震資料的分析，盆地邊界及內部構造單元不明確。圖1中白色方框區(qū)域內為貝克盆地大致范圍，在貝克島油氣苗露頭較多，可能具有較大的勘探潛力。因此需要對盆地的分布范圍、盆地主要結構樣式和規(guī)模、主要斷裂展布進行研究。

限于目前資料情況，此次研究主要借助最新的高精度衛(wèi)星重力資料(V 18．1)，以區(qū)域地質、地震及鉆井資料為約束，充分挖掘衛(wèi)星重力資料信息與地震資料相結合，開展盆地構造區(qū)劃研究。

衛(wèi)星重力具有高精度、高覆蓋和高重復率等優(yōu)點，能給出大尺度的重力變化。目前，衛(wèi)星測高重力分辨率已經(jīng)接近10～8 km(半波長)的觀測精度，基本滿足海洋地質與資源調查精度要求［2］。世界范圍內的衛(wèi)星測高重力數(shù)據(jù)逐步被石油公司針對其研究海域進行生產(chǎn)，并應用于前期構造研究，比如盆地結構分析、構造單元劃分，較好的解決了地震資料稀缺所造成的研究困難［3－5］。同時，在缺乏船測重力的研究區(qū)可以借用衛(wèi)星重力資料進行替代，縮短研究周期，節(jié)約研究成本，達到研究目的。國內自20世紀90年代以來在海域如南黃海［6］、南海［7］和陸域的新疆天山［8］及塔里木盆地［9］等地區(qū)做了較好的探索及應用;國外有石油公司利用高精度衛(wèi)星重力資料在南亞的印度洋中央盆地［10］和柬埔寨磅遜盆地［11－12］、南美的桑托斯盆地［13］、西非的科特迪瓦盆地［14］和印尼東部 bone盆地［15］等油氣快速普查與早期評價中取得了較好的應用效果。

圖1 貝克盆地區(qū)域構造背景Fig．1 Tectonic map of Biak Basin

1 貝克盆地構造背景及演化模式

貝克盆地呈近東西走向，為受該區(qū)大型走滑斷裂帶Sorong斷裂作用形成的走滑拉分盆地。從構造位置上看，貝克盆地位于向西南俯沖的太平洋板塊與向東北漂移的澳大利亞板塊及西側的歐亞板塊共3個大型板塊匯聚區(qū)(圖1)。貝克盆地位于構造應力釋放區(qū)，主要以走滑和擠壓構造活動為主，相關學者(Valerio Memmo等，2013)提出一個在左旋走滑斷裂作用下形成演化模型。如圖2所示該模型由一個寬脆性剪切帶(Cladouhos，1999a)，海陸過渡區(qū)與北部分隔，Yapen斷裂帶與南部隔離。這些板塊尺度的構造輪廓從宏觀上較好的表現(xiàn)了左行走滑斷裂，在貝克島和南貝克島的橫向運動產(chǎn)生內部微巖層獨立向西移動［16］。

與南貝克島相對應，貝克島向西高速度運動，產(chǎn)生的貝克島開口，形成拉分盆地。在這種構造背景下，貝克島可能是一個S形擠壓結構逐步由Riedel斷裂系統(tǒng)旋轉移位形成的。該模型是符合該區(qū)域的區(qū)域構造背景以及Yapen島右旋走滑斷層(Charlton，2008)，因此可以被解釋為“多米諾式”的 Yapen 斷裂結構［17］。

圖2 貝克盆地形成模式（據(jù)Valerio Memmo，2013）Fig．2 Structural model explaining the formation of Biak Basin

2 貝克盆地布格重力場特征

此次研究首次采用的海底深度資料和衛(wèi)星重力數(shù)據(jù)來自美國2位教授 David T．Sandwell和Walter H．F．Smith共同維護的全球衛(wèi)星重力數(shù)據(jù)V 18．1，數(shù)據(jù)測網(wǎng)密度為1'×1'，其中重力數(shù)據(jù)精度可以達到 3．03×10－5m/s2，滿足區(qū)域研究精度需求。布格重力異常計算選用統(tǒng)一的中間層密度，地殼平均密度取2．30×103kg/m3，海水密度取1．03×103kg/m3，用雙界面模型重力場快速正演方法(王萬銀和潘作樞，1993)計算海水的重力影響值，進而對衛(wèi)星測高重力異常做海水影響校正，并由此計算得到布格重力異常［18－19］。

在對衛(wèi)星測高重力數(shù)據(jù)去除海水影響之后獲得研究區(qū)的布格重力異常，從圖3的等值線圖可以看出，研究區(qū)布格重力異常受南北2個島嶼夾持，呈近東西走向分布，在重力高、低過渡帶發(fā)育有重力梯級帶。布格重力異常值的大小反映了基底的起伏特征，重力異常值大，反映高密度的基底埋深相對較淺;重力異常值小，反映基底埋深相對較深，上覆地層較厚［20］。布格重力異常等值線圈閉的重力高、重力低反映局部隆起與坳陷構造，布格重力異常等值線密集的重力梯級帶反映斷裂較發(fā)育區(qū)。由研究區(qū)布格重力場值中西部發(fā)育的重力低值帶及中東部發(fā)育重力高值帶特征及分析研究區(qū)基底為中西部較深，中東部較淺;由布格重力異常等值線圈閉的重力高、重力低推測研究區(qū)基底起伏較大，隆起與坳陷發(fā)育。

圖3 貝克盆地布格重力異常Fig．3 Bouguer gravity anomaly of Biak Basin

3 貝克盆地主要斷裂展布特征

控盆斷裂的重力異常有明顯的特征，利用重力異常識別斷裂位置的方法有很很多，文中利用近年來比較先進且有效的位場邊緣識別技術獲取貝克盆地布格重力異?？偹剑瓜驅?shù)，根據(jù)其極大值位置或者極大值的錯斷位置識別相對應的斷裂。據(jù)圖4所示，貝克盆地南部的Sorong走滑斷裂帶具有明顯的異常極大值，同時由其派生南部(F1)和北部(F2)2組大型斷裂控制了盆地的南北端，西部斷裂F3為后期形成，控制盆地西端。從圖5中可以看出，近東西向斷裂F1和F2表現(xiàn)為明顯的極大值連線，并且連續(xù)性較好，異常幅值明顯，具有分段性特征。F1斷裂為Sorong斷裂Yapen島部分，主要位于海域，在海陸分界處，斷裂發(fā)生明顯錯動，極大值表現(xiàn)為明顯的錯斷，整體幅值較窄，在東南段幅值較為寬緩。F2斷裂貫穿貝克島，西北端極大值發(fā)生錯斷，異常幅值較窄，極大值向東南方向延伸連續(xù)較好，幅值逐漸變寬緩;F3斷裂呈北東－南西走向，主要位于海域，為布格重力異?？偹剑瓜驅?shù)極大值的連線，異常幅值較窄，極大值向東北方向延伸連續(xù)較差，幅值逐漸變寬緩。

4 貝克盆地構造單元劃分

剩余布格重力異常是在布格重力異常的基礎上，去除區(qū)域重力異常后的剩余部分，主要反映局部地質構造的重力異常特征。此次運用位場分離方法較好的分離出了局部密度體的重力異常，借助區(qū)域研究和地質剖面(圖5)，識別了主要大型坳陷和隆起，并結合盆地斷裂展布特征劃分了構造單元(圖6)。

圖4 貝克盆地布格重力異?？偹剑瓜驅?shù)Fig．4 NVDR_THDR of the Bouguer gravity anomaly

圖5 貝克盆地區(qū)域剖面（測線位置如圖6所示）Fig．5 Geosection of the Biak Basin

圖6 貝克盆地構造單元劃分Fig．6 Tectonic units in Biak Basin背景值為剩余布格衛(wèi)星重力異常，方框為研究區(qū)，紅線為圖5測線位置

利用剩余布格重力異常的零值線進一步明確了貝克盆地邊界(圖6中紅線所示)，根據(jù)異常所對應的隆起和坳陷劃分貝克盆地二級構造單元，整體從西到東呈“西高東低、三坳三隆、隆坳相間、東西分塊”的構造格局，分別為西部坳陷、西部隆起、中央坳陷、中南隆起、東南坳陷和東北隆起共6個二極構造單元，與區(qū)域地質剖面具有較好的對應關系。

5 結論

1)根據(jù)布格重力異常NVDR_THDR共推斷出貝克盆地一級斷裂3條(F1，F(xiàn)2，F(xiàn)3)，斷裂對貝克盆地構造單元的分布具有控制作用;

2)利用剩余衛(wèi)星測高重力異常的零值線劃分了貝克盆地邊界，根據(jù)異常相對高低識別了貝克盆地的6個構造單元，整體從西到東呈“西高東低、三坳三隆、隆坳相間、東西分塊”的構造格局;

3)研究表明貝克盆地發(fā)育大型坳陷，面積大于6 000 km2，沉積地層厚度大于7 000 m，具備較大生烴潛力，具備較大的勘探價值，為該盆地下一步研究指明了方向。

［1］ Valerio Memmo．Deposition and deformation in the recent Biak Basin［J］．Proceedings，Indonesian:Petroleum Association，2013(37):122 －131．

［2］ 焦孟梅，張景發(fā)，姜文亮，等．衛(wèi)星重力發(fā)展及應用［C］//地殼構造與地殼應力文集，2007:95－108．JIAO Meng-mei，ZHANG Jing-fa，JIANG Wen-liang，et al．Progress and application of satellite gravity measurement techniques and the theory［C］//Collection of Tectonic and Crustal Stress，2007:95 －108．

［4］ 姜 林，薄冬梅，周 波，等．塔里木盆地輪南地區(qū)油氣成藏期次分析［J］．西安科技大學學報，2015，35(4):444－449．JIANG Lin，BAO Dong-mei，ZHOU Bo，et al．Hydrocarbon accumulation periods in Lunnan Area of Tarim Basin［J］．Journal of Xi’an University of Science and Technology，2015，35(4):444 －449．

［5］ 楊 婕，占 惠．衛(wèi)星重力探測技術的發(fā)展［J］．國際地震動態(tài)，2008(5):23－27．YANG Jie，ZHAN Hui．Development of the satellite gravity surveying technology［J］．Recent Developments in World Seismology，2008(5):23 －27．

［6］ 王 佟，馮 帆，王慶偉，等．中國西北地區(qū)柴北緣賦煤區(qū)構造單元新厘定［J］．西安科技大學學報，2016，36(1):24－30．WANGTong，F(xiàn)ENGFan，WANGQing-wei，et al．A new division of the structural units of the north Qaidam coal zone in Northwest China［J］，Journal of Xi’an University of Science and Technology，2016，36(1):24 －30．

［7］ 王萬銀，劉金蘭，邱之云，等．利用衛(wèi)星重力異常研究南黃海地區(qū)中生界厚度［J］．中國海上油氣，2004，16(3):151－157．WANG Wan-yin，LIU Jin-lan，QIU Zhi-yun，et al．A research on mesozoic thickness using satellite gravity anomaly in the southern Yellow Sea［J］．China Offshore Oil and Gas，2004，16(3):151 －157．

［8］ 王懋基，宋正范，尹春霞．南海衛(wèi)星重力研究［J］．物探與化探，1998，22(5):329 －336．WANG Mao-ji，SONG Zheng-fan，YIN Chun-xia．Satellite gravity study on South China Sea［J］．Geophysical and Geochemical Exploration，1998，22(5):329 －336．

［9］ 樓 海，王椿鏞，王 飛．衛(wèi)星重力資料揭示的新疆天山地區(qū)構造動力學狀態(tài)［J］．地震學報，2000，22(5):482－490．LOU Hai，WANG Chun-yong，WANG Fei．Dynamic features of the Tianshan orogen deduced from satellite gravity data［J］．Acta Seismologica Sinica，2000，22(5):482－490．

［10］ 張 健，周國藩，徐忠祥．塔里木盆地衛(wèi)星重力異常與油氣資源的相關性研究［J］．地球科學，1997，22(6):643－647．ZHANG Jian，ZHOU Guo-fan，XU Zhong-xiang．The relativity of hydrocarbon resources and the satellite gravity anomaly in Tarim Basin［J］．Earth Science，1997，22(6):643－647．

［10］ Gopala Rao D．Satellite gravity anomalies and crustal features of the Central Indian Ocean Basin［J］．Current Science，2004，86(7):947 －956．

［11］ 張 科，呂 棟．重震聯(lián)合研究在磅迅盆地中的應用［J］．物探與化探，2013，37(6):983 －987 ZHANG Ke，LV Dong．The application of seismic and gravity integrated method to the Kampong Basin，Cambodia［J］．Geophysical and Geochemical Exploration，2013，37(6):983 －987．

［12］ 張 科，趙汝敏，李愛山，等．衛(wèi)星重力資料在柬埔寨磅遜盆地中的應用研究［J］．地球物理學進展，2014，29(1):339－345．ZHANG Ke，ZHAO Ru-min，LI Ai-shan，et al．The satellite gravity data applied in the Cambodia Kampong Basin［J］．Progress in Geophysics，2014，29(1):339－345．

［13］ Rader Abdul Fattah．Use of GOCE satellite gravity data for basin and petroleum system analysis in the Arabian Peninsula［C］//In:10thMiddle East Geosciences Conference and Exhibition，Qatar，2012:1 －18．

［14］ Fainstein R F，Roy L．Seismic constrained gravity inversion for sub-basalt exploration in West Coast，India［M］．Amsterdam:71stEAGE Conference and Exhibition，2009．

［15］ 楊 敏，王萬銀，杜向東，等．利用衛(wèi)星重力數(shù)據(jù)研究印尼Bone盆地構造區(qū)劃［J］．地球物理學進展，2016，31(3):999 －1009．YANG Min，WANG Wan-yin，DU Xiang-dong，et al．Tectonic districts by satellite gravity data of the Bone Basin，Sulawesi，Indonesia［J］．Progress in Geophysics，2016，31(3):999 －1009．

［16］ Mark Thompson，Cholldy Reminton．Detection of liquid hydrocarbon seepage in Indonesian offshore frontier basin using airborne laser fluotosensor(ALF)the result of a pertamina/BP joint study［C］//20thAnnual Convention Proceedings，1991:76 －95．

［17］ Mayall M J，Cox M．Deposition and diagenesis of miocene limestones，Senkang Basin，Sulawesi，Indonesia［J］．Sedimentary Geology，1988，59:77 －92．

［18］ 劉元龍，鄭建昌，武傳珍．利用重力資料反演三維密度界面的質面系數(shù)法［J］．地球物理學報，1987，30(2):186－196．LIU Yuan-long，ZHENG Jian-chang，WU Chuan-zhen．Inversion of three dimensional density discontinuity by use of a method of“compressed mass plane coefficient”based on the gravity data［J］．Chinese Journal of Geophysics，1987，30(2):186 －196．

［19］ 王萬銀，潘作樞．雙界面模型重力場快速正演問題［J］．石油物探，1993，32(2):81 －87．WANG Wan-yin，PAN Zuo-Shu．Fast solution of forward and inversion problems for gravity field in a dual interface model［J］．Geophysical Prospecting for Petroleum，1993，32(2):81 －87．

［20］ 王萬銀，邱之云，楊 永，等．位場邊緣識別方法研究進展［J］．地球物理學進展，2010，25(1):196 －210．WANG Wan-yin，QIU Zhi-yun，YANG Yong，et al．Some advances in the edge recongation of the potential field［J］．Progress in Geophys，2010，25(1):196－210．

［21］ 馮 娟，孟小紅，陳召曦，等．三維密度界面的正反演研究和應用［J］．地球物理學報，2014，57(6):1934－1945．FENG Juan，MENG Xiao-hong，CHEN Zhao-xi，et al．The investigation and application of three-dimensional density interface［J］．Chinese Journal of Geophysics，2014，57(6):1934 －1945．