河北地區(qū)地殼一維速度模型的新近研究

王莉嬋,朱元清,李雪英,王 鐸,蔡玲玲,李冬圣,楊 銳

(1.河北省地震局，石家莊 050021；2.上海市地震局，上海 200062；3.唐山市地震局，河北 唐山 063000)

?

河北地區(qū)地殼一維速度模型的新近研究

王莉嬋1,朱元清2,李雪英1,王 鐸3,蔡玲玲1,李冬圣1,楊 銳1

(1.河北省地震局，石家莊 050021；2.上海市地震局，上海 200062；3.唐山市地震局，河北 唐山 063000)

為建立更適合河北地區(qū)地震速報及編目的地殼一維速度模型，收集河北數(shù)字地震臺網(wǎng)2009年1月至2014年10月定位臺站數(shù)N≥10，震級M≥1.0的1 113條地震以及精選的48條地震；利用滑動窗擬合、折合走時等方法及結(jié)合前人研究成果確定一維速度模型中各參數(shù)的范圍，并利用Hyposat批處理程序，對精選的48條地震進行了58 902次組合計算，對殘差最小的10組數(shù)據(jù)利用與編目結(jié)果震中位置偏差大小為原則進行篩選，由此建立了適應河北地區(qū)地震速報及編目的地殼一維速度新模型；然后利用批處理殘差大小比較、PTD深度定位、爆破檢驗、典型地震檢驗等方法對其驗證。結(jié)果表明，相較于華南模型，新模型更適用于河北地區(qū)。

速度模型；滑動窗擬合；折合走時；Hyposat批處理

0 引言

河北地區(qū)位于華北地震帶北部，地質(zhì)構(gòu)造復雜，東側(cè)為太行山脈，北部為陰山—燕山隆起，中部是冀中平原新生代沉降帶，陰山—燕山隆起與冀中平原新生代沉降帶之間是中、新生帶燕山沉降帶。受地質(zhì)構(gòu)造影響，河北地區(qū)斷裂帶十分發(fā)育，北部有張家口—渤海地震活動帶，西部為山西地震帶[1]，中南部地區(qū)地震帶主要分為NNE-NE和NW向2組，其中又以NNE-NE向為主[2]，特定的地質(zhì)構(gòu)造決定了河北省為地震多發(fā)地區(qū)。根據(jù)近年河北省地震臺網(wǎng)的數(shù)字觀測資料可知，河北省大震多發(fā)，中小地震更是相當頻繁，密集度高。因此，加強河北地區(qū)地震活動觀測研究具有十分重要意義。

