MSDP定位程序的討論

1、關(guān)于關(guān)于MSDP定位程序的討論定位程序的討論Jopens/msdp中的地震定位程序中的地震定位程序注：注：2為一種改進(jìn)型的遺傳算法為一種改進(jìn)型的遺傳算法序號(hào)序號(hào)定位方法定位方法適用范圍適用范圍理論走時(shí)理論走時(shí)地方震地方震近震近震遠(yuǎn)震遠(yuǎn)震速度結(jié)構(gòu)速度結(jié)構(gòu)模型模型區(qū)域走時(shí)區(qū)域走時(shí)表表全球走時(shí)全球走時(shí)表表1單純型2自適應(yīng)3hypo20004hyposat5locsat6ISCloc7ISCloc28單臺(tái)9Hyposat 3D10Loc3dSB（川滇3D）單純型和自適應(yīng)單純型和自適應(yīng) 全局搜索算法，需要配區(qū)域走時(shí)表和全球走時(shí)表，在定位時(shí)需根據(jù)地震類型（遠(yuǎn)震還是近震）選擇配置好的走時(shí)表（nc格式）HYP

2、2000 適用于地方震和近震的地震定位和定震級(jí)的適用于地方震和近震的地震定位和定震級(jí)的Fortran程序。程序。 HYPOINVERSE-2000用戶指南用戶指南, 2002年年4月版本月版本,F.Klein,USGS Open File Report 02-171，ver1.0 “Hypoinverse是一個(gè)計(jì)算機(jī)程序，它處理一個(gè)地震的地震臺(tái)是一個(gè)計(jì)算機(jī)程序，它處理一個(gè)地震的地震臺(tái)數(shù)據(jù)文件（如數(shù)據(jù)文件（如P波到時(shí)和地震圖上的振幅和持續(xù)時(shí)間），得到波到時(shí)和地震圖上的振幅和持續(xù)時(shí)間），得到地震位置和震級(jí)。它是美國(guó)地質(zhì)調(diào)查局（地震位置和震級(jí)。它是美國(guó)地質(zhì)調(diào)查局（USGS）一系列類似）一系列類似程序

3、之一，這些程序包括程序之一，這些程序包括HYPOLAYR（Eaton，1969），），HYPO71（Lee and Lahr，1971），），HYPOELLIPSE（Lahr，1980）” 在利用水平分層模型時(shí)，采取措施補(bǔ)償?shù)貧M向不均勻性對(duì)在利用水平分層模型時(shí)，采取措施補(bǔ)償?shù)貧M向不均勻性對(duì)走時(shí)的影響（見下面幾張幻燈片）。走時(shí)的影響（見下面幾張幻燈片）。HYP2000：地殼速度模型：地殼速度模型 所有模型都是平地球模型并假定臺(tái)站位于地球表所有模型都是平地球模型并假定臺(tái)站位于地球表面。臺(tái)站高程沒(méi)有被使用，但高程的延時(shí)效應(yīng)大面。臺(tái)站高程沒(méi)有被使用，但高程的延時(shí)效應(yīng)大部分能以臺(tái)站延時(shí)來(lái)補(bǔ)償。部分能

4、以臺(tái)站延時(shí)來(lái)補(bǔ)償。 這一平地球假設(shè)意味著震源深度是相對(duì)于由附近這一平地球假設(shè)意味著震源深度是相對(duì)于由附近的那些地震臺(tái)站確定的平均的局部表面。這在計(jì)的那些地震臺(tái)站確定的平均的局部表面。這在計(jì)算上很簡(jiǎn)單。算上很簡(jiǎn)單。 對(duì)于每個(gè)地殼模型，速度只隨深度改變。只能通對(duì)于每個(gè)地殼模型，速度只隨深度改變。只能通過(guò)根據(jù)位置選擇模型或?qū)τ谕坏卣鸩煌_(tái)站使過(guò)根據(jù)位置選擇模型或?qū)τ谕坏卣鸩煌_(tái)站使用兩個(gè)不同的模型來(lái)補(bǔ)償。用兩個(gè)不同的模型來(lái)補(bǔ)償。圖圖1，定義兩個(gè)地殼模型區(qū)和它們之間過(guò)渡區(qū)的節(jié)點(diǎn)舉例，定義兩個(gè)地殼模型區(qū)和它們之間過(guò)渡區(qū)的節(jié)點(diǎn)舉例對(duì)不同地區(qū)給定不同的地殼模型對(duì)不同地區(qū)給定不同的地殼模型圖圖2.多地殼模

