震后快速鉆探:歷史、現(xiàn)狀與未來

Emily E.Brodsky Kuo-FongMa Jim Mori Demian M.Saffer and the participants of the ICDP/SCEC International Workshop on Rap id Response FaultDrilling

引言

大地震發(fā)生后立即在斷層帶上進行鉆探可以幫助我們獲取更多的地震信息,尤其是可獲得決定斷層動態(tài)破裂的摩擦水平和強度,觀測斷層的愈合過程及可能觸發(fā)余震的應(yīng)力變化,并可獲取控制破裂過程的重要的物理和化學(xué)屬性。在 2008年 11月由國際大陸科學(xué)鉆探項目 (ICDP)和南加州地震中心

(SCEC)聯(lián)合在日本東京舉辦的為期 3天的“震后快速鉆探:歷史、現(xiàn)狀與未來”研討會上,共有代表 10個國家的 44位科學(xué)家針對震后快速鉆探的科學(xué)和技術(shù)問題進行了深入的研討。會議梳理了震后快速鉆探可能解決的科學(xué)問題,主要議題包括震后快速鉆探和斷層帶鉆探所取得的主要成果、先進的經(jīng)驗技術(shù)以及觀測策略,詳細探討了如何權(quán)衡震后快速鉆探與滿足多領(lǐng)域科學(xué)目標(biāo)的能力之間的關(guān)系,并對適宜進行鉆探的場所進行了評估。本文概要性地介紹了會議探討的重點科學(xué)問題、觀測策略和建議,供大地震后快速鉆探項目立項參考,詳細的報告參見http:∥www.pm c.ucsc.edu/～rap id/。

1 科學(xué)動機

科學(xué)家一直在研究地震的物理過程以期能在一定程度上減輕地震災(zāi)害。目前主要的研究成果包括大致勾畫出地震可能發(fā)生的區(qū)域,用地震波測定地震的大小和地點,模擬理想狀態(tài)下的斷層破裂和地面震動,還可設(shè)計出足以抵御近震的房屋等,這些成果已經(jīng)為保護地震災(zāi)害頻發(fā)區(qū)域中生存的上億人作出了直接的貢獻。

但是,要發(fā)揮更有效的減災(zāi)作用,還需要許多與地下深部斷層的受力情況和演化進程相關(guān)的基礎(chǔ)數(shù)據(jù)。只有知道了斷層滑動過程中的應(yīng)力分布特征,才有可能預(yù)測斷層的破裂時間和滑動量。對于地震后斷層如何演化以及變化后的斷層是否會再繼續(xù)造成地震等問題我們尚一無所知,甚至連最有預(yù)報前景的余震,我們也不了解其機理?？茖W(xué)家也一直在努力尋求一些特殊的地質(zhì)特征來將孕震斷層與通過簡單蠕滑作用釋放應(yīng)力的斷層區(qū)分開來。大地震發(fā)生后立即對斷裂進行監(jiān)測,為填補這些認(rèn)識上的空白提供了有利的機會。過去的兩年中,我們已經(jīng)獲取了斷層變化的信息,通過斷層深鉆也推進了對地震物理特征的認(rèn)識。具體來說,深鉆可能解決如下問題:

1.1 是什么控制了地震的大小?

當(dāng)區(qū)域應(yīng)力超過斷層本身的摩擦強度時,就會有地震發(fā)生,而當(dāng)斷層的摩擦強度普遍較低或斷層大范圍內(nèi)積累的應(yīng)力較高時,則有大地震發(fā)生。相對而言,區(qū)域累積應(yīng)變比較容易被監(jiān)測到,而斷層的強度及其隨時間和滑動的變化特征則無法被直接觀測到。大地震往往伴隨著較大的位移和較高的滑動速度,雖然目前對斷層的摩擦強度及其變化知之甚少,但大地震發(fā)生時斷層的弱化現(xiàn)象卻非常顯著。要認(rèn)識斷層的運動過程,必然首先要了解斷層摩擦強度的變化及其對地震的控制機理。

1.2 斷層是如何愈合并開始孕育下一次地震的?

在兩次地震之間,斷層的強度恢復(fù)得非常緩慢。地震后,雖然絕大部分的斷層恢復(fù)過程尚不可監(jiān)測,但最近的研究表明,地震波速、地表變形、滲透率以及余震速率等這些間接參數(shù)均發(fā)生了快速的變化[1-2]。實驗室中則記錄到滑動之后出現(xiàn)的與速率-狀態(tài)相關(guān)的斷層愈合過程以及物理化學(xué)過程。盡管實驗室的斷層尺度與天然狀態(tài)有很大的差距[3-4],但直接的野外觀測資料將有助于進行更深入的實驗室模擬。

1.3 地震之間是如何相互作用的?