地震活動觀測依賴于測震臺網(wǎng)精確的基礎數(shù)據(jù)，而精準的地震定位不僅與定位臺站的數(shù)量、分布以及定位軟件、方法等有關，而且也受到地殼速度模型的制約。目前，河北數(shù)字地震臺網(wǎng)處理地震時在用速度模型有2個——華南速度模型及華北速度模型，前者用于日常編目工作，后者用于地震速報工作。其中華南模型是范玉蘭等[3]利用華南地區(qū)大量的天然地震和人工爆破資料研究得到的雙層平均地殼模型，與河北地質(zhì)情況不符；后者是中國地震局地球物理研究所利用唐山和邢臺老震區(qū)的地震資料計算所得，不符合河北整體區(qū)域地質(zhì)情況。2種模型的不統(tǒng)一性給河北數(shù)字地震臺網(wǎng)進行地震精準定位帶來了困難。雖然前人針對河北地區(qū)地下介質(zhì)速度結(jié)構(gòu)也開展過很多工作，但是由于河北地質(zhì)情況復雜，北部山區(qū)與南部平原地區(qū)地殼厚度差異很大。大部分學者是針對不同的區(qū)域開展不同的研究：如于利民等[4]、張學民等[5]均基于深源遠震體波資料并利用Haskell矩陣傳遞方法計算理論地震圖，與實際測資料擬合對比，前者得到河北北部部分臺站下方的速度結(jié)構(gòu)，后者得到陡河、灤縣等臺站下方的速度模型；張成科等[6]利用地震測深剖面資料，進行射線追蹤走時擬合計算，得到了文安、蔚縣等地的一維速度模型；李強等[7]、于湘?zhèn)サ萚8]、王志鑠等[9]通過層析成像反演方法，利用P波到時數(shù)據(jù)，得到了華北區(qū)域地殼三維速度結(jié)構(gòu)；王帥軍等[10]利用張家口—渤海構(gòu)造帶的測震剖面資料，得到鹽山—大興—延慶、北京—懷來—豐鎮(zhèn)等不同測震剖面的地下速度結(jié)構(gòu)；田曉峰等人[11]利用廣義反演理論，對滄縣及唐山地區(qū)的人工地震測深勘探資料進行反演，得到該區(qū)域的速度結(jié)構(gòu)和莫霍面深度；田寶峰等人[12]利用接收函數(shù)方位變化研究方法，進行地殼各向異性反演計算出4個流動臺站的地下結(jié)構(gòu)；趙博等人[1]利用雙差定位方法分析華北地區(qū)的地震活動過程中，針對華北東、西部分別給出了2個不同的速度模型。通過分析可以看出，以上速度模型結(jié)果或是區(qū)域性的、或是三維速度結(jié)構(gòu)，難以直接應用于河北臺網(wǎng)的實際工作中。所以為提高河北臺網(wǎng)地震定位精度，建立一個符合河北實際地質(zhì)情況，且可直接應用于河北地震臺網(wǎng)編目和速報工作的速度模型十分重要，具有很強的現(xiàn)實意義。

1 研究思路及方法

本文通過對大量地震震相資料進行處理分析，建立適用于河北地區(qū)地震速報及編目的新地殼速度模型。首先，收集河北地震臺網(wǎng)“十五”數(shù)字化工程項目以來至本文研究期間的地震目錄及震相資料，利用不同的處理方法判定速度模型中各參數(shù)的變化范圍；其次，基于上述研究成果設定模型參數(shù)范圍及各個參數(shù)的變化步長，利用Hyposat批處理程序，對選取的地震數(shù)據(jù)進行批處理，根據(jù)結(jié)果殘差大小及震中偏差大小的合理性分析，建立新的速度模型；最后，利用批量定位殘差比較、PTD深度定位、人工爆破數(shù)據(jù)、典型地震事件對新建立模型進行驗證，分析新模型的合理性。

1)線性擬合方法

由地震波走時方程可知，Pg、Pb、Pn3種震相走時分別為雙曲線和直線，其漸近線斜率或直線斜率即為各震相的擬合速度。所以通過提取各震相的震中距及走時數(shù)據(jù)，利用近似走時方程計算得到各種震相的擬合速度：

T=Δ/V+b。

(1)

式中：T為震相走時，V為擬合速度，Δ為震中距，b為常數(shù)。

2)折合走時方法

Pb、Pn的理論走時方程為

(2)

(3)

結(jié)合Pn、Pb走時路徑(圖1)，對公式(1)、(2)、(3)進行探討可知：震相走時可分為2部分，一部分與震中距及波速相關(與地震波橫向傳輸距離相關——Δ/V)；一部分與震中距無關(與介質(zhì)厚度、波速、震源深度有關——b)。所以震相走時=橫向走時+豎向走時(折合走時)，對公式(1)進行變換，得到公式

TZ=T-Δ/V 。

(4)

式中：TZ為折合走時；T為理論/觀測震相走時；Δ為震中距；V為波速。

圖1 Pn、Pb走時路徑

因為折合走時只與介質(zhì)厚度(康拉德界面或莫霍面深度)、震源深度和波速相關，所以同一震源深度的情況下，理論折合走時應為常數(shù)值b。本文中利用設置不同的波速和康拉德界面及莫霍面深度，通過觀察理論折合走時與觀測折合走時的分布對速度模型進行調(diào)整。