5、型區(qū)域舉例多地殼模型區(qū)域舉例兩類地殼模型兩類地殼模型: 1.均勻?qū)拥貧つＰ途鶆驅(qū)拥貧つＰ途鶆驅(qū)拥貧つＰ臀募母袷绞蔷鶆驅(qū)拥貧つＰ臀募母袷绞切行?格式格式數(shù)據(jù)數(shù)據(jù)行行1 A30 模型名稱模型名稱行行2及以及以后各行后各行 2F5.2 層的速度和到層頂?shù)纳疃葘拥乃俣群偷綄禹數(shù)纳疃?每層使用一行，首先是頂層。到第一每層使用一行，首先是頂層。到第一層頂?shù)纳疃缺仨殲閷禹數(shù)纳疃缺仨殲?.0，最后一層是，最后一層是下伏的半空間。下伏的半空間。Hyp2000運(yùn)行時(shí)，臨時(shí)用給定地殼模型計(jì)算走時(shí)運(yùn)行時(shí)，臨時(shí)用給定地殼模型計(jì)算走時(shí)/Jopens/msdp/work/hyp2000_crhPARKFIELD MOD

6、EL 1.42 0.00 3.24 0.25 4.82 1.50 5.36 2.50 5.60 3.50 5.87 6.00 6.15 9.00 6.60 15.00 8.00 25.00兩類地殼模型兩類地殼模型: 2.使用走時(shí)表的線性梯度模型使用走時(shí)表的線性梯度模型 模型由模型由2-15個(gè)點(diǎn)構(gòu)成，用戶指定在每個(gè)點(diǎn)速度和個(gè)點(diǎn)構(gòu)成，用戶指定在每個(gè)點(diǎn)速度和深度。假定由線性速度梯度連接相鄰點(diǎn)。最后一深度。假定由線性速度梯度連接相鄰點(diǎn)。最后一點(diǎn)給出模型下方均勻半空間的速度，它必須是各點(diǎn)給出模型下方均勻半空間的速度，它必須是各速度值中最大的，以保證沿半空間頂面能生成繞速度值中最大的，以保證沿半空間頂面能

7、生成繞射波。射波。 使用使用hyp2000程序包帶有的程序程序包帶有的程序TTGEN，由給定的，由給定的線性梯度模型計(jì)算出相應(yīng)的走時(shí)表，其后線性梯度模型計(jì)算出相應(yīng)的走時(shí)表，其后hyp2000直接使用走時(shí)表，不必每次用地殼模型計(jì)算走時(shí)。直接使用走時(shí)表，不必每次用地殼模型計(jì)算走時(shí)。HYP2000：臺(tái)站延時(shí)校正：臺(tái)站延時(shí)校正：在在Hypoinverse（HYP2000）中有如下關(guān)系：）中有如下關(guān)系：TOBS = SEC + CCOR - OTRES = TOBS - TCAL DLY= TOBS (TCAL + DLY)其中，其中，SEC = 觀測(cè)到時(shí)，觀測(cè)到時(shí)，CCOR = 時(shí)鐘校正時(shí)鐘校正 (只

8、在傳統(tǒng)格式震相只在傳統(tǒng)格式震相文件中使用文件中使用)，OT = 發(fā)震時(shí)刻發(fā)震時(shí)刻TOBS = 觀測(cè)走時(shí)，觀測(cè)走時(shí)，TCAL = 計(jì)算走時(shí)，是根據(jù)給定地殼模型，取所有臺(tái)站的高程計(jì)算走時(shí)，是根據(jù)給定地殼模型，取所有臺(tái)站的高程均為均為0km計(jì)算出來(lái)的計(jì)算出來(lái)的DLY = 臺(tái)站延時(shí)，包含臺(tái)站高程和臺(tái)站局部地殼結(jié)構(gòu)的貢獻(xiàn)臺(tái)站延時(shí)，包含臺(tái)站高程和臺(tái)站局部地殼結(jié)構(gòu)的貢獻(xiàn)RES =走時(shí)殘差走時(shí)殘差（見（見Hypoinverse使用指南，使用指南，2002年版）年版）HYP2000：臺(tái)站延時(shí)校正：臺(tái)站延時(shí)校正臺(tái)站延時(shí)的輸入方式臺(tái)站延時(shí)的輸入方式 在臺(tái)站文件中寫入。在在臺(tái)站文件中寫入。在HYP2000中的臺(tái)站文件