現(xiàn)有的數(shù)據(jù)顯示,近 80%的地震是被另外的地震事件觸發(fā)的[5],甚至有一部分的后續(xù)地震不論是震級還是破壞力都比前一個地震大。雖然余震是所有地震中最可能被預(yù)報出來的地震事件,但我們對于可激發(fā)下一個地震究竟是地震波的動態(tài)應(yīng)力還是靜態(tài)應(yīng)力尚一無所知。應(yīng)力變化是如何在主震之后幾周甚至幾個月引發(fā)一個新的地震的?認(rèn)識地震觸發(fā)機理是實現(xiàn)地震預(yù)報的重要一步。

1.4 發(fā)震斷層有哪些物性和化學(xué)特征?

最近在活斷層上的鉆探表明,流體、礦物和變形構(gòu)造均可以影響斷層的行為[3,6-7],某些特殊的物性組合可能會導(dǎo)致斷層發(fā)生周期性的快速滑動,而非漸變性的蠕動[8]。在斷層上取樣并測量,將有助于分析是何種物性的組合最終導(dǎo)致了大破裂的發(fā)生。而如在震后的斷層愈合過程中發(fā)生了地震,新的表面破裂區(qū)可能發(fā)生快速化學(xué)作用,進而破壞已有的滑動痕跡,因此在破裂帶上的取樣越快越好。

2 解決方案:觀測什么?為什么?

為了研究斷層的摩擦、愈合、相互作用以及成核過程,需要對熱、應(yīng)力、地質(zhì)構(gòu)造、摩擦特征、滲透性、地震波速、流體化學(xué)以及區(qū)域地面運動進行觀測。

2.1 溫度測量

穿過斷層做若干個溫度剖面是定量研究同震摩擦的最直接方法[9]。由于大部分的摩擦阻力轉(zhuǎn)換為熱,因此地震時斷層上任何的溫度增加都可以用于斷層滑動中摩擦生熱的累積量測量。為測定地震的摩擦水平,應(yīng)在地震后盡快測量發(fā)震斷層的溫度,而斷層深部剪應(yīng)力比較大、足以產(chǎn)生可觀測到的溫度異常的區(qū)域是進行溫度測量的首選地點,其中剪應(yīng)力是有效正應(yīng)力與有效摩擦系數(shù)的函數(shù) (圖 1)。在進行鉆孔溫度數(shù)據(jù)分析時,還必須考慮震后流體的平流作用對溫度觀測的影響[10]。Fulton等[11]利用多個二維數(shù)值模型模擬了瞬時流體與熱流,結(jié)果顯示,熱增壓或剪切擠壓等作用導(dǎo)致斷層局部孔隙壓力增大,由此被擠出斷層的流體直接影響了斷層上的溫度變化。這些模擬結(jié)果表明,如果真的存在高滲透等平流作用,那么盡管摩擦熱異常的形態(tài)可能受到影響,這種異常也將在兩年左右的時間中消失殆盡。另外,采用重復(fù)測井的方式監(jiān)視熱異常隨時間的變化則有可能將摩擦熱信號與因水壓變化導(dǎo)致的平流作用分開。

2.2 直接應(yīng)力測量

直接觀測斷層的絕對應(yīng)力可以定量地給出斷層的強度。震后快速鉆探時通過在不同深度獲取三維應(yīng)力的方向與大小則可以給出斷層破裂時的應(yīng)力變化。盡管估算三維應(yīng)力的大小和方向比較困難,但通過密度測井可以估算垂向應(yīng)力,通過泄漏試驗、小型預(yù)壓裂測試及水壓致裂試驗等可以估算最小水平應(yīng)力,通過斷裂寬度可以估算最大水平應(yīng)力[12]。利用非彈性應(yīng)變恢復(fù)法 (ASR)分析巖芯也可給出互補的應(yīng)力數(shù)據(jù)[13]。震后鉆探得越晚,就越難獲得應(yīng)力數(shù)據(jù)。

圖 1 當(dāng)擴散率α取 10-6 m 2·s-1時滑動量為 5 m的逆沖型地震產(chǎn)生的摩擦溫度異常[11]。(a)鉆孔深度為 1 km時的溫度異常。圖中曲線是有效摩擦系數(shù)為 0.6和0.1時的摩擦熱,分別代表強度較大和較小的斷層。實線和虛線分別表示震后 1年和 2年的摩擦熱異常。垂直虛線為假設(shè)可以被探測到的最低溫度變化,即0.2℃。(b)鉆孔深度為 2 km時的溫度異常。(c)曲線為可以觀測到 0.2℃溫度異常的最小鉆孔深度,頂部直線表示不同斷裂上鉆孔的深度和完鉆時間

2.3 巖芯地質(zhì)分析

