2023年2月5日三水ML4.0級(jí)礦山地震震相特征與發(fā)震構(gòu)造

楊 選，馬曉靜

（廣東省地震局，廣州 510070）

0 引言

2023 年2 月5 日，廣東省佛山市三水區(qū)樂(lè)平鎮(zhèn)附近發(fā)生了ML4.0 級(jí)地震，中國(guó)地震臺(tái)網(wǎng)中心（CENC）正式目錄的發(fā)震時(shí)間是11∶38∶15.27，震源位置是經(jīng)度112.953°E，緯度23.279°N，深度11 km。佛山、廣州、深圳等地區(qū)有感，有感范圍廣泛，且震中位于自l988年3月開始注水采鹽的隔坑鹽礦礦區(qū)[1]，引起了廣大群眾及政府部門的高度關(guān)注。這次地震部分臺(tái)站記錄波形中出現(xiàn)規(guī)則正弦形面波，震中周邊20 km 范圍內(nèi)，震前5 年僅發(fā)生過(guò)3 次地震，分別是2020 年1 月20 日ML1.0、2022 年5 月13 日ML1.6 和2023 年1 月23 日ML0.9 地 震。震后一個(gè)月內(nèi)沒(méi)有記錄到余震活動(dòng)，地震波形及地震活動(dòng)類型均較為獨(dú)特。利用現(xiàn)有觀測(cè)資料深入剖析此次地震的特征和發(fā)震構(gòu)造背景，對(duì)于認(rèn)識(shí)本區(qū)域特殊類型地震、地震活動(dòng)性以及淺源地震災(zāi)害風(fēng)險(xiǎn)防治等具有一定意義。

1 地震記錄資料及特征

為研究這次地震的特征，本文收集了震中800 km 范圍內(nèi)廣東（包括廣州和佛山等地方臺(tái)網(wǎng)以及新建成尚未驗(yàn)收的40 個(gè)加密臺(tái)）、廣西、湖南、福建、江西和海南等臺(tái)網(wǎng)[2]的2000 多個(gè)臺(tái)站的記錄，從中挑選出質(zhì)量較高的342 個(gè)測(cè)震記錄和410 個(gè)強(qiáng)震記錄，共計(jì)752個(gè)（見圖1）。經(jīng)過(guò)人工仔細(xì)分析，識(shí)別出了967個(gè)震相（見圖2），包括Sg震相280個(gè)、Pg 震 相264 個(gè)、Pn 震 相199 個(gè)、Sn 震 相123 個(gè)、PmP震相23個(gè)、SmS震相34個(gè)和Pb震相24個(gè)。

圖2 震相走時(shí)曲線Fig.2 Travel time curve of seismic phase

震中距30 km 內(nèi)的8 個(gè)臺(tái)均有多組反射波，Sn與Sg 走時(shí)曲線交點(diǎn)對(duì)應(yīng)的震中距在150 km 左右，相對(duì)較大（據(jù)華南走時(shí)表中的數(shù)據(jù)推測(cè)：震源深度為0、5、10、15和20 km時(shí)，此交點(diǎn)對(duì)應(yīng)的震中距分別為164.5、151.7、136.4、119.2 和100.0 km），是典型的上地殼事件。

所有臺(tái)站記錄都在Sg 波之前有豐富的高頻成份，有75 個(gè)臺(tái)站記錄中在Sg 波之后出現(xiàn)了周期變化（周期在1.1～1.5s 之間）、毛刺少、振動(dòng)規(guī)則的短周期面波，其持續(xù)時(shí)間大致在5 至20 s 之間；面波振幅振幅整體隨震中距增大而減小，但部份臺(tái)存在突然增大數(shù)倍的現(xiàn)象（見附圖2），最大峰值的視速度約為2.85 km/s，在垂直向比較明顯，符合Rg 面波的特征[3]。Rg波是短周期的地殼瑞利波，因?yàn)樘烊坏卣鹩幸欢ǖ纳疃?，所以極少出現(xiàn)Rg 波。正是因?yàn)槌霈F(xiàn)Rg 波，同時(shí)又有較多的高頻成份，所以三水地震的類型爭(zhēng)議較大，圖3給出了DOG臺(tái)記錄到的本次地震與2016 年10 月28 日廣東連平M3.0塌陷的波形，兩者差異較為明顯，塌陷滑坡高頻成份不發(fā)育，爆破通常具有各向同性，S 波不發(fā)育的特征，據(jù)此可排除此次事件為塌陷、滑坡或爆破的可能。