9、有中的臺(tái)站文件有3種格式：種格式：HYPO71臺(tái)站數(shù)據(jù)格式，臺(tái)站數(shù)據(jù)格式，HYPOINVERSE臺(tái)站數(shù)據(jù)格式臺(tái)站數(shù)據(jù)格式#1和和HYPOINVERSE臺(tái)站數(shù)據(jù)格式臺(tái)站數(shù)據(jù)格式#2，在，在JOPENS/MSDP套用時(shí)采套用時(shí)采用的是用的是HYPOINVERSE臺(tái)站數(shù)據(jù)格式臺(tái)站數(shù)據(jù)格式#2，見，見/jopens/msdp/work/hyp2000-sta 單獨(dú)建立臺(tái)站延時(shí)文件。單獨(dú)建立臺(tái)站延時(shí)文件。 由于在由于在jopens/msdp中是在每次進(jìn)行地震定位時(shí)臨時(shí)生成中是在每次進(jìn)行地震定位時(shí)臨時(shí)生成臺(tái)站文件臺(tái)站文件hyp2000-sta，所以不便于在臺(tái)站文件中引入臺(tái)站，所以不便于在臺(tái)站文件中引入臺(tái)站

10、延時(shí)。延時(shí)。HYP2000：關(guān)于臺(tái)站高程：關(guān)于臺(tái)站高程 在臺(tái)站文件中有臺(tái)站高程的位置，但程序中沒(méi)在臺(tái)站文件中有臺(tái)站高程的位置，但程序中沒(méi)有使用臺(tái)站高程。臺(tái)站高程對(duì)走時(shí)的影響在臺(tái)站延有使用臺(tái)站高程。臺(tái)站高程對(duì)走時(shí)的影響在臺(tái)站延時(shí)中反映。因此，如果將臺(tái)站延時(shí)設(shè)為時(shí)中反映。因此，如果將臺(tái)站延時(shí)設(shè)為0，則表明，則表明不考慮臺(tái)站高程對(duì)走時(shí)的影響。不考慮臺(tái)站高程對(duì)走時(shí)的影響。在在BLOC86、HYPO81BJ程序中程序中對(duì)每個(gè)臺(tái)站計(jì)算走對(duì)每個(gè)臺(tái)站計(jì)算走時(shí)時(shí)用臺(tái)站高程修時(shí)時(shí)用臺(tái)站高程修改第一層的層頂高改第一層的層頂高度度LocSat：用大尺度稀疏臺(tái)網(wǎng)進(jìn)行地震定位：用大尺度稀疏臺(tái)網(wǎng)進(jìn)行地震定位 LocSat的

11、主要目的是利用來(lái)自若干小型地震臺(tái)陣和若干單的主要目的是利用來(lái)自若干小型地震臺(tái)陣和若干單個(gè)地震臺(tái)站的數(shù)據(jù)進(jìn)行地震定位。被定位的事件為個(gè)地震臺(tái)站的數(shù)據(jù)進(jìn)行地震定位。被定位的事件為“區(qū)域區(qū)域地震地震”（到各臺(tái)站或臺(tái)（到各臺(tái)站或臺(tái)zhen的距離在的距離在2000km左右為區(qū)域左右為區(qū)域地震）地震） 臺(tái)陣（以及臺(tái)陣（以及3分向單臺(tái)）給出的震中方位角也用于地震定分向單臺(tái)）給出的震中方位角也用于地震定位，為此引入了震中方位角對(duì)震源位置坐標(biāo)的偏導(dǎo)數(shù)位，為此引入了震中方位角對(duì)震源位置坐標(biāo)的偏導(dǎo)數(shù) STEVEN R. BRATT AND THOMAS C. BACHE, LOCATING EVENTS WITH