巖芯的微尺度構(gòu)造中保留了斷層活動與摩擦損耗的詳細信息,例如假玄武巖玻璃就是斷層滑動初期具有高摩擦系數(shù)的直接證據(jù),因為這類物質(zhì)是巖石在高溫下熔融產(chǎn)生的。根據(jù)熔融組份和母巖的滲透率就可以推斷地震中被熔融物質(zhì)浸染的斷層面具有較低的摩擦系數(shù)。巖芯中微尺度的構(gòu)造還可用于認(rèn)識斷層泥流化的動力學(xué)特征以及流變特征[14]。如果在巖芯發(fā)現(xiàn)了假玄武巖玻璃或流化痕跡,就需要通過實驗手段分析熔融物的流變特征或粒子流。顆粒的尺度也包含了擴散過程的關(guān)鍵信息。在物質(zhì)開裂和表面生成過程中吸收的能量不會因摩擦等被消耗掉。用地層微電阻率掃描成像測井 (FMI)和鉆頭電阻率測井 (RAB)等井內(nèi)成像方法測定鉆孔和巖芯中裂隙的密度,則可以估算包括主滑動面在內(nèi)的能量損耗。因此,巖芯中的微構(gòu)造反映了斷層的摩擦特征。斷層帶上的礦物會在斷層愈合過程中發(fā)生快速的變化,因而需要在震后盡快取樣 (圖 2)。

圖 2 科林斯巖芯滑動面上的擦痕

2.4 實驗室測量

多數(shù)的斷層摩擦工作是基于室內(nèi)巖石有限速度滑動實驗進行的。但是天然環(huán)境中斷層的礦物特性、顆粒分布、水文環(huán)境以及破裂構(gòu)造等都使區(qū)域摩擦特性變得非常復(fù)雜。在實驗室可通過精確控制有效應(yīng)力和滑動條件來模擬實驗材料的可觀測應(yīng)力,此外,直接的實驗室摩擦測量還可采用大型巖體來逼近真實的斷層環(huán)境。為了解在地震中何種物理機制具有相對較強的控制能力,必須采集到帶有摩擦痕跡的斷層泥,并將實驗室結(jié)果與溫度測量和巖芯應(yīng)力測量結(jié)果對比。巖芯斷層泥高速實驗顯示斷層在滑動過程中出現(xiàn)了明顯的弱化現(xiàn)象[15]。

2.5 地震波速測量

地震后斷層上的波速會有所降低[16]。噪聲互相關(guān)研究表明,2004年帕克菲爾德地震之后波速的變化與區(qū)域應(yīng)變的緩慢衰減有關(guān)[1]。但目前尚不知道地震波速的變化僅僅是淺表現(xiàn)象還是由于斷層流變引起的深部過程,因而深井觀測尤為重要。而震后時間越長,這種波速變化就越不明顯。在斷層上近場觀測地殼深部的垂向地震波速變化(VSPs)可能會捕捉到斷層的愈合過程,并得到強余震期內(nèi)的彈性模量變化。在進行鉆孔觀測的同時,還應(yīng)該在地表布設(shè)短周期地震儀并定期復(fù)測。

2.6 水文地質(zhì)分析

熱或多孔彈性過程會導(dǎo)致斷層上的孔隙流體壓增高,而動態(tài)斷層弱化機制與孔隙壓的瞬時變化密切相關(guān)。斷層滲透性測量將有助于認(rèn)識斷層的愈合進程。斷層中流體的受禁程度可以用于估算地震時的剪切應(yīng)力以及下一次地震前的斷層閉鎖強度。另外,斷層的滲透性與斷層的破裂程度直接相關(guān)。

地震后斷層的滲透性很快就發(fā)生了變化[2]?？紫秹涸黾涌梢允箶鄬佑行д龖?yīng)力減小、剪切強度降低,因此孔隙壓的測量結(jié)果可用于估算斷層的靜摩擦強度[17]。圣安德烈斯斷層深部觀測鉆探 (SAFOD)[18]和臺灣車籠埔斷裂[19]鉆探顯示,鉆孔穿過的斷層帶并沒有處于過壓狀態(tài),這表明用孔隙壓的長期衰減無法解釋這些斷層帶上觀測到的弱化現(xiàn)象。

2.7 流體化學(xué)分析

化學(xué)沉淀作用會導(dǎo)致裂隙消失和粘連,因此,參照破碎帶和斷層核心區(qū)域中角礫巖給出的時間序列,通過流體化學(xué)采樣可以為斷層的愈合過程建立物理-化學(xué)模型?，F(xiàn)場還應(yīng)進行實時的氣體和水文監(jiān)測。

3 以往震后快速鉆探項目帶來的啟示

震后快速鉆探工作已經(jīng)在日本的野島斷裂、臺灣地區(qū)的車籠埔斷裂進行過。未來還會有相關(guān)的項目不斷開展。此次會議的全文版報告中還提及了一些非地震后的斷層鉆探項目,如圣安德烈斯斷層深部觀測鉆探(SAFOD)項目、科林斯斷裂的埃永 (A igion)正斷層鉆探項目及南非地震天然實驗室實施的鉆探項目。

3.1 野島斷裂鉆探項目