3)Hyposat批處理方法

Hyposat定位方法是基于傳統(tǒng)蓋格基本原理進行地震定位的方法[13]，Hyposat批處理軟件通過調(diào)用MSDP中的Hyposat定位程序，針對一維速度模型中每個參數(shù)，設置各自的參數(shù)范圍以及不同的變化步長，不斷更改速度模型，對大量地震數(shù)據(jù)同時進行Hyposat定位處理，并對不同模型結(jié)果中的殘差進行比較，選擇較為合理的速度模型。

4)PTD方法

PTD方法是朱元清等人[14]提出來的，是將不同臺站的初至震相到時(包括Pg、Pn2種震相)做變換后，根據(jù)Pg與Pn的到時差組對計算地震的震源深度。該方法對速度模型要求較高，當速度模型越接近實際情況，PTD程序計算結(jié)果中的有效震相數(shù)量也就越多，且圖形結(jié)果越符合高斯分布。

2 模型中各參數(shù)范圍確定

河北臺網(wǎng)目前使用的速度模型主要分為上地殼、下地殼以及上地幔頂層等3層結(jié)構(gòu)，包括以下8個參數(shù)：上地殼Pg波、Sg波速度及康拉德界面(conrad)深度Hcon，下地殼Pb波、Sb波速度及莫霍面(moho)深度Hmoho，莫霍面頂層Pn波、Sn波速度。由于S波是后續(xù)震相，在河北臺網(wǎng)日常編目中Sn震相數(shù)量很少、Sb震相無識別，無法直接利用Sb、Sn震相進行速度求取。本文首先利用P波震相求取一維P波速度模型，在后續(xù)批處理過程中，利用波速比求得S波速度，并確定最終速度模型。

2.1 數(shù)據(jù)資料選取

統(tǒng)計了2009年1月1日至2014年10月31日期間，河北省數(shù)字地震臺網(wǎng)記錄到的臺站數(shù)N≥10(包括鄰省臺站)，且M≥1.0的河北省內(nèi)地震(不包括北京、天津)，共計1 113條。按震級分類， M1.0～2.0的地震有1 001條，M2.1～3.0的地震有96條，M3.1～4.0的地震有13條，M4.1～5.0的地震有3條。按震相統(tǒng)計，1 113條地震記錄到32 911個Pg震相，6 198個Pn震相，44 279個Sg震相，632個Sn震相。

2.2 VPg、VPn范圍的確定

將收集到的1 113條地震資料，提取每個地震的Pg、Pn震相資料的震中距及走時，求取不同震相在不同震中距范圍內(nèi)的擬合速度和穩(wěn)定性。針對Pg、Pn出現(xiàn)震中距不同，分別自50km、150km開始，根據(jù)以下2種情況進行滑動擬合，分析不同震中距范圍內(nèi)2種震相的速度變化特征。

1)線性擬合

表1VPg、VPn隨震中距增大時的變化為VPg、VPn以50km為步長，不斷擴大震中距的擬合結(jié)果：①當震中距在50～150km范圍內(nèi)，Pg速度最小為6.04km/s，隨著震中距范圍不斷增大，速度也不斷增大，最大達到6.21km/s，平均速度為6.17km/s。②當震中距在150～300km范圍內(nèi)，Pn速度最小為7.96km/s，隨著震中距增大，Pn速度也在不斷增大，最大達到8.05km/s，平均速度為8.02km/s。

表1 VPg、VPn隨震中距增大時的變化

2)滑動窗擬合