圖3 DOG臺(tái)記錄到的本次地震與一次塌陷波形Fig.3 Waveforms of this earthquake and a collapse recorded by the DOG station

圖4 視出射角與真出射角示意圖Fig.4 Schematic diagram of apparent and true dip angles

從能觀測(cè)到Rg 波的臺(tái)站空間分布看：多數(shù)臺(tái)站分布在江河流域或出海口，這些區(qū)域沉積層相對(duì)較厚，由此推斷Rg 波的發(fā)育程度似乎還與沉積層厚度有一定關(guān)聯(lián)。

2 地震深度及精確確定位

在沒(méi)有較精細(xì)地殼模型的情況下，即使有較多的近臺(tái)記錄，深度搜索仍然還是不穩(wěn)定，特別是JOPENS-MSDP 提供的幾種定位算法沒(méi)有對(duì)近臺(tái)資料特別增加權(quán)重設(shè)置，使得精確測(cè)量震源深度更困難。但是，在有較近的臺(tái)站（震中距小于兩倍震源深度）記錄時(shí)，通過(guò)增加測(cè)量視出射角，進(jìn)行單臺(tái)射線追蹤，能對(duì)深度起到一定的約束作用，可改善對(duì)震源深度的分辨效果[4]。

2.1 視出射角法

地震波入射到地面時(shí)，總位移矢量和地面之間的夾角稱為視出射角，常用ē表示；而入射波與地表的夾角稱為真出射角，常用e表示。由于地震波在地面上形成了反射波，所以視出射角不等于入射波的真出射角。若將地球表面視為自由表面，則視出射角和真出射角的關(guān)系與地殼介質(zhì)的泊松比有關(guān)[5]。在有三分向記錄到清晰初動(dòng)的近臺(tái)記錄時(shí)，可根據(jù)記錄到的地面位移的垂直分量AUD和兩個(gè)水平分量AEW與ANS按（1）式先確定視出射角，再根據(jù)真出射角和視出射角的關(guān)系（2）式計(jì)算真出射角，測(cè)量出直達(dá)P 和S 波的到時(shí)tS和tP即可按（3）式，估計(jì)出震源深度[5-6]，式中tS-tP是直達(dá)S與直達(dá)P波的到時(shí)差，U為虛波速度，e為真出射角。

在收集到的記錄中，Pg 較早達(dá)到HWG（tP=11∶38∶16.590）和SSSE（tP=11∶38∶17.160）兩臺(tái)，Sg波的到時(shí)使用質(zhì)點(diǎn)運(yùn)動(dòng)軌跡分析，較精確地找出質(zhì)點(diǎn)運(yùn)動(dòng)軌跡明顯改變的起跳點(diǎn)（圖5），分別是GZ/HWG（tS=11∶38∶17.290）和FS/SSSE（tS=11∶38∶17.970），仿真成位移記錄測(cè)量三分向初動(dòng)，據(jù)范玉蘭等人的研究，華南地區(qū)直達(dá)P 波速度VPg=6.01 km/s，直達(dá)S 波速度VSg=3.55 km/s，虛波速度U=8.673 km/s，可估算出本次地震的震源深度在4.3 km左右，見表1。

表1 視出射角法估算本次地震震源深度結(jié)果Table 1 The focal depth results of this earthquake estimated by the apparent dip angle method