12、A SPARSE NETWORK OF REGIONAL ARRAYS，BSSA, 78, 2, pp. 780-798, 1988：文章重點(diǎn)是如何：文章重點(diǎn)是如何在線性化觀測(cè)方程組中引入方位角的空間偏導(dǎo)數(shù)以及誤差在線性化觀測(cè)方程組中引入方位角的空間偏導(dǎo)數(shù)以及誤差估計(jì)算法。估計(jì)算法。 LocSat使用說(shuō)明（使用說(shuō)明（LocSat程序包（程序包（1991年版）年版）確定方位角確定方位角相對(duì)于相對(duì)于直角坐標(biāo)的偏導(dǎo)數(shù)直角坐標(biāo)的偏導(dǎo)數(shù)：a)x,向東，向東，(b) y,向北向北LocSat:走時(shí)計(jì)算走時(shí)計(jì)算 使用多種震相的走時(shí)表，每個(gè)震相對(duì)應(yīng)一個(gè)走時(shí)使用多種震相的走時(shí)表，每個(gè)震相對(duì)應(yīng)一個(gè)走時(shí)數(shù)據(jù)文件，其

13、中按給定深度給出以等間隔遞增震數(shù)據(jù)文件，其中按給定深度給出以等間隔遞增震中距的走時(shí)。中距的走時(shí)。走時(shí)表配置：增加了區(qū)域走時(shí)表的選擇(./lib/locsat/tables/)Hyposat “HYPOSAT是一個(gè)為震源定位而開發(fā)的程序包。是一個(gè)為震源定位而開發(fā)的程序包。它利用了走時(shí)數(shù)據(jù)、反方位角（即臺(tái)站到事件的它利用了走時(shí)數(shù)據(jù)、反方位角（即臺(tái)站到事件的方位角）值和射線參數(shù)值。所考慮的震相是那些方位角）值和射線參數(shù)值。所考慮的震相是那些包含于包含于IASP91震相類型表中的震相，方模型如果震相類型表中的震相，方模型如果使用地，還包括從康拉德和莫霍不連續(xù)界面反射使用地，還包括從康拉德和莫霍不連續(xù)界

14、面反射的震相。程序遵循的震相。程序遵循IASPEI標(biāo)準(zhǔn)震相命名工作組所標(biāo)準(zhǔn)震相命名工作組所推薦的震相名。推薦的震相名?！?P. Bormann主編，中文版新地震觀測(cè)實(shí)踐手主編，中文版新地震觀測(cè)實(shí)踐手冊(cè)，地震出版社，其中含有冊(cè)，地震出版社，其中含有Hyposat的使用說(shuō)明的使用說(shuō)明HYPOSAT: hyposat-parameterJOPENS/MSDP中的中的Hyposat參數(shù)設(shè)置界面截圖參數(shù)設(shè)置界面截圖HYPOSAT: hyposat-parameter以下是在以下是在MSDP中可由用戶設(shè)置的參數(shù)中可由用戶設(shè)置的參數(shù)GLOBAL MODEI（全球模型）：計(jì)算所有與走時(shí)相關(guān)的理論全球模型）：計(jì)