表2為Pg、Pn震相分別以150km、200km為窗口，50km為步長的震中距變化范圍內(nèi)的變化特征。由表2看出，Pg、Pn震相分別在震中距50～400km、150～500km范圍內(nèi)變化較為穩(wěn)定，即近距離范圍內(nèi)較穩(wěn)定，之后由于震相數(shù)量減少，波動較大。針對不同滑動窗的擬合結(jié)果，舍棄數(shù)據(jù)量使用率低于10%的數(shù)據(jù)結(jié)果后進行統(tǒng)計：第1種滑動窗結(jié)果中VPg的變化范圍為：6.10～6.38km/s，平均值為6.26km/s，VPn的變化范圍為：7.96～8.16km/s，平均值為8.05km/s；第2種滑動窗結(jié)果中VPg的變化范圍為：6.16～6.35km/s，平均值為6.26km/s，VPn的變化范圍為：7.96～8.21km/s，平均值為8.08km/s。結(jié)合線性擬合及滑動窗擬合結(jié)果，確定VPg、VPn的變化范圍分別為：6.0～6.4km/s，7.9～8.3km/s。

2.3 VPb范圍的確定

選取河北省內(nèi)M≥2.0，定位過程中本省所屬臺站使用率在40%以上，且震中分布能覆蓋全省大部分地區(qū)的地震，共計48條(圖2a)。其中河北北部的承德地區(qū)選用的地震為震級M≥1.7，本省臺站使用率在17%以上(圖2b)。根據(jù)張惠[15]的方法，明確本區(qū)域Pb震相特征(圖3a)，對選定的48條地震重新進行震相標定，共標定出96個Pb震相。對此96個震相做線性擬合，擬合結(jié)果VPb=6.56km/s(圖3b)。

表2 VPg、VPn在不同震中距內(nèi)的變化

a 研究選用的地震震中分布 b 研究選用的地震射線分布圖2 河北省選用的48條地震震中分布及選用的地震射線分布圖

a 密云臺Pb震相 b 96個Pb震相的線性擬合圖3 本區(qū)域典型Pb震相圖與Pb震相擬合結(jié)果

根據(jù)折合走時理論，即當速度參數(shù)合適時，同一震源深度下，實際震相應均勻分布在理論折合走時兩側(cè)。統(tǒng)計歷年地震定位結(jié)果可知，河北地區(qū)的地震震源深度大部分位于5～20km之間，所以本文將折合走時中的震源深度設置為5km、10km、15km、20km4個深度值?；谌A南模型、華北模型與擬合結(jié)果中VPb值，利用識別的96個Pb震相做折合走時，不斷調(diào)整VPb大小，當VPb=6.68km/s時，實際震相均勻?qū)ΨQ的分布在理論折合走時兩側(cè)，說明VPb應在6.68km/s左右(圖4c)。結(jié)合華南模型、華北模型以及Pb震相擬合速度結(jié)果和折合走時調(diào)整結(jié)果，認為VPb的變化范圍為：6.5～6.9km/s。

圖4 VPb調(diào)整前后折合走時

2.4 Hcon、Hmoho范圍的確定

康氏界面深度與莫霍面深度可以利用爆破資料的折合走時進行調(diào)整。

是以華南模型為基礎，結(jié)合上文中VPb的調(diào)整結(jié)果，對Hcon、Hmoho進行調(diào)整前后的折合走時。圖4a、4b分別是以爆破震相資料基礎，對華南模型中Hcon、Hmoho進行調(diào)整前后所做的折合走時；圖4c、4d分別是以48條地震震相資料為基礎，對華南模型中Hcon、Hmoho進行調(diào)整前后所做的折合走時。通過對比可以看出，4幅圖中實際震相的折合走時與理論的基本平行，但華南模型的結(jié)果沒有均勻的分布在理論折合走時的兩側(cè)，對Hcon、Hmoho進行調(diào)整，Hcon由21km增加為22km，Hmoho由33km增加為35km，調(diào)整后理論與實際震相的折合走時分布較為均勻，所以Hcon、Hmoho的范圍分別應該在22、35km左右。

根據(jù)查閱文獻可知，由于康拉德界面存在不連續(xù)性，前人針對其研究較少，所以利用爆破折合走時及華南模型、華北模型等結(jié)果，將Hcon的參數(shù)變化范圍定為18～23km。