圖5 HWG臺(tái)（左）和SSSE臺(tái)（右）的質(zhì)點(diǎn)運(yùn)動(dòng)軌跡圖（Pg波到達(dá)時(shí)刻為時(shí)間零點(diǎn)，軌跡突變起點(diǎn)對(duì)應(yīng)時(shí)刻為Sg到時(shí)）Fig.5 Trajectory plots of particle motion for HWG station（left）and SSSE station（right）（time zero corresponds to the arrival time of Pg wave，and the starting point of the abrupt change in trajectory corresponds to the arrival time of Sg wave）.

由于地殼巖石各層介質(zhì)密度不同，波速亦會(huì)不同。波在不同介質(zhì)的分界面上按折射定理產(chǎn)生折射，使地震波射線向上彎曲，呈曲線狀[7]。因此，視出射角法估計(jì)出的震源深度多數(shù)會(huì)偏深，特別是隨著震中距增大，估計(jì)的誤差也會(huì)越大。

2.2 深度震相

對(duì)于地殼結(jié)構(gòu)相對(duì)簡(jiǎn)單的區(qū)域，近震深度震相sPg、sPmP 和sPn，以及它們的參考震相Pg，PmP 和Pn 在近震記錄上通?？梢郧宄赜^測(cè)到。在近震地震圖上識(shí)別出這些深度震相及其參考震相后，可以利用測(cè)量它們之間的到時(shí)差或通過(guò)波形對(duì)比的方法相對(duì)精確地確定震源深度[8]。1956年蘇聯(lián)學(xué)者A.A.維琴斯卡雅，提出在近震中可觀測(cè)sPn 震相。之后，H.B.康道爾斯卡婭給出雙層地殼中近震的sPn-P 的走時(shí)方程[9]。Saikia 在實(shí)際數(shù)據(jù)和理論波形中都發(fā)現(xiàn)在確定的震源深度下sPn相對(duì)Pn的到時(shí)差很穩(wěn)定的特征，因此可以用來(lái)約束震源深度[10]。國(guó)內(nèi)自上世紀(jì)70 年代初開始張誠(chéng)、房明山、任克新、王登偉等人先后研究了sPn震相，認(rèn)為采用sPn 和Pn 波的走時(shí)差來(lái)測(cè)定地震深度的方法簡(jiǎn)便準(zhǔn)確[11－16]。

sPn 波是指發(fā)生在地殼內(nèi)的地震產(chǎn)生的S 波入射到地表經(jīng)地面反射時(shí)，其中SV 成分會(huì)轉(zhuǎn)換為P波，之后轉(zhuǎn)換P波入射到莫霍面，當(dāng)入射角為臨界角時(shí)，在地幔頂部沿著Moho 界面向前滑行，形成Pn 波。拾取sPn 的方法基本上是在同一時(shí)間窗內(nèi)，將不同震中距的Pn 震相對(duì)齊，則sPn 震相也大致對(duì)齊了，再根據(jù)其動(dòng)力學(xué)特征，即周期和振幅都比Pn大的特征，可拾取sPn。洪星，蔡杏輝，吳國(guó)瑞等人給出了華南地區(qū)的sPn 震相動(dòng)力學(xué)的初步特征，深度與sPn-Pn 的關(guān)系滿足（4）式[15-18]。依據(jù)這些特征我們識(shí)別出12 個(gè)臺(tái)的sPn（見表2 及圖6a），sPn 和Pn 的到時(shí)差最大1.55 s，最小1.17 s，平均1.36 s，對(duì)應(yīng)的深度為3.75 km，標(biāo)準(zhǔn)偏差是0.39 km（見表2）。

表2 本次地震識(shí)別出的sPn震相Table 2 sPn seismic phases identified in this earthquake

圖6 原始記錄圖中識(shí)別出的sPn（a）及臺(tái)站（震中距200～300 km）垂向Pn的疊加圖（b）Fig.6 sPn identified in the original record（a）and the overlay map of vertical Pn of stations（epicentral distance between 200 and 300 km km）（b）