1995年日本神戶地震 (MW6.9)后,全球首個震后快速鉆探項目在野島斷裂上開始實施,鉆孔附近地震造成的地表破裂為 1～2m。震后 14個月開始陸續(xù)完成了 7口深度為 500～1800m的鉆井。同時進行了標(biāo)準(zhǔn)測井、鉆孔電視 (BHTV)、成像測井 (FMI)、偶極橫波聲波成像 (DSI)等地球物理探測工作。在識別斷層帶和認(rèn)識其他鉆井物理性質(zhì)時還使用了電阻率、地震波速及多種成像技術(shù)。測井結(jié)果與斷層泥、碎屑巖等巖芯地質(zhì)分析結(jié)果吻合[20]。標(biāo)準(zhǔn)測井過程中進行的溫度觀測精度較低,沒有捕捉到與摩擦熱有關(guān)的信號,但斷層附近發(fā)現(xiàn)了與摩擦熱有關(guān)的現(xiàn)象。

在斷裂帶上的花崗巖層連續(xù)取芯后發(fā)現(xiàn)其地質(zhì)結(jié)構(gòu)復(fù)雜。同時還發(fā)現(xiàn)地震活動周期與強烈的熱液蝕變有相伴關(guān)系[6,21-22],在地表 10 km以下深度形成的超級碎屑巖和假玄武玻璃與歷史上的 6～7級地震有關(guān)[23-24]。通過對巖芯進行水壓致裂試驗和應(yīng)力測量,得到了區(qū)域應(yīng)力場的方向[25-27],結(jié)果給出最大壓應(yīng)力方向與斷層垂直,這說明斷層的靜摩擦系數(shù)較低。

為研究斷層滲透率的變化,還在1800m的鉆孔中重復(fù)進行了注水試驗。從地震發(fā)生后一年開始,多年多次對斷層注水以誘發(fā)小地震。誘發(fā)地震的位置和時間以及電阻率變化表明地震發(fā)生兩年后該地點的滲透率為10-16m2～10-14m2[28],此后三年,滲透率則減少了至少 50%[2]。

3.2 臺灣地區(qū)車籠埔斷裂鉆探項目

1999年臺灣地區(qū)發(fā)生的集集MW7.6地震在車籠埔斷裂造成了較大的地表破裂,地表可見破裂達 100 km,位移為 1～12m。較淺的斷層滑動為開展孕震斷裂鉆探提供了良好的機會。在地表破裂 12m的地方,臺灣車籠埔斷裂鉆探項目 (TCDP)于震后 65個月和73個月在逆沖斷層上鉆探完成了相距 40m、深度約 1 km的 2口鉆井。此前,震后 18個月時,在斷層的北側(cè)和南側(cè)分別完成了330m和 180m的淺鉆井各 1口[29]。

TCDP項目選定地震中地表破裂和地表滑動速率較大的大坑鎮(zhèn)研究地震能量補償和斷層滑動弱化現(xiàn)象[30]。對斷層巖芯的綜合分析顯示,此次地震造成的滑動位于 A號鉆孔 1111m和 B號鉆孔 1136m深處[7,31-35]。鉆探中發(fā)現(xiàn)的黑色斷層泥從其粘土礦物組成、具有明顯的各向異性結(jié)構(gòu)以及不含有更新的微構(gòu)造等方面可以確認(rèn)是在此次地震中斷層面滑動生成的,而其中所含的細小顆粒則表明此地區(qū)在地震中釋放了較多的破裂能[7]。

在滑動量較大的斷層北側(cè)對深井和淺井均進行了溫度測量[10,29],均發(fā)現(xiàn)了斷層附近有較低水平的殘余熱,推測局部摩擦系數(shù)較低,約為 0.1～0.2(圖 3)。在溫度觀測受到了水流和熱擴散的影響,但這也從另外一個角度說明了震后快速鉆探測溫的重要性。

圖 3 臺灣鉆孔中觀測到的溫度異常及其與最大摩擦系數(shù)的關(guān)系[10]

H irono等[36]采用非破壞性方法研究了巖芯的物理特征。鉆孔A中發(fā)現(xiàn)的主滑動帶上出現(xiàn)了粘土碎屑以及異常的粒子尺度分布,這說明摩擦產(chǎn)生的熱量以及熱增壓造成了斷層泥液化[37-38]。磁化系數(shù)異常[39-40]、粘土礦物組份[31,41]、巖石化學(xué)特性[42-43]和斷層泥注入[44,38]等分析均證實了鉆探中發(fā)現(xiàn)的主滑動帶是集集地震造成的,和歷史地震產(chǎn)生的滑動帶一樣,局部溫度曾達到 900℃,且均出現(xiàn)了熱增壓和液化現(xiàn)象[42]。流體的循環(huán)似乎只局限在車籠埔斷裂上,但在地震帶和破壞區(qū)中起到控制作用的是分解-擴散-沉積過程產(chǎn)生的熱變化和同震變形[38]。