15、算所有與走時(shí)相關(guān)的理論數(shù)據(jù)時(shí)所使用的參考模型的類型。這個(gè)軟件包包括了下面數(shù)據(jù)時(shí)所使用的參考模型的類型。這個(gè)軟件包包括了下面的模型：的模型： akl35 AKl35模型模型(Kennett et al., 1995) iasp91 IASP9l模型模型(Kennett, 1991; Kennett & Engdahl, 1991) jb Jeffreys-Bullen模型模型(Jeffreys & Bullen, 1940 and later) prem PREM模型模型(Dziewonski & Anderson, 1981) sp6 SP6模型模型(Morelli &

16、amp; Dziewonski, 1993) 在輸入數(shù)據(jù)文件在輸入數(shù)據(jù)文件hyopsat-in中給出到時(shí)、方位角、慢中給出到時(shí)、方位角、慢du數(shù)數(shù)據(jù)，每個(gè)數(shù)據(jù)一行，每行中有一組控制字，其中給出采用據(jù)，每個(gè)數(shù)據(jù)一行，每行中有一組控制字，其中給出采用第幾號(hào)全球模型。如不指出（現(xiàn)在如此），則默認(rèn)取第幾號(hào)全球模型。如不指出（現(xiàn)在如此），則默認(rèn)取1號(hào)號(hào)模型。模型。HYPOSAT: hyposat-parameter LOCAL OR REGIONAL MODEL（ 地方或區(qū)域模型地方或區(qū)域模型 （將估算（將估算下下 速度模型文件名速度模型文件名 或區(qū)域或區(qū)域 地方地方 列列震相震相的的走時(shí)（走時(shí)（ 只要

17、它們只要它們?cè)谝欢ǖ木嚯x和震源深度能夠觀測(cè)到）：在一定的距離和震源深度能夠觀測(cè)到）：Pg、Pb、Pn、P、pPg、pPb、pPn、pP、PbP（即在本程序中由（即在本程序中由“康拉德康拉德”上側(cè)的反射震相）、上側(cè)的反射震相）、PmP、PgPg、PbPb、PnPn、PP及轉(zhuǎn)及轉(zhuǎn)換震相換震相sPg、sPb、sPn、sP、SbP，和，和SmP。相同的震相設(shè)。相同的震相設(shè)置分別用于置分別用于S型震相。型震相。 PHASE INDEX FOR LOCAL MODEL（地方模型的震相指數(shù)）：地方模型的震相指數(shù)）：這個(gè)參數(shù)允許用戶指明一組震相，在地方和（或）區(qū)域模這個(gè)參數(shù)允許用戶指明一組震相，在地方和（或）

18、區(qū)域模型中將計(jì)算這組震相的走時(shí)及其偏導(dǎo)數(shù)型中將計(jì)算這組震相的走時(shí)及其偏導(dǎo)數(shù)。 該參數(shù)是一個(gè)該參數(shù)是一個(gè)4位數(shù)字的數(shù)位數(shù)字的數(shù)。只要不被設(shè)置為負(fù)值，直達(dá)只要不被設(shè)置為負(fù)值，直達(dá)震相震相Pg、Pb、Pn、P（或相同的或相同的S震相）總會(huì)被計(jì)算。震相）總會(huì)被計(jì)算。Hyposat的區(qū)域地殼模型：的區(qū)域地殼模型：/jopens/msdp/work/hyposat-crh10. 給出適用的震中距范圍（度）給出適用的震中距范圍（度） 0.000 5.400 3.100 10.000 5.800 3.200 20.000 5.800 3.200CONR 20.000 6.500 3.600 30.000 6.

19、800 3.900MOHO 30.000 8.100 4.500 77.500 8.050 4.400 120.000 8.100 4.500HYPOSAT: hyposat-parameterCRUST 5.1 PATH（CRUST 5.1路徑）：存放路徑）：存放CRUST 5.1數(shù)據(jù)文件數(shù)據(jù)文件（Mooney et a1.,1998）的目錄的路徑。的目錄的路徑。CRUST 5.1（CRUST 5.1參數(shù)參數(shù)）：）：該參數(shù)控制著該參數(shù)控制著CRUST 5.1模型的使用：模型的使用：=0 根本不使用根本不使用CRUST 51； =1 CRUST 5.1用于計(jì)算與臺(tái)站下方地方地殼結(jié)構(gòu)有關(guān)的臺(tái)站用