由于sPn 震相易受噪音和尾波的影響，單臺(tái)挑選的可靠性不高。為了提高sPn-Pn測(cè)定的準(zhǔn)確性，采用滑動(dòng)時(shí)窗相關(guān)法讀取sPn 震相?；瑒?dòng)時(shí)窗相關(guān)法常用來(lái)估計(jì)信號(hào)之間相關(guān)性大小隨時(shí)間變化的程度。移動(dòng)時(shí)窗相關(guān)法的具體方法如下：對(duì)于一個(gè)地震事件的兩條波形記錄，從某一起始時(shí)間，計(jì)算一定時(shí)窗范圍內(nèi)其相關(guān)性的大小。然后時(shí)窗向前推移（重疊），計(jì)算新時(shí)窗內(nèi)對(duì)應(yīng)波形的相關(guān)值。對(duì)于一個(gè)地震的N條波形記錄，每?jī)蓚€(gè)波形記錄進(jìn)行互相關(guān)計(jì)算，然后再進(jìn)行總體疊加?；ハ嚓P(guān)系數(shù)的大小反映了波形的相似度。把Pn 波初至對(duì)齊后，可利用互相關(guān)峰值的時(shí)移值來(lái)確定sPn 和Pn的到時(shí)差[19]。

本次地震短周期面波比較發(fā)育，初步判斷不可能太深，結(jié)合實(shí)際的波形看：震中距在200～300 km的臺(tái)站，信噪比較高，且這個(gè)范圍內(nèi)Pb和Pg等續(xù)至震相與初至震相Pn 的到時(shí)差已超過(guò)3 s，結(jié)合視出射角法估計(jì)的結(jié)果初步推測(cè)，此范圍的資料可用于搜索本次地震的sPn 震相。對(duì)符合條件的臺(tái)站資料進(jìn)行0.1 至1.5 Hz 帶通濾波，以其中Pn 到時(shí)最小的臺(tái)站記錄為基準(zhǔn)，此到時(shí)為零時(shí)，取0.2 秒為窗長(zhǎng)的數(shù)據(jù)為參考波形，其余的每個(gè)臺(tái)站記錄以相同窗長(zhǎng)按1 個(gè)采樣點(diǎn)滑移與參考波形做互相關(guān)，使Pn 對(duì)齊之后疊加即獲得高信噪比的波形記錄；同時(shí)滑移相關(guān)系數(shù)也對(duì)齊疊加，疊加互相關(guān)系數(shù)在零時(shí)刻附近兩個(gè)峰值的時(shí)差就是所要提取的sPn和Pn 的到時(shí)差。最終搜索到sPn 和Pn 的到時(shí)差為1.53 s，對(duì)應(yīng)深度是4.2 km，詳細(xì)結(jié)果見圖6（b）。

利用人工分析識(shí)別出的sPn 和Pn 到時(shí)差確定的本次地震震源深度為：3.8±0.39 km，通過(guò)移動(dòng)時(shí)窗相關(guān)方法搜索獲取的sPn與Pn的到時(shí)差，估算出震源深度是4.2 km，兩種方法基本一致，因此可確定這次地震的震源深度為3.8 km，相對(duì)誤差約10%。

2.3 微觀震中位置的精確測(cè)定

在震源深度由深度震相確定的情況，采用固定深度，只反演發(fā)震時(shí)刻和震中經(jīng)緯度，減少反演參數(shù)提高定位的穩(wěn)定性和精度。將震源深度固定為3.8 km，使用華南走時(shí)表，利用單純型、hyposat 和LocSAT 進(jìn)行重新定位，3 種方法的定位結(jié)果間兩兩相距約0.7 km，與國(guó)家臺(tái)網(wǎng)中心CENC正式目錄發(fā)震時(shí)刻分別提前0.50 s、0.43 s 和0.36 s，水平位置沿東南向分別偏移1.1 km、0.8 km 和1.1 km，都在麥家村與隔坑村連線上，更靠近隔坑村一側(cè)。為減小誤差和不確定性，提高結(jié)果的可靠性取3種定位方法的定位結(jié)果的平均值作為最終結(jié)果。即：此次地震的發(fā)震時(shí)刻11∶38∶14.84，經(jīng)度112.958 3°N，緯度23.271 6°E，深度3.8 km，震級(jí)ML4.0，詳見表3 和圖7。從圖7 可以看出本次地震發(fā)生在隔坑礦區(qū)北東邊緣處的F7斷層上。