跨鉆孔的注水實驗顯示,斷層的滲透率為 10-16m2～10-18m2,實時監(jiān)測數(shù)據(jù)表明斷層處于過壓狀態(tài)[19]。淺孔巖芯樣品的實驗室分析則表明,斷層北側(cè)的滲透率比南側(cè)低一個數(shù)量級[15]。

TCDP項目在采用軟巖石進行摩擦實驗并進一步分析滑動弱化過程等方面廣泛地引起了國際地學(xué)界的關(guān)注。研究表明,主滑動帶構(gòu)造及地殼厚度有明顯的橫向變化,而斷層泥中的粘土和礦物在同震滑動的弱化過程中發(fā)揮了重要的作用。成像測井 (FMI)和偶極橫波聲波成像 (DSI)測井也顯示,主滑動帶附近應(yīng)力方向發(fā)生了變化,這表明最大水平主應(yīng)力與最小水平主應(yīng)力方向交換時斷層上應(yīng)力降最大[45]。

3.3 中國汶川斷裂科學(xué)鉆探項目

本次研討會期間,也就是 2008年 5月12日中國四川發(fā)生 MW7.9級地震后的第178天,汶川斷裂科學(xué)鉆探項目即付實施。這是目前震后鉆探實施最快的一個項目,計劃在汶川地區(qū)穿過斷層帶完成 4口深鉆。

4 地震鉆探項目推薦方案

震后鉆探關(guān)鍵之處在于快。本次會議力圖制定一個類似藍圖的規(guī)劃以供將來的快速鉆探項目參照。所有的參與者均認(rèn)識到現(xiàn)實工作中可能會遇到各種難題,且沒有發(fā)生地震時無論什么樣的設(shè)計均無法涵蓋可能出現(xiàn)的細節(jié)。但仔細探討震后快速鉆探可能解決哪些科學(xué)問題,或許可以為震后快速實施鉆探工程提供一些有益的參考。

震后快速鉆探項目適于在造成地表 1m以上滑動的大陸型地震后進行,以期能在鉆探可及的深度上獲取因較大滑動而造成的地球物理和地質(zhì)異常。具體多大的地震后可實施快速鉆探項目視科學(xué)目標(biāo)而定。在大陸上每隔 2～3年就會發(fā)生造成地表滑動超過1m的 7級以上地震,例如 1992年的蘭德斯地震、1999年的伊茲米特地震、1999年赫克托礦地震、2002年德納里 (Denali)地震都比較適宜進行震后快速鉆探,當(dāng)然還有 1995年神戶地震、1999年集集地震和 2008年汶川地震。比較適宜進行震后快速鉆探的國家包括中國、伊朗、日本、新西蘭、土耳其和美國。有部分國家已經(jīng)開始在斷層帶上進行鉆探,這將大大利于未來可能進行的震后快速鉆探項目的實施。

地震鉆探項目的主要目的是研究地震能量補償和斷層帶活動過程。首先,應(yīng)該選擇已知的和相對簡單的斷層和地質(zhì)構(gòu)造,因為在單一的主斷層比在斷裂系更容易獲取地震學(xué)資料。其次,應(yīng)該在地震反演和地表變形數(shù)據(jù)給出的斷層上同震位移較大的地區(qū)鉆探,因為地表位移量大于 1m時能捕獲到摩擦信息。理想狀態(tài)下,鉆探應(yīng)盡量接近或達到地震成核的深度,并將科學(xué)問題鎖定在摩擦、增溫及地震相互作用等方面。第三,盡量獲取結(jié)晶巖石,因為與沉積巖相比,結(jié)晶巖的滲透率較低且?guī)в懈嗟哪Σ列畔?另外,松散沉積物也不利于對孕震斷層進行深入研究。第四,應(yīng)盡量選擇傾向斷層,以便垂直鉆孔可以高角度穿過斷層面。最后,鉆探地點最好選擇在有地震觀測臺和 GPS觀測臺的地區(qū)以便獲取充足的地震信息。

為有效利用井內(nèi)空間并避免相互干擾,在設(shè)計鉆探計劃時還應(yīng)對綜合鉆孔測量和井中觀測系統(tǒng)做好規(guī)劃。為在時間和空間上盡量實現(xiàn)眾多科學(xué)目標(biāo),還應(yīng)注意觀測次序,優(yōu)選觀測手段,例如盡量采用對井孔擾動較小的應(yīng)力測試和水壓實驗,或者僅在某一段上進行某種實驗。在現(xiàn)場觀測時,溫度測量、地震觀測、應(yīng)力觀測、水壓試驗、流體采樣等在時間上或空間上都可能發(fā)生沖突。

圖 4 快速鉆探方案示例