20、于計(jì)算與臺(tái)站下方地方地殼結(jié)構(gòu)有關(guān)的臺(tái)站校正值；校正值； =2 CRUST 5.1用于定義距離在用于定義距離在6以內(nèi)的震源與臺(tái)站之間的地以內(nèi)的震源與臺(tái)站之間的地方和（或）區(qū)域速度模型：方和（或）區(qū)域速度模型： =3 CRUST 5.1用于定義地方和（或）區(qū)域速度，并對(duì)地球表用于定義地方和（或）區(qū)域速度，并對(duì)地球表面臺(tái)站處和反射點(diǎn)處的地方地殼結(jié)構(gòu)加以校正。如果該參數(shù)被設(shè)面臺(tái)站處和反射點(diǎn)處的地方地殼結(jié)構(gòu)加以校正。如果該參數(shù)被設(shè)置為大于置為大于0的任意值，且模型的任意值，且模型CRUST 5.1可用，則可計(jì)算地球表面可用，則可計(jì)算地球表面反射震相的反射點(diǎn)處地殼結(jié)構(gòu)的時(shí)間校正（例如反射震相的反射點(diǎn)處地殼

21、結(jié)構(gòu)的時(shí)間校正（例如PnPn，sS，pp，）。CURST5.1: Cntype_keycrust_code crust_name v1 v2 v3 v4 v5 v6 v7 v8 - P wavevs1 vs2 vs3 vs4 vs5 vs6 vs7 vs8 - S wavet1 t2 t3 t4 t5 t6 t7 t8 total - ThicknessE1 Jeffreys-Bullen (standard model)0 0 0 0 5.56 0 6.5 7.780 0 0 0 3.37 0 3.74 4.360 0 0 0 20.0 0 15.0 inf. 33E3 IASP91 (st

22、andard model)0 0 0 0 5.8 0 6.5 8.040 0 0 0 3.36 0 3.75 4.470 0 0 0 20.0 0 15.0 inf. 35E5 SP6 (standard model).E6 AK135 (standard model).D7 Platform 4 km seds.3.811.52.546.26.67.38.21.940.01.12.13.63.74.04.70 00.53.514149 inf.41CRT 5.1: CNelevatio-180-175-170-165-160-155-150-145-140-135-130-125-120-1

23、15-110-105-100-95-90-85-80-75-70-65-60-55-50-45-40-35-30-25-20-15-10-50510152025303540455055606570758085909510010511011512012513013514014515015516016517017540-5218-5562-5682-5639-5582-5466-5206-5257-5100-4811-4460-9651533174123261231414245198432139-2390-4310-4862-5095-5201-4995-4601-3741-2288-2462-4

24、215-4773-4174-443-8-974-856-921-2104-747-1334338901210129550974144614943837300924223536341837612341145595373-12-108-1536-587-2792-5600-5717-5249-4964-5355-5137-4754CRUST 5.1: CNtype-180 -175-170-165-160-155-150-145-140-135-130-125-120-115-110-105-100 -95 -90 -85 -80 -75 -70 -65 -60 -55 -50 -45 -40 -

25、35 -30 -25 -20 -15 -10 -5 0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 105 110 115 120 125 130 135 140 145 150 155 160 165 170 17540 A0 A0 A0 A0 A0 A0 A0 A0 A0 A0 A0 S1 N3 N3 H5 D5 D5 DE D5 QA Q9 A8 A4 A2 A2 A2 A0 A1 A0 A0 B0 A0 A0 A2 N4 N4 B3 B8 S4 B8 S3 T7 P1 P1 P1 P1 D2 D4 D4 DB

26、P1 P1 Q0 R0 R0 P4 Q9 D7 GB GC C4 GB A3 J1 K2 A0 A0 A0 A0 A0 A0 A0Hyposat: 走時(shí)計(jì)算走時(shí)計(jì)算 利用利用CRUST5.1 進(jìn)行射線到達(dá)地面點(diǎn)（臺(tái)站或反射進(jìn)行射線到達(dá)地面點(diǎn)（臺(tái)站或反射波的反射點(diǎn)）局部地殼的走時(shí)校正。從而使計(jì)算出波的反射點(diǎn)）局部地殼的走時(shí)校正。從而使計(jì)算出的走時(shí)比使用全球平均模型計(jì)算的走時(shí)更準(zhǔn)確。的走時(shí)比使用全球平均模型計(jì)算的走時(shí)更準(zhǔn)確。HYPOSAT: hyposat-parameter中中已由已由MSDP設(shè)定的部分參數(shù)設(shè)定的部分參數(shù) P-VELOCITY TO CORRECT ELEVATION（校正高程的