表3 2023年2月5日ML4.0三水地震重新修定的發(fā)震時(shí)刻和微觀震中Table 3 Revised origin time and micro-epicenter of the Sanshui ML4.0 earthquake on February 5th，2023

圖7 2023年2月5日ML4.0三水地震新修定震中與F7斷層位置Fig.7 Location of the revised epicenter of the Sanshui ML4.0 earthquake and F7 fault

3 震源機(jī)制及發(fā)震構(gòu)造

3.1 震源機(jī)制

從收集到的波形資料中識(shí)別出P 波初動(dòng)79 個(gè)，這些初動(dòng)呈典型的四象限分布，加入16 個(gè)震幅比資料和16個(gè)Sg波初動(dòng)，使用focmec方法搜索到了該地震的震源機(jī)制解，見圖8。該地震的震源機(jī)制解為：節(jié)面Ⅰ走向334.90°，傾角54.04°，滑動(dòng)角-18.84°；節(jié)面Ⅱ走向74.84°，傾角76.23°，滑動(dòng)角-142.53°；P 軸方位301.68°，仰角36.88°；T 軸方位201.45°，仰角13.32°；B軸方位95.06°，仰角50.00°。顯然這是一次典型的走滑型構(gòu)造活動(dòng)，可能的最大水平壓力方向?yàn)?90°，這與康英等人給出的該區(qū)域的主壓力σ1較一致[20]，表明區(qū)域應(yīng)力場(chǎng)在本次地震孕育過(guò)程中起主導(dǎo)作用。

圖8 三水ML4.0級(jí)地震的震源機(jī)制解Fig.8 Focal mechanism solution of the Sanshui ML4.0 earthquake

3.2 地質(zhì)構(gòu)造環(huán)境及發(fā)震構(gòu)造確定

這次地震發(fā)生在珠江三角洲西北部的三水盆地北部，三水盆地是珠江三角洲的主體盆地，走向北北西，呈菱形，面積約3300 平方公里?；子稍沤?、下古生界變質(zhì)巖、上古生界砂頁(yè)巖和碳酸鹽巖、中生界砂頁(yè)巖和零星出露的燕山期花崗巖組成。上覆沉積厚約3600 m 的白堊系地層和厚約3200 余米的老第三紀(jì)地層。據(jù)物探調(diào)查，三水盆地老第三紀(jì)早期曾是殼下活動(dòng)和下地幔隆起最活躍的地方[21]。在晚第三紀(jì)時(shí)的三水盆地仍有張性構(gòu)造活動(dòng)，其活動(dòng)性亦明顯減弱。早第三紀(jì)的噴出巖比較廣泛地見于三水獅山、南海市王借崗等地，其地表露頭或地質(zhì)鉆探常顯示屬于基性的玄武巖被酸性或中性的粗面巖、安山巖所披覆，說(shuō)明受上地幔上沖影響的火山巖活動(dòng)已明顯地隨時(shí)間推移而減弱[22]。

震中區(qū)主要控制斷層處于北西向西江斷裂和白坭—沙灣斷裂之間，被東西的廣三斷裂切割，見圖9（a）。

圖9 三水ML4級(jí)地震震中附近的主要斷裂分布Fig.9 Distribution of main faults near the epicenter of the Sanshui ML4 earthquake