圖 4是一個現(xiàn)場如何觀測避免沖突的例子。設(shè)計中主要考慮的因素是用最小的鉆探深度獲取有意義的摩擦熱信息、探測到儀器精度可識別的信號以及與用平均剪應(yīng)力和滑動總量推算出的熱異常量級相當(dāng)?shù)臒嵝盘枴Ｉ鲜鰧W(xué)術(shù)目標(biāo)與鉆探時間總時長是需要均衡考慮的。在震后 6～12個月鉆進 2 km對于觀測熱過程非常有利,雖然鉆孔越深越好,但時間和經(jīng)費均會隨著深度增加而呈指數(shù)增加。上述深度也正好適于地球物理觀測設(shè)備的安置。圖中的案例是以孔深2 km、距地面1750m處穿過斷層而設(shè)計的。該設(shè)計的各項細節(jié)可以根據(jù)目標(biāo)深度而調(diào)整。

5 建議

此次研討會的主要建議是,大陸地區(qū)發(fā)生地表位移超過 1m的地震之后,在有條件的地區(qū)應(yīng)盡快實施鉆探工作。鉆探的優(yōu)選地點是具有明顯地質(zhì)特征的傾向斷層上的結(jié)晶巖體。鉆探應(yīng)在震后 6個月內(nèi)開始,鉆探深度至少應(yīng)達 2 km。為實現(xiàn)此目標(biāo),與會學(xué)者建議 ICDP和其他相關(guān)學(xué)術(shù)組織進行如下工作。

5.1 對 ICDP的建議

(1)ICDP應(yīng)成立斷層鉆探分委員會。

(2)形成震后快速鉆探的評價與資助機制。

(3)斷層鉆探項目中增加應(yīng)對大地震的臨時項目。

(4)在新成立的鉆探分委員會指導(dǎo)下,建立各主要斷裂的鉆探數(shù)據(jù)資料庫。相關(guān)資料包括已經(jīng)完成和正在進行的鉆孔信息、地震與大地形變測量信息及古地震信息,以及已完成井孔的地球物理測井信息。

(5)配備應(yīng)急組合工具,包括樣品庫、鉆孔設(shè)備及特殊的井下工具等。

(6)對適于深井熱環(huán)境觀測的傳感器的研發(fā)和改進予以資助。

5.2 對其他學(xué)術(shù)組織的建議

(1)各國和地區(qū)均應(yīng)制定針對區(qū)域特征的鉆探計劃。事先擁有斷層帶附近詳實的地質(zhì)學(xué)和地球物理學(xué)背景資料對于震后快速鉆探的選址非常重要。

(2)在可能發(fā)生大地震的地區(qū)開展地震學(xué)和地質(zhì)學(xué)研究,并對已有鉆孔進行測井。這些基礎(chǔ)資料將有助于震后資料分析。

致謝

本次研討會受到國際大陸科學(xué)鉆探項目、南加州地震中心、加州大學(xué)圣克魯茲分校及東京大學(xué)的大力支持,特此感謝。感謝Emily Brodsky教授為本刊提供了原稿。

[1]Brenguier F,Campillo M,Hadziioannou C,et al.Postseismic relaxation along the San Andreas fault at Parkfield from continuous seismological observations.Science,2008,321:1478-1481,doi:10.1126/science.1160943

[2]Kitagawa Y,FujimoriK,Koizum iN.Temporal change in perm eability of the rock estim ated from repeatedwater in jection experim entsnear theNo jim a fau lt in Awaji Island,Japan.Geophys.Res.Lett.,2002,29:1483,doi:10.1029/2001GL014030

[3]H ickm an S,Sibson R,B ruhn R.In troduction to special section:Mechanical invo lvem en tof fluids in faulting.J.Geophys.Res.,1995,100:12831-12840,doi:10.1029/95JB01121

[4]Yasuhara H,Marone C,Elsworth D.Faultzone restrengthening and frictionalhealing:the role of pressure solution.J.Geophys.Res.,2005,110:B06310,doi:10.1029/2004JB0033272005

[5]Marsan D,Lengline I.Extending earthquakes’reach through cascading.Science,2008,319:1076-1079,doi:10.1126/science.1148783

[6]OhtaniT,Fujimoto K,Ito H,etal.Fau lt rocksand past to recent fluid characteristics from the borehole survey of the Nojima fault ruptured in the 1995 Kobe earthquake,southwest Japan.J.Geophys.Res.,2000,105(B7):16161-16171,do i:10.1029/2000JB900086

[7]Ma K F,Tanaka H,Song SR,etal.Slip zone and energeticsof a large earthquake from the Taiwan Chelungpu-fault D rilling Pro ject.Nature,2006,444:473-476,doi:10.1038/nature05253

[8]Dixon TH,Moore JC.The Seismogenic Zone of Subduction Thrust Faults.Irvington N Y(Columbia University Press),2007:680

[9]Lachenb ruch A H,Sass J H.Heat flow and energetics of the San And reas fau lt zone.J.Geophys.Res.,1980,85:6185-6222

[10]Kano Y,MoriJ,Fu jio R,etal.Heat signature on the Chelungpu fault associated with the 1999 Chi-Chi,Tai-wan earthquake.Geophys.Res.Lett.,2006,33:L14306,doi:10.1029/2006GL026733