前人針對莫霍面深度研究較多，其中羅艷等[16]利用接收函數(shù)得到了河北地區(qū)67個臺站下方的莫霍面深度，其分布位于30～42km之間，與大多數(shù)研究成果相符，且上文爆破折合走時結(jié)果正好位于該范圍內(nèi)，故確定河北地區(qū)的Hmoho參數(shù)變化范圍為30～42km。

3 河北區(qū)域一維速度模型建立

根據(jù)一維速度模型中各參數(shù)變化范圍，對選定的地震數(shù)據(jù)進行Hyposat批處理，根據(jù)Hyposat定位結(jié)果殘差大小，以及與編目結(jié)果震中偏差的大小等來確立河北區(qū)域最終的速度模型。本次批處理過程分為2步：首先利用固定的波速比進行定位，尋找殘差較小且較為合理的P波速度模型為最優(yōu)；其次，在此模型基礎上，通過調(diào)整波速比的方法，確定最終速度模型。

3.1 P波速度模型確立

本次批處理過程使用的是上文的48條地震，震相數(shù)量分別為：3 206個Pg震相，1 258個Pn震相，3 380個Sg震相，155個Sn震相，計算過程中將波速比固定為理論波速比1.73。

Hyposat批處理過程分為2步：

第1步，將上文得到的VPg、VPb、VPn、Hcon、Hmoho的變化范圍設置為批處理各參數(shù)的變化范圍，并對速度參數(shù)和厚度參數(shù)分別取步長0.1km/s、1km進行網(wǎng)格搜索運算，尋找殘差最小的速度模型。以下為第一步的取值范圍及步長：

VPg取6.00～6.40km/s，步長為0.1km/s；

VPb取6.50～6.90km/s，步長為0.1km/s；

VPn取7.90～8.30km/s，步長為0.1km/s；

Hcon取18～23km，步長為1km；

Hmoho取30～42km，步長為1km。

第1步殘差最小的結(jié)果為VPg=6.20km/s、VPb=6.50km/s、VPn=7.90km/s、Hcon=23km、Hmoho=35km，殘差為0.417。

第2步：以第1步結(jié)果為基礎，將速度參數(shù)步長調(diào)整為0.01km/s，厚度參數(shù)步長0.5km，各參數(shù)取不同的變化范圍，進行進一步的批處理：

VPg取6.10～6.25km/s，步長為0.01km/s；

VPb取6.45～6.60km/s，步長為0.01km/s；

VPn取7.85～8.00km/s，步長為0.01km/s；

Hcon取21～24km，步長為0.5km；

Hmoho取34～36km，步長為0.5km。

表3為殘差較小的前10組模型和華南模型、華北模型的Hyposat定位結(jié)果以及每組結(jié)果與編目結(jié)果的震中偏差分布。經(jīng)過比較，可以看出第4組模型雖然殘差略大，但是其震中偏差小于2km的地震條數(shù)最多，為84.8%，平均震中偏差為最小值1.9km，故最終認為第4組為最優(yōu)P波速度模型：VPg=6.11km/s，VPb=6.60km/s，VPn=7.96km/s，Hcon=22km，Hmoho=35km。

表3 Hyposat批處理殘差最小的前10組模型

3.2 新模型確立

一般情況下而言，隨著深度的增加，壓力和溫度增大，導致縱波速度增加會快于橫波速度[17]，因而縱、橫波波速比會隨深度增加而增大。根據(jù)這一原則，在P波最優(yōu)模型基礎上，以0.01步長調(diào)整波速比進行批處理計算，利用Sg、Sn震相確立河北新模型，得到殘差最小的3層波速比分別為1.72、1.73、1.76。經(jīng)計算，本次結(jié)果與編目結(jié)果的震中偏差為1.9km，與P波最優(yōu)模型結(jié)果的震中偏差基本一致，所以最終認為表4為河北新模型參數(shù)。

表4 河北2層地殼新模型參數(shù)

4 新模型檢驗

4.1 批量定位殘差檢驗