西江斷裂基本沿西江下游的北西向河道發(fā)育，南東段由磨刀門入南海，北西段可能沿綏江斷續(xù)延伸至廣寧、懷集一帶，全長(zhǎng)近200 km。斷裂以九江為界，北西段基本沿西江右岸（西岸）發(fā)育，南東段則沿左岸（東岸）發(fā)育，總體產(chǎn)狀是走向NW30°～50°，傾向NE，傾角大于70°[23]。基于淺層地震探測(cè)，西江斷裂未發(fā)現(xiàn)切穿第四系地層，該區(qū)已有第四系底部沉積物年齡40 000 a左右，可見，該斷裂晚更新世中期以來(lái)（40 000.BP）未發(fā)生明顯活動(dòng)[23]。

白坭—沙灣斷裂部分地控制了三水盆地及珠江三角洲第四紀(jì)沉積的北東邊界。斷裂北起花縣白坭向南東經(jīng)南海縣官窯、松崗至番禺沙灣沿蕉門出伶仃洋斷裂走向NW30°～40°，傾向SW，傾角大于50°，是一條第四紀(jì)以來(lái)活動(dòng)的斷裂[24-25]。

廣州—三水隱伏斷裂東段稱羅浮山斷裂（含瘦狗嶺斷裂）[26]。該斷裂部分地控制了晚更新世以來(lái)珠江三角洲的北界，沿?cái)鄬訋r石硅化、糜棱巖化、擠壓破碎等。在廣州以西至三水附近該斷裂帶雖然被第四紀(jì)沉積隱伏，但在MSS6 波段的衛(wèi)星像片上，該斷裂以明顯的東西向深色線性條帶清晰地顯示出來(lái)，并與三水以西的東西向斷裂相連，說(shuō)明了該斷裂全新世以來(lái)的活動(dòng)[25-26]。

這次地震發(fā)生在三水市以北約20 km 處的南邊鎮(zhèn)麥家村與隔坑村之間，靠近隔坑村一側(cè)，從圖7和圖8 看震中位于隔坑鹽礦礦區(qū)北東邊緣與F7 斷層交匯處。據(jù)廣東省地質(zhì)工程公司的技術(shù)報(bào)告，F(xiàn)7 斷層與白坭—沙灣斷裂相距不足2 km，是一條正斷層，其切割深度超過(guò)2 km，走向?yàn)镹NW345°，傾向SW，傾角約50°～60°[27]。與震源機(jī)制解的節(jié)面Ⅰ吻合較好，據(jù)此推斷本次地震應(yīng)該是F7 斷層錯(cuò)動(dòng)引起的，結(jié)合反演出滑動(dòng)角認(rèn)為其錯(cuò)動(dòng)類型為：左旋走滑為主，兼少量正斷。

4 結(jié)論及討論

（1）1997年9月23日和26日隔坑村附近曾分別發(fā)生了ML3.3 和ML4.4 的地震，李純清等人分析認(rèn)為其震相特征與一般構(gòu)造地震不同[28]，廖華康認(rèn)為這兩次三水地震，不可能是天然地震，極有可能是一場(chǎng)塌陷地震[29]。之后該區(qū)域的地震類型一直存在爭(zhēng)議。直至本次地震發(fā)生后，利用近年來(lái)新建的密集觀測(cè)站網(wǎng)收集到了豐富的波形資料，綜合分析判斷排除此次事件為塌陷、滑坡或爆破的可能。比較兩次ML4 以上地震在XNH、ZHQ 和HUD 三臺(tái)的記錄，發(fā)現(xiàn)相似度較高，相關(guān)系數(shù)分別是0.90、0.89和0.73，它們應(yīng)該是同類型地震。

（2）依據(jù)Sn與Sg走時(shí)曲線交點(diǎn)對(duì)應(yīng)的震中距較大、近臺(tái)視出射角以及較遠(yuǎn)臺(tái)記錄的sPn與Pn的走時(shí)差變化等震相特征，修訂了2023 年2 月5 日廣東三水ML4.0 級(jí)地震的發(fā)震時(shí)刻和震源位置為：11∶38∶14.84，經(jīng)度112.958°N，緯度23.270°E，深度3.8 km，震級(jí)ML4.0。震中位于隔坑鹽礦礦區(qū)北東邊緣與F7 斷層交匯處，由深度震相確定的深度比采鹽井深1 km左右，距最近的隔4號(hào)開采井水平距離約500 m。從空間范圍上可以確定地震發(fā)在開采區(qū)，這是典型的與采礦活動(dòng)密切相關(guān)的礦山地震[27]。