[11]Fulton PM,Saffer D M,B rodsky E E.The role of advection on fau lt zone temperature after an earthquake:imp lications for rap id response d rilling.EOS Trans.AGU,Fall Meeting Supp l.,2008,89(52),Abstract T53F-05

[12]ZobackMD,Barton C A,B rudyM,et al.Determination of stressorientation andm agnitude in deep wells.Int.J.Rock Mech Min.Sci.,2003,40:1049-1076,doi:10.1016/j.ijrmm s.2003.07.001

[13]L in W,Kwasniew skiM,Imam ura T,et al.Determination of three-dim ensional in situ stresses from anelastic strain recovery measurement of cores atgreat dep th.Tectonophysics,2006,426:221-238,doi:10.1016/j.tecto.2006.02.019

[14]B rodsky E E,Rowe CD,Manighenni F,etal.A geo logical fingerp rintof extrem ely low viscosity fault fluids during an earthquake.J.Geophys.Res.,2009,114:B01303,do i:10.1029/2008JB005633

[15]Tanikaw a W,Shim amoto T.Frictional and transportp ropertiesof the Chelungpu fault from shallow borehole data and their correlationw ith seism ic behavior during the 1999 Chi-Chiearthquake.J.Geophys.Res.,2009,114:B01402,doi:10.1029/2008JB005750

[16]L i Y G,V idale J E,A k iK,etal.Evidence of shallow fau ltzone strengthening after the 1992M7.5 Landers,Califo rnia,earthquake.Science,1998,279(5348):217-219,doi:10.1126/science.279.5348.217

[17]R ice JR.Fault stress states,pore p ressure distributions,and the weakness of the San And reasfault.In:Evans B, Wong T F,(Eds.),FaultMechanicsand Transport Properties of Rocks.San Diego,Calif.,(A cademic Press),1992:475-503

[18]ZobackMD,H ickm an SH.Prelim inary observationsof stress and fluid p ressu re in and near the San Andreas faultat depth in the SAFOD boreholes.EOS Trans.AGU,FallMeet.Supp l.,2005,86(52),Abstract T21A-0438

[19]DoanML,B rodsky E E,Kano Y,etal.In situm easurem entof the hyd raulic diffusivity of the active Chelungpu Fau lt,Taiwan.Geophys.Res.Lett.,2006,33:L16317,doi:10.1029/2006GL026889

[20]Ikeda R.Outline of the fault zone drilling p ro jectby N IEDin the vicinity of the 1995 Hyogo-ken Nanbu earthquake,Japan.Island A rc,2001,10:199-205,do i:10.1046/j.1440-1738.2001.00318.x

[21]Bou llierAM,Fu jimoto K,Ito H,etal.Structuralevo lution of theNo jim a fau lt(Awaji Island,Japan)revisited from the GSJ drill hole atH irabayashi.Earth,Planets Space,2004,56(12):1233-1240

[22]Tanaka H,Om ura K,Matsuda T,etal.A rchitecturalevo lution of theNo jima faultand identification of the activated slip layer by Kobe earthquake.J.Geophys.Res.,2007,112:B07304,do i:10.1029/2005JB003977[23]BoullierAM,OhtaniT,Fujimoto K,etal.Fluid inclusions in p seudotachylytes from the Nojim afau lt,Japan.J.Geophys.Res.So lid Earth,2001,106(B10):21965-21977,doi:10.1029/2000JB000043

[24]Murakam iM,Tagam iT.Dating pseudotachylyteof theNo jim a fau ltusing the zircon fission-trackm ethod.Geophys.Res.Lett.,2004,31:L12604,doi:10.1029/2004GL020211

[25]Ikeda R,Iio Y,Om uraK.In situ stressm easurem ents in N IED bo reho les in and around the faultzone near the 1995 Hyogo-ken Nanbu earthquake,Japan. Island A rc,2001,10:252-260,doi:10.1046/j.1440-1738.2001.00323.x

[26]Tsukahara H,Ikeda R,Yam amoto K.In situ stressmeasurem ents in a boreho le close to theNojim a Fau lt.Island A rc,2001,10:261-265,doi:10.1046/j.1440-1738.2001.00324.x

[27]Yam amoto K,Yabe Y.Stresses at sites c lose to the No jim a faultm easured from core samp les.Island A rc,2001,10:266-281,doi:10.1046/j.1440-1738.2001.00325.x

[28]Tadokoro K.Repeated water injection experim ents at the No jim a fau lt zone,Japan:induced earthquakes and temporal change in perm eability.Posterp resentation.ICDP Workshop on Drilling the Chelungpu Fault.Investigating the physicsof faulting for a recent large earthquake,27-30 September2001,Taiwan,Taipei

[29]Tanaka H,Chen W M, W ang C Y,et al.Frictional heat from fau lting of the 1999 Chi-Chi,Taiwan,earthquake.Geophys.Res.Lett.,2006,33:L16316,doi:10.1029/2006GL026673