27、校正高程的P波速度）：波速度）：如果在如果在STATION CORRECTION FILE中沒(méi)有設(shè)置該參數(shù)，那么該中沒(méi)有設(shè)置該參數(shù)，那么該參數(shù)就是對(duì)臺(tái)站高程進(jìn)行校正的地方參數(shù)就是對(duì)臺(tái)站高程進(jìn)行校正的地方P波速度（波速度（Vpl，默認(rèn)值，默認(rèn)值為為5.8krn/s）。如果）。如果Vpl=99，就不采用臺(tái)站高程校正，且不，就不采用臺(tái)站高程校正，且不使用使用STATION CORRECTION FILE。 S-VELOCITY TO CORRECT ELEVATION（校正高程的（校正高程的S波速度）：波速度）：如果在如果在ELEVATION CORRECTION FILE中沒(méi)有給出該參數(shù)，那么中沒(méi)

28、有給出該參數(shù)，那么該參數(shù)就是對(duì)臺(tái)站高程進(jìn)行校正的地方該參數(shù)就是對(duì)臺(tái)站高程進(jìn)行校正的地方S波速度（波速度（Vsl，默認(rèn)，默認(rèn)值為值為Vpl/3 S-VELOCITY TO CORRECT ELEVATION : 0. （取默認(rèn)值（取默認(rèn)值Vpl/3） P-VELOCITY TO CORRECT ELEVATION : 4.5 （說(shuō)明進(jìn)行了高程校正）（說(shuō)明進(jìn)行了高程校正）HYPOSAT: hyposat-parameter中中已由已由MSDP設(shè)定的部分參數(shù)設(shè)定的部分參數(shù) STARTING SOURCE TIME（初始震源時(shí)間）如果沒(méi)有設(shè)置這初始震源時(shí)間）如果沒(méi)有設(shè)置這個(gè)值，將使用個(gè)值，將使用Wad

29、ati（和達(dá)）（和達(dá)）法根據(jù)直達(dá)法根據(jù)直達(dá)S型與型與P型震相的型震相的觀測(cè)結(jié)果之間的走時(shí)差估算初始震源時(shí)間。對(duì)此，程序分觀測(cè)結(jié)果之間的走時(shí)差估算初始震源時(shí)間。對(duì)此，程序分別為每種震相類型計(jì)算別為每種震相類型計(jì)算vp/vs關(guān)系，并估算了震源時(shí)間。關(guān)系，并估算了震源時(shí)間。初始震源時(shí)間就是所有估算震源時(shí)間的平均值。只有一個(gè)初始震源時(shí)間就是所有估算震源時(shí)間的平均值。只有一個(gè)S-P觀測(cè)值的情況下，觀測(cè)值的情況下，vp/vs被設(shè)置為被設(shè)置為 。如果沒(méi)有。如果沒(méi)有S-P時(shí)間的時(shí)間的觀測(cè)值，震源時(shí)間就設(shè)置為最早觀測(cè)到的起始時(shí)間。觀測(cè)值，震源時(shí)間就設(shè)置為最早觀測(cè)到的起始時(shí)間。 0代表代表1970年年1月月1日，

30、也得到了正確結(jié)果（昨天為什么不日，也得到了正確結(jié)果（昨天為什么不對(duì)？）對(duì)？）HYPOSAT: hyposat-parameter中中已由已由MSDP設(shè)定的部分參數(shù)設(shè)定的部分參數(shù) CONFIDENCE LEVEL(置信水平置信水平)：以百分比表示的不確定性：以百分比表示的不確定性輸出的可信度，默認(rèn)的不確定性為輸出的可信度，默認(rèn)的不確定性為一個(gè)標(biāo)準(zhǔn)偏差一個(gè)標(biāo)準(zhǔn)偏差(即約即約68.3)。在這里已固定為在這里已固定為90%，在參數(shù)文件中的賦值不起，在參數(shù)文件中的賦值不起作用。作用。 CONFIDENCE LEVEL (68.3 - 99.99 %) : 90. Confidence level of given uncertainties: 90.00 %Source time : 2009 05 04 09 16 9.473 +/- 0.154 s or 1241428569.473 +/- 0.154 sEpicenter lat: 23.5326 +/- 0.0090 degE