（3）結(jié)合震源機(jī)制和已探明的小區(qū)域地質(zhì)構(gòu)造資料，推斷這次地震的發(fā)震構(gòu)造走向?yàn)?34.90°，傾角為54.04°，滑動(dòng)角為-18.84°，可能是F7 斷層沿?cái)鄬用姘l(fā)生左旋走滑正斷錯(cuò)動(dòng)而引發(fā)的。長(zhǎng)期的抽取鹵水采鹽，采空區(qū)應(yīng)力變化引起礦體的應(yīng)力分布會(huì)發(fā)生調(diào)整，沿F7 斷層的應(yīng)力超過(guò)了其承載能力，就使上盤發(fā)生破裂和滑動(dòng)，從而引發(fā)地震。魏柏林等人研究1997 年9 月2 次三水地震認(rèn)為地震的最大主壓應(yīng)力軸近于直立，與重力方向一致，產(chǎn)生的是塌陷地震[30]。我們研究結(jié)果似乎不全支持其結(jié)論，雖然都是沿F7 斷層滑動(dòng)，但最大主壓應(yīng)力軸并非與重力方向一致，而是與該區(qū)域的主壓力σ1較一致，即是區(qū)域應(yīng)力主控制下的礦體應(yīng)力調(diào)整。重力主導(dǎo)下礦區(qū)應(yīng)力調(diào)整產(chǎn)生的礦震應(yīng)該多數(shù)是孤立的，與周邊的中強(qiáng)震時(shí)空關(guān)系不明顯，但區(qū)域應(yīng)力場(chǎng)控制下的調(diào)整則會(huì)與周邊的中強(qiáng)震存在某種看似巧合的時(shí)空關(guān)系。

（4）深度為3.8 km 是一次罕見的淺源地震。據(jù)廣州市地震活斷層探測(cè)與地震危險(xiǎn)性評(píng)價(jià)時(shí)深反射地震勘探專題報(bào)告給出白坭—沙灣斷裂附近的剖面淺部（4 km 以上）成層性較好，層位較多，地層結(jié)構(gòu)簡(jiǎn)單，總體趨勢(shì)為西高東低的單斜形態(tài)，巖性可能是沉積巖層。震源在沉積層產(chǎn)生的P波與SV 波大于臨界角入射，由P 波與SV 波耦合在自由地表形成的波動(dòng)即為短周期的瑞利面波Rg，從觀測(cè)點(diǎn)的空間分布看Rg 波的發(fā)育程度似乎還與震中距和臺(tái)站附近的沉積層厚度相關(guān)。本文初步分析廣東地區(qū)Rg 波的最大峰值的視速度（周期在1.1～1.5 s 間近似的群速度）約為2.85 km/s，峰值振幅隨震中距的衰減關(guān)系為：log（Vmax）=-1.8004*log（Δ）+5.0812。從三水地區(qū)的震例研究認(rèn)為極淺源地震激發(fā)的Rg 波是造成震級(jí)小、震中烈度高、破壞性大和選擇性破壞明顯為其特征[21]的主要因素，在采礦區(qū)附近或有定沉積層建設(shè)場(chǎng)地的抗震設(shè)防工作中應(yīng)該給予足夠的重視。

附2023 年2 月5 日三水ML4.0 級(jí)地震獲取廣東地區(qū)Rg波初步特征

附圖1 廣東三水地震的Rg波峰值速度對(duì)應(yīng)走時(shí)隨震中距變化Fig.1 Variation of travel time corresponding to peak velocity of Rg wave with epicentral distance of Sanshui earthquake，Guangdong