批量定位檢驗是利用河北新模型對河北區(qū)域2009年1月至2014年10月M≥2.0的地震進行Hyposat批處理定位，并與華南模型、華北模型的Hyposat定位結(jié)果作對比，以此驗證新模型的合理性。

2009年1月至2014年10月河北地區(qū)發(fā)生M≥2.0事件共132次，去除邊界地震以及單臺記錄地震和爆破、可疑地震事件，本次Hyposat批處理共使用地震條數(shù)為125條?；?種模型對125條地震分別進行Hyposat批處理定位，并對3個結(jié)果進行定位殘差對比。新模型的殘差范圍在0.195～0.846之間，平均殘差為0.420；華南模型的殘差范圍在0.22～1.462，平均殘差為0.589；華北模型的誤差范圍在0.240～1.922之間，平均殘差為0.603。故認為新模型優(yōu)于華南模型與華北模型。

4.2 PTD深度定位檢驗

為檢驗新模型在河北整個地區(qū)的適用性，在河北西部、東部、南部、中部4個地區(qū)分別選取了發(fā)生在張家口、唐山、贊皇、廊坊等地的4個地震進行PTD深度定位檢驗(表5，圖6)。由表5看出，新模型的深度略微加大，但較為合理，且新模型的有效震相使用率較華南模型分別提高了36.6%、20.85%、29.57%、35.07%，另外,圖6中符合高斯分布的圖形結(jié)果也證明了新模型定位結(jié)果的可靠，由此證明了新模型在河北地區(qū)的合理性。

表5 4個地震的PTD處理結(jié)果

圖6 PTD處理的圖形結(jié)果

4.3 人工爆破檢驗

中國地震局地球物理研究所于2007年12月12日03點00分在河北懷來地區(qū)進行科學探測爆破，河北省測震臺網(wǎng)完整地記錄了此次爆破波形。運用MSDP軟件中Hyposat定位方法，針對河北新模型與原華南模型分別對此次爆破進行定位，通過定位殘差比較以及2種結(jié)果與實測結(jié)果對比，分析新模型的有效性。2種模型對比結(jié)果(表6)顯示，新模型與華南模型的定位結(jié)果與實測結(jié)果震中差分別為1.02km與3.13km，定位精度提高了2km，并且新模型的定位殘差略小于華南模型的殘差。由此證明新模型優(yōu)于華南模型。

表6 懷來爆破實測結(jié)果與華南模型及河北新模型定位結(jié)果比較

4.4 典型地震事件

利用典型地震事件進行模型驗證過程中，采用了2種對比方法：不同模型(河北新模型、華南模型)、同一定位方法對同一地震進行定位；同一模型(河北新模型)、不同定位方法對同一地震進行定位。本文中選擇的典型地震為2009年以來，河北地區(qū)發(fā)生的最大的2次地震以及天津發(fā)生在下地殼的一次地震：2012年5月28日河北唐山M4.8地震、2014年9月6日河北涿鹿M4.3地震以及2016年5月1日天津北辰區(qū)M2.8地震。

1)不同模型、同一定位方法結(jié)果比較

表7為2012年唐山地震與2014年涿鹿地震基于不同模型(河北新模型、華南模型)，運用Hyposat定位方法分別進行定位的結(jié)果。從中可以看出，河北新模型較華南模型的定位結(jié)果中，發(fā)震時刻與震中位置相差不大，但河北新模型的殘差更小。其中，唐山地震與涿鹿地震的新模型定位結(jié)果的震源深度分別為14.3km、15.9km，而李冬圣等[18]針對這2個地震的矩張量反演結(jié)果中震源深度分別為11km、14km，二者相近。

表7 唐山、涿鹿地震在河北新模型與華南模型結(jié)果中的比較

2)同一模型、不同定位方法結(jié)果比較