[30]Ma K F,B rodsky E E,Mori J,et al.Evidence for fault lubrication during the 1999 Chi-Chi,Taiwan earthquake(MW7.6).Geophys.Res.Lett.,2003,30(5):1244,doi:10.1029/2002GL015380

[31]Kuo L W,Song SR,Chen H Y.Characteristicsof claym ineralogy in the fau ltzoneof the TCDP and its imp lication.EOS Trans.AGU,FallMeet.Supp l.,2005,86(52),Abstract T43D-05

[32]Hung JH, Wu Y H,Yeh EC,etal.Subsurface structure,physicalp roperties,and faultzone characteristics in the scientific drill ho lesof Taiwan Chelungpu-faultD rilling Project.Terr.A tmos.Ocean.Sci.,2007,18(2):271-293,doi:10.3319/TAO.2007.18.2.271(TCDP)

[33]Sone H,Yeh EC,Nakaya T,etal.Mesoscop ic structuralobservationsof cores from the Chelungpu fau lt system,Taiwan Chelungpu-Fau ltD rilling Pro ject ho le-A,Taiw an.Terr.A tmos.Ocean.Sci.,2007,18(2):359-377,doi:10.3319/TAO.2007.18.2.359(TCDP)

[34]Yeh EC,Sone H,Nakaya T,etal.Core descrip tion and characteristicsof fau ltzones from Ho le-A of the Taiwan Chelungpu-Fau ltD rilling Pro ject.Terr.A tmos.Ocean.Sci.,2007,18(2):327-357,doi:10.3319/TAO.2007.18.2.327(TCDP)

[35]H irono T,SakaguchiM,O tsuk iK,etal.Characterization of slip zone associatedwith the 1999 Taiwan Chi-Chi earthquake:X-ray CT im age analysesandm icrostructural observations of the Taiwan Chelungpu fault.Tectonophysics,2008,449:63-84,doi:10.1016/j.tec to.2007.12.002

[36]H irono T,Yeh EC,L in W,etal.Nondestructive continuousphysical property measurements of core samples recovered from ho le B,Taiwan Chelungpu-Fau ltD rilling Pro ject.J.Geophys.Res.,2007,112:B07404,doi:10.1029/2006JB004738

[37]Boutareaud S,Bou llierA M,Beck P,et al.C last-c lay aggregates as new indicatorof shallow crustal seism ic slips.Geophys.Res.Abstr.,2008,10:1607-7962/gra/EGU 2008-A-04694

[38]Bou llierA M,Yeh EC,Boutareaud S,etal.Microscale anatom y of the 1999 Chi-Chiearthquake faultzone.Geochem.Geophys.Geosyst.,2009,10:Q03016,doi:10.1029/2008GC002252

[39]H irono T,L in W,Yeh EC,etal.H ighmagnetic suscep tibility of fau ltgougew ithin Taiwan Chelungpu fau lt:nondestructive continuousm easurem entsof physical and chem icalproperties in fault rocks recovered from Ho le B,TCDP.Geophys.Res.Lett.,2006,33:L15303,doi:10.1029/2006GL026133

[40]Mishim a T,H irono T,Soh W,etal.Therm alhistory estim ation of the Taiwan Chelungpu fau ltusing rock-m agneticm ethods.Geophys.Res.Lett.,2006,33:L23311,doi:10.1029/2006GL028088

[41]Hashimoto Y,TadaiO,Tanim izuM,etal.Characteristicsof ch lorites in seismogenic fau ltzones:The Taiwan Chelungpu Fau ltD rilling Pro ject(TCDP)co re sam p le.e-Earth,2008,3:1-6,http:∥www.electronic-earth.net/3/issue1.htm l

[42]H irono T,Yokoyam a T,Ham ada Y,etal.A chem icalkinetic app roach to estim ate dynam ic shear stressduring the 1999 Taiwan Chi-Chi earthquake. Geophys. Res. Lett.,2007,34:L19308, doi: 10.1029/2007GL030743

[43]Ishikawa T,Tanim izuM,NagaishiK,etal.Coseism ic fluid-rock interactionsathigh temperatures in the Chelungpu fault.Nature Geosci.,2008,1:679-683,doi:10.1038/ngeo308

[44]O tsukiK,Uduk iT,Monzaw aN,etal.C layey injection veinsand p seudo tachylyte from two bo reho lespenetrating the Chelungpu fau lt,Taiwan:their imp lications for the contrastive seism ic slip behaviors during the 1999 Chi-Chiearthquake.Island A rc,2005,14:22-36

[45] W u H Y,Ma K F,ZobackM,etal.Stressorientationsof Taiwan Chelungpu-Fau ltD rilling Pro ject(TCDP)ho le-A as observed from geophysical logs. Geophys. Res.Lett.,2007,34:L01303,doi:10.1029/2006GL028050