地球物理學

青藏高原地震活動特征及當前地震活動形勢

鄧起東，程紹平，馬冀，等

摘要：青藏高原是我國現代構造活動和地震活動最強烈的地區，在高原內已記錄到多達 19次 8級以上巨大地震和100余次7～7.9級地震，它們均發生在喜馬拉雅板塊邊界構造帶和板內斷塊區及其次級斷塊的邊界活動構造帶上。高原南緣的喜馬拉雅構造帶是印度板塊與歐亞板塊俯沖碰撞邊界，由MCT（主中央斷裂）、MBT（主邊界斷裂）和MFT（主前緣斷裂）3個巨大逆斷層組成，這是一個由北向南發展的前展式逆斷層系，其中以MFT活動性最強，是一條強烈活動的地震帶，已先后發生了9次8級以上強震，2015年尼泊爾廓爾喀（Gorkha）8級地震是其中最新的一次地震。青藏高原內部是一個由多個次級斷塊組成的斷塊區，主體包括有拉薩斷塊、羌塘斷塊、巴顏喀喇斷塊、東昆侖—柴達木斷塊和祁連山斷塊，這些走向北西西向的條狀斷塊邊界被活動斷裂所控制，以不同的速率自西向東滑動，高原內部10次8級地震和絕大多數7級和7級以上地震也都發生在斷塊邊界構造帶上。早在 1900年以前的地震記錄即已存在某種叢集起伏現象，1786—1850年和1863—1905年的地震活動即存在兩次叢集。自1900年有地震儀器記錄以來，青藏高原明顯存在 3次地震活動叢集高潮，即 1920—1937年、1947—1976年和1995年—現在，在每次地震活動叢集期都形成以8級地震為核心的7級以上地震活動系列，它們分別是20～30年代的海原—古浪地震系列，50～70年代的察隅—當雄地震系列和20世紀末期以來的昆侖—汶川地震系列。每個地震系列都有自己的主體活動區，最新的昆侖—汶川地震系列的主體活動區即為巴顏喀喇斷塊。每一個地震系列都有自己的核心地震，如海原—古浪，察隅—當雄和昆侖—汶川地震系列的核心地震分別為海原和古浪地震、察隅和當雄地震及昆侖和汶川地震。巴顏喀喇斷塊南北兩邊界上 1997年瑪尼7.5級地震，2010年玉樹7.1級地震和2001年昆侖山口西8.1級地震均為走滑型地震，斷塊東端2008年汶川8級地震和2013年蘆山7.0級地震反映了巴顏喀喇斷塊前端為擠壓邊界，后方2008年7.3級于田地震為正斷層，反映了該塊體后端為拉張邊界，2014年的于田7.3級地震由于更靠近阿爾金斷裂而具有較大的走滑分量，以上共同說明巴顏喀喇斷塊的東向滑動。自1950年以來，全球已經經歷了兩期強震叢集峰值期，第一期出現在1950—1965年，發生了3次Mw≥9.0特大地震和許多Mw≥8.0強震。第2期出現在2001年以來，已發生了2次Mw≥9.0的特大地震。正如前述，青藏高原的地震活動性與這些高水平全球大地震在時間和空間分布上的一致性，充分反映了它們兩者間的動力學聯系。目前，雖然又已經歷了數年短暫的大地震相對平靜，我們仍不能掉以輕心，仍需注意全球板塊邊界構造帶 8～9級地震和板內大陸斷塊區 7～8級地震的發生，仍需注意未來地震危險性，加強地震監測工作。

來源出版物：地球物理學報, 2014, 57(7): 2025-2042

入選年份：2016

基于不等式約束的最小二乘法三維電阻率反演及其算法優化

劉斌，李術才，李樹忱，等

摘要：目的：隨著直流電阻率法應用范圍的拓展，對三維電阻率反演方法的反演精度和反演速度等的要求越來越高，但仍存在一些未能解決的問題，主要有（1）非唯一性問題是三維電阻率反演的固有問題，有時表現得較為嚴重；（2）反演耗時較長，成為制約三維電阻率反演在實際應用中推廣的瓶頸。針對以上問題，研究了攜帶不等式約束的三維電阻率反演方程，基于Cholesky分解算法和預條件共軛梯度算法（PCG）的特點提出了反演計算效率優化方案。方法：為減小三維電阻率反演問題的多解性，將模型參數的變化范圍（實際介質電阻率取值范圍）作為先驗信息，將表征模型參數變化范圍的不等式約束施加到反演方程中，通過構建攜帶光滑約束和不等式約束的三維電阻率反演目標函數，以降低反演的多解性。針對某一具體的觀測數據而言，可以根據一般物理常識或者鉆孔等方式獲得的不等式約束范圍，也可以根據視電阻率數據的特點來確定不等式約束。針對三維電阻率反演計算效率優化的問題，在整個反演流程中，最為耗時的是敏感度矩陣的計算和反演方程的求解，提高二者的求解速度對于優化三維電阻率反演的計算效率至關重要。據此，提出了三維電阻率反演計算效率優化方案：Cholesky分解法被用來求解敏感度矩陣計算中的多個點源場的正演問題，該方法只需對總體系數矩陣進行一次分解，然后對不同的右端向量進行回代即可；預條件共軛梯度法被引入到反演方程的求解中，將雅可比迭代中的對角陣作為預處理矩陣，具有求逆方便、無需內存空間的特點，能夠加快收斂速度。結果：為了評價最小二乘法三維電阻率反演的效果，利用合成觀測數據以及煤礦突水模型試驗中巖層斷裂及突水通道形成過程電阻率監測的試驗數據進行了反演計算研究。合成算例中反演模型網格數量為 15×15×40＝9000個，模型的反演迭代次數為10次，耗時約160 min（計算機配置為主頻 1.83 GHz，內存 1.0 GB）。得到的反演結果中兩個低阻異常體在位置、規模、形狀、電阻率值等特征方面均與原模型基本一致，尤其是提取出的常體與原模型中的低阻體完全一致。煤礦突水模型試驗算例中三維反演模型在x，y，z方向上的網格數分別為70×2×20＝2800個，達到收斂所需的迭代次數為7次，總耗時為520 s（計算機主頻為1.83 GHz，內存為1.0 GB），表明反演方法具有迭代次數少，收斂速度快的優點。與實際試驗情況對比，反演結果中的高阻裂隙位置與試驗中裂隙發生的位置基本一致，表明所提出的三維電阻率反演方法為核心技術的電阻率層析成像監測系統可將裂隙擴展、巖層破斷等重要突涌水災害前兆過程以成像的方式直觀表達出來，能夠捕捉到突水前兆信息。結論：反演算例表明，借助不等式約束和反演效率優化方案，三維電阻率最小二乘反演方法可得到較為精確的反演結果，反演效率得到提高。煤礦突水模型試驗反演算例表明，三維電阻率反演方法的優點有利于直流電阻率勘探方法在地質災害實時監測這一新領域推廣。

來源出版物：地球物理學報, 2012, 55(1): 260-268

入選年份：2016