本次對比選擇的是2013年唐山M4.8地震及2016年天津北辰M2.8地震，根據(jù)上文定位結(jié)果以及地震目錄可知，二者震源分別位于上、下地殼中。根據(jù)波形傳播規(guī)律可知，發(fā)生在下地殼的地震缺少Pb震相，據(jù)此檢查2次地震波形，發(fā)現(xiàn)唐山地震具有明顯的Pb震相，而天津北辰地震則無，所以也證明2次地震分別發(fā)生在上、下地殼中?；诤颖毙履Ｐ?，應用Hyposat定位方法及PTD定位方法對2次地震分別進行定位。其中，前者2種方法的定位結(jié)果中，震源深度分別為14.3km、15.0km；后者的2種震源深度結(jié)果分別為26.5km、26.8km。2種定位方法震源深度結(jié)果相近且符合地震波形特征。由此可見，河北新模型定位結(jié)果穩(wěn)定、準確。

綜合上述結(jié)果，認為河北新模型相較于原華南模型在地震定位過程中殘差更小，精度更高，與河北實際地質(zhì)情況更為相符。

5 結(jié)論與討論

本文在利用數(shù)字化震相資料建立河北區(qū)域一維速度模型過程中，主要得到以下幾點認識：

1)通過滑動窗擬合、折合走時、Hyposat批處理等程序和方法建立河北地區(qū)新的一維速度模型，并利用新舊模型批處理殘差比較、PTD震源深度定位、人工爆破、典型地震等數(shù)據(jù)和方法對新模型進行驗證。結(jié)果表明，相較于河北臺網(wǎng)現(xiàn)用的華南速度模型，新的速度模型定位精度更高，殘差更小。

2)在標定Pb震相過程中，發(fā)現(xiàn)北部山區(qū)地震的Pb震相較為清晰，東南部平原地區(qū)Pb震相較難識別，且整體Pb震相并不發(fā)育，由此說明河北地區(qū)康拉德界面并不連續(xù)，尤其是東南部地區(qū)，地殼結(jié)構(gòu)更為復雜。

3)在建立河北地區(qū)新模型的過程中，發(fā)現(xiàn)河北南部平原與北部山區(qū)地質(zhì)構(gòu)造有所差異，其速度模型也應不同；且在課題研究過程中，基于新模型對唐山地區(qū)部分地震重新做了PTD定位，發(fā)現(xiàn)新模型定位結(jié)果的震源深度較原華南模型結(jié)果偏深，后續(xù)可針對以上2點做進一步研究。

致謝 感謝河北省地震局、上海市地震局以及一維速度模型相關工作人員提供的支持和幫助；感謝審稿老師提出的寶貴意見及建議。

[1] 趙博，高原，石玉濤.用雙差定位結(jié)果分析華北地區(qū)的地震活動[J].地震，2013，33(1):12-21.

[2] 胡幸平，崔效鋒.華北地區(qū)中部地震精定位與構(gòu)造應力場研究[J].震災防御技術(shù)，2013，8(4):351-360.

[3] 范玉蘭，林紀曾，胡瑞賀，等.華南地區(qū)近震走時表的研制[J].華南地震，1990，10(2):1-16.

[4] 于利民，刁桂苓，李欽祖，等.由深源遠震體波記錄反演華北北部地殼上地幔速度結(jié)構(gòu)[J].華北地震科學，1995，13(3):11-20.

[5] 張學民，束沛鎰，刁桂苓，等.利用數(shù)字地震記錄研究唐山震區(qū)臺下的P、S波速度結(jié)構(gòu)[J].華北地震科學， 2001，19(1):10-17.

[6] 張成科，張先康，蓋玉杰，等.文安-蔚縣-察右中旗剖面地殼上地幔速度結(jié)構(gòu)與構(gòu)造研究[J].華北地震科學， 1997，15(3):18-28.

[7] 李強，王椿鏞，劉瑞豐，等.應用層析成像技術(shù)研究華北地殼速度結(jié)構(gòu)[J].地震地磁觀測與研究，1999，20(5) :88-97.

[8] 于湘?zhèn)?，陳運泰，王培德.京津唐地區(qū)中上地殼三維P波速度結(jié)構(gòu)[J].地震學報，2003，25(1):1-14.

[9] 王志鑠，王椿鏞，曾融生，等.華北及鄰區(qū)地殼上地幔三維速度結(jié)構(gòu)的地震走時層析成像[J].CT理論與研究應用， 2008，17(2):15-27.

[10] 王帥軍，張先康，劉寶峰，等.張渤地震構(gòu)造帶中西段及鄰區(qū)深部構(gòu)造探測[J].大地測量與地球動力學， 2005，25(3):110-115.

[11] 田曉峰，張先康，嘉世旭，等.華北北部不同構(gòu)造單元地殼速度模型——均勻分層介質(zhì)反射波走時反演的結(jié)果[J]，華北地震科學，2006，24(3):1-5.

[12] 田寶峰，李娟，姚振興.利用接收函數(shù)方位變化研究太行山區(qū)地殼各向異性[J].地震學報，2008，30(4):355-366.

[13] 陳貴美，楊選，劉錦.廣東數(shù)字地震臺網(wǎng)“十五”系統(tǒng)的幾種地震定位方法的定位效果分析[J].華南地震，2009，29(1):69-78.

[14] 朱元清，夏從俊，李平.確定震源深度的PTD方法及其應用[J].地震地磁觀測與研究，1997，18(3):21-29.

[15] 張蕙.淺析太原盆地地震波形特征[J].山西地震，2012 (2):1-4.

[16] 羅艷，崇加軍，倪四道，等.首都圈地區(qū)莫霍面起伏及沉積層厚度[J].地球物理學報，2008，51(4):1135-1145.

[17] 盧琳，閆桂京，陳建文.地層溫度和壓力對地震波速的影響[J].海洋地質(zhì)動態(tài)，2015，21(9):13-16，21.

[18] 李冬圣，李文軍，王艷榮，等.2014年河北涿鹿M4.3地震矩張量反演[J].華北地震科學，2016，34(1):36-40.

Study on One-dimensional Crustal Velocity Model in Hebei Area

WANG Li-chan1,ZHU Yuan-qing2,LI Xue-ying1,WANG Duo3, CAI Ling-ling1,LI Dong-sheng1,YANG Rui1

(1.Earthquake Administration of Hebei Province,Shijiazhuang 050021,China;2.Earthquake Administration of Shanghai,Shanghai 200062,China;3.Earthquake Administration of Tangshan,Tangshan 063000,China)

In order to establish a more appropriate one-dimensional velocity model for the earthquake quick reporting and compiling catalogue in Hebei Seismometry Network,1113 earthquakes with the positioning station number over 10 and the M magnitude over 1were collected and 48 earthquakes were selected from January 2009 to October 2014.Firstly,the range of parameters in one-dimensional velocity model was got by the methods of the sliding window fitting,the travel time reducing and referring to some previous researches.Secondly,to establish the new model of Hebei by the program of Hyposat Batching with the 48 selected earthquakes for 58902 calculations.The new velocity model was established according to the criterion,that the residual of the calculation was minor and the deviation of epicenter between the calculation and the compiling catalogue was minimum.Thirdly,the new velocity model was examined by the comparison of residual results,the depth positioning by the PTD,the blast events and the typical earthquakes and so on.It was shown that comparing to the South-China mode the new velocity model is more in conformity with the geological feature in Hebei area.

velocity model; sliding window fitting；travel time reducing；Hyposat batching

王莉嬋，朱元清，李雪英,等．河北地區(qū)地殼一維速度模型研究[J]．華北地震科學，2016，34(4):1-10.

2016-08-12

2014年河北省地震科技星火計劃項目(DZ20150424088)

王莉嬋(1988—)，女，河北石家莊人，助理工程師，主要從事地震監(jiān)測工作．E-mail：wlc84872583@163.com

P315.2

A

1003-1375(2016)04-0001-10

10.3969/j.issn.1003-1375.2016.04.001