Redetermination of Focal Depth in Hutubi Region of the Middle and East of North Tianshan in Xinjiang
-
摘要: 联合Hyposat法、PTD法和gCAP矩张量反演法,重新测定新疆北天山中东段呼图壁地区2010—2017年502个地震的震源深度,并对震源深度剖面进行初步分析。结果表明,重新测定的震源深度优势分布为15—20km,平均震源深度为16km,呼图壁MS 6.2地震的震源深度为20km;研究区南部和中部的震源深度集中分布在20km左右,与北天山壳内低速体的层位相当,可能是上地壳和下地壳之间的韧性剪切带存在的部位,起到滑脱层的作用,研究区北部的震源深度则向浅部扩展;呼图壁MS 6.2地震的发震断裂可能在清水河子断裂下方的1条隐伏反冲断层上,可能是霍尔果斯断裂向前沿断坡冲断受阻而在相反方向上发育分支反冲断层的结果。
-
关键词:
- 呼图壁地区 /
- 震源深度 /
- Hyposat方法 /
- PTD方法 /
- gCAP矩张量反演法
Abstract: Combined with Hyposat method, PTD method and gCAP moment tensor inversion method, the focal depth of 502 earthquakes in Hutubi region during 2010—2017 are re-determined, and the focal depth profile is preliminarily analyzed. The results show that the dominant distribution of focal depth is 15—20km, the average focal depth is 16km, and the focal depth of Hutubi MS 6.2 earthquake is 20km. The focal depth in the southern and central regions of the study area is concentrated around 20km, the depth is similar to that of the low velocity layer in the 20km crust of the North Tianshan Mountains, which may be the ductile shear zone between the upper crust and the lower crust. The location of the shear zone acts as a slip layer, and the focal depth in the northern part of the study area extends to the shallow part. The causative fault of the Hutubi MS 6.2 earthquake may be on a hidden recoil fault beneath the Qingshuihezi fault and may be the result of the thrust of the Horgos fault to the front slope and the development of the branch recoil fault in the opposite direction.-
Key words:
- Hutubi region /
- Focal depth /
- Hyposat method /
- PTD method /
- gCAP moment tensor inversion method
-
引言
陆地与大气紧密相连,地震前大气的异常情况受到地震学者们的重视。前人进行了大量探索性研究工作,以寻求地震前兆(徐国钧等,1993;李贵福等,1996;曹新来等,1997;张铁宝等,2013)。在这些地震研究中,长期连续完整且均一性较好的气象资料是研究地震前气温异常的数据基础。但是,由于各种原因(如环境干扰、硬件故障等),历史气温资料缺测现象时有发生,导致气温观测资料缺测,造成历史资料的不连续(王建国等,2010;2013;姚会琴等,2012)。
中国许多学者开展了对日、月、年时间尺度的气象资料缺测插补研究,并利用一维车贝雪夫多项式展开、线性回归、标准序列法、基于SVD的迭代等方法对中国部分地区的气象日、月、年值资料进行了恢复性实验(张秀芝等, 1996a, 1996b;涂诗玉等,2001;张永领等,2006;王海军等,2008;余予等,2012),但对整点气温值进行缺测插补鲜见文献报道。在国外,Huth等(1995)建立回归模型来插补缺测的日气温数据,Eischeid等(2000)采用空间差值法,插补后建立美国西部40年逐日气温和降水数据集,但这些插补模型只用于1个或数个缺测日数据的插补,不合适用于连续几个月数据缺测的情况。DeGaetano等(1995)引用改进的标准序列法,对美国东北部近400个站的日最高、最低气温缺测值进行了插补。标准序列法和线性回归法解决了插补长期连续缺测数据的问题。整点气温值的缺测插补方法可借鉴日平均、最高、最低气温值的缺测插补方法。由于线性回归法具有更好的统计性能和稳健性,本文对线性回归法进行改进,考虑了距离因素,采用线性回归模型参数求解法,解决了连续数日甚至数月造成的气温缺测问题,为地震前后气温变化特征的研究提供长期连续完整可靠的数据资料,也为今后开展卫星遥感红外亮温与卫星过境时刻气温的对比研究奠定了数据基础。
1. 研究区及数据资料
1.1 研究区介绍
本文选取的研究区范围为37°—42°N、113°—119°E,在研究区内收集到15个地震观测站(气温观测站)的气温整点值数据,观测站分布情况如图 1所示。
1.2 气温资料
本文收集了15个地震观测站气温数据,气温指地面以上1.5m处百叶窗测得的空气温度。太阳的热能被地面吸收后,地面再通过辐射、传导和对流把热传给空气,这是空气中热量的主要来源。气温的观测范围-30℃—70℃,精度0.1℃。由于各观测站安装时间、停测时间不同,导致观测时间的长度不等,多数观测数据起止时间为2007年1月1日,截止时间为2014年12月31日。
1.3 气温数据的缺测情况
由于观测环境、仪器设备故障等原因,部分观测站的观测数据缺测,包括仅缺测1个值、缺测1天的值(即24个观测值)、连续缺测几天甚至1个多月导致长达数千个观测值连续缺测等情况。基于震例研究,本文仅讨论唐山站气温数据的缺测情况(表 1),可以看出该站缺测情况较严重,如2008年8月18日—9月15日连续缺测696个整点值,2009年4月2日—5月13日连续缺测1000个整点值,2012年6月2日—24日连续缺测552个整点值,2013年1月21日—2月5日连续缺测384个整点值。唐山站数据连续缺测时间较长,连续缺测几十个值的频率较高,同时该站还存在错误值,如2010年11月22日18时的观测气温整点值是59.701,类似的错误值在每年均有出现。2010、2011年的观测数据较完整。
表 1 唐山站整点气温数据的缺测统计Table 1. Missing data in integral point temperatures from the Tangshan site数据起止时间 年份 数据缺测情况 2007-01-01 — 2013-12-18 2007 7月14—15日、7月25日 2008 3月8—11日、11月17—24日、11月26日、12月3日、12月7—11日 2009 3月18日、3月26日、4月2日—5月13日、6月27—29日、7月25日、11月10日—11日8时、11月19—21日、11月25—27日、11月30日、12月1—5日、12月13日 2012 6月1—24日、8月30日—9月1日、9月8日、12月24—29日 2013 1月21日—2月5日、6月25日、8月1—3日、7月27—30日、9月5—6日、9月23—24日 1.4 气温数据的日变情况
图 2为唐山站2012年12月17日—22日气温整点值变化曲线,可以看出气温日变明显,最高温一般出现在14时左右,最低温出现在8时左右,符合日变规律。
1.5 气温数据的月变情况
对唐山站2012年1、4、7、10月的日均值进行分析,绘制了相关曲线,如图 3所示。由图可以看出,气温变化整体表现出明显的夏高冬低的年变规律;1月、7月气温变化平稳,4月气温表现为升温过程,10月气温表现为明显的降温过程。
1.6 气温数据的年变情况
唐山站2008—2013年连续的气温整点值数据变化曲线如图 4所示。由图可以看出,气温表现出明显的夏高冬低的年变特征,夏季最高日平均温度接近31℃,冬季最低日平均温度接近-14℃;气温在冬季和夏季处于稳定的状态,气温从3月开始上升,1月—7月处于升温过程,气温上升快,7月、8月气温达到最高,9月后气温开始迅速降低,进入降温阶段,符合季节变化规律;此外,数据缺测明显。
2. 研究方法
应用线性回归法解决连续几日甚至数月的整点气温值缺测问题,并采用交叉验证方法对插补结果进行误差分析。
2.1 线性回归模型的建立
本文对线性回归方法进行改进,考虑了参考站和缺测站之间的距离。改进后该方法更科学,可以更好地去除距离因素的影响。
利用邻近站资料对距离进行加权,建立回归模型,插补缺测站资料的方程式为:
$$ {\hat y_i} = \frac{{{a_{1i}}{x_{1i}}}}{{{d_{1i}}}} + \frac{{{a_{2i}}{x_{2i}}}}{{{d_{2i}}}} + \cdots + \frac{{a{}_{mi}{x_{mi}}}}{{{d_{mi}}}} + {a_{m + 1}} $$ (1) 其中,${\hat y_i}$为插补值,${x_{mi}}$为临近站数据,${a_{mi}}$为回归模型的参数,${d_{mi}}$为邻近站与缺测站之间的距离,m为临近站站数。利用最小二乘原理求解回归模型的参数,即使观测值和插补值之间差值的平方和(Q)最小:
$$ \min Q = {\rm{min}}\sum\limits_{i = 1}^n {{{({y_i} - \hat y{}_i)}^2}} $$ (2) 其中${y_i}$为插补站观测值。
因历史同期各要素时空变化规律通常比较相似,选择缺测整点值前后若干整点值的历史同期(不包括缺测值所在的年份)数据,作为拟合回归模型的样本数据,建立线性回归模型,并利用附近站资料,计算缺测记录插补值。
2.2 误差检验方法
本文采用交叉验证的方法对缺测记录的插补结果进行分析,即假设某个站的记录缺测,首先利用插补模型插补整点气温数据,然后对插补值与实际观测资料进行对比和误差分析,并用平均绝对误差(Mean Absolute Error,MAE)代表插补精度(王海军等,2008),其表达式为:
$$ MAE = \frac{1}{N}\sum\limits_{i = 1}^N {\left| {{x_{io}} - {x_{ei}}} \right|} $$ (3) 其中,${x_{oi}}$为第i个实际整点观测值,${x_{ei}}$为第i个插补整点值,N为插补整点值的个数。
3. 插补结果分析
3.1 邻近参考站及时间窗选择
参考站的选择不仅与观测站密度有关,也与插补站及其邻近站所处的地理环境有关(如平原、丘陵、山区等)。同时,时间窗大小也对缺测数据的插补精度有影响。本文采用滑动优选法确定时间窗,时间窗的宽度为气温整点值个数,高度为年数。以选择缺测整点值所在的年份为中心,其前后若干点值历史同期若干年的数据作为样本数据,对于前后无资料的年份,则使用靠近插补年份的资料。
唐山站位于华北平原,周围地势平坦,气温变化相近,故采用最短距离的原则选取临近参考站。根据距离及地形因素,选取了该站周边的北京、昌黎、蓟县、宁河、青光和徐庄子6个观测站,年数为7年。采用15个整点值作为样本资料,建立线性回归模型,插补唐山站的缺测值(包括连续和不连续的单点缺测值)。插补站及其邻近站的基本信息见表 2。
表 2 唐山插补站及其邻近台站信息Table 2. Information of interpolation site and its neighboring sites台站 经度/°E 纬度/°N 海拔/m 与唐山站距离/km 唐山 116.6 40.4 30 0 北京 116 40 60 156 昌黎 119 39.7 18.4 83 蓟县 117.5 40.1 65 82 宁河 117.7 39.4 2.5 45 青光 117 39.2 3 110 徐庄子 117.2 38.7 2 142 3.2 误差检验
采用交叉验证的方法评估上述插补方法,统计了2010年3月的31天实际观测整点值与相应插补值的相关系数,统计结果见表 3。从表中可以看出,3月8日、14日、15日、19日和20日的相关系数相对较低。王海军等(2008)经过对比研究,在平原地区选取了4个参考站,选取年数为8年、天数为15天,插补误差最小。唐山站也位于平原地区,考虑与唐山站的距离及地形因素,选取昌黎、宁河、青光、徐庄子4个参考站,并选取年数为7年、15个整点值的优化模型。通过对比发现,4个参考站的相关系数偏高(表 3)。
表 3 观测值与插补值相关系数Table 3. Correlation coefficients between interpolated and observed data日期 优化模型相关系数 相关系数 2010-03-01 0.922018 0.916277 2010-03-02 0.962184 0.967497 2010-03-03 0.934722 0.94439 2010-03-04 0.986854 0.973136 2010-03-05 0.97029 0.931489 2010-03-06 0.925587 0.958141 2010-03-07 0.962133 0.96263 2010-03-08 0.72234 0.652517 2010-03-09 0.896098 0.915444 2010-03-10 0.980709 0.977074 2010-03-11 0.928618 0.952106 2010-03-12 0.892358 0.900113 2010-03-13 0.950986 0.959182 2010-03-14 0.799879 0.753849 2010-03-15 0.881494 0.792166 2010-03-16 0.974627 0.968444 2010-03-17 0.972681 0.970844 2010-03-18 0.976031 0.973557 2010-03-19 0.596366 0.618699 2010-03-20 0.846827 0.745552 2010-03-21 0.951358 0.951009 2010-03-22 0.932557 0.951905 2010-03-23 0.901255 0.919988 2010-03-24 0.897624 0.920017 2010-03-25 0.974534 0.981077 2010-03-26 0.958726 0.979361 2010-03-27 0.976428 0.948801 2010-03-28 0.987332 0.993851 2010-03-29 0.927046 0.962993 2010-03-30 0.841163 0.902429 2010-03-31 0.968156 0.977347 利用式(3)对优化后的模型得到的插补结果计算平均绝对误差,并统计其误差的比例分布(表 4)。从表 4可以看出,插补误差在±0.5℃范围内的比例为60.5%,在±0.8℃范围内的比例为80.6%,其误差绝对值大于1.0℃的为9.4%,平均绝对误差为0.82℃。
表 4 唐山站整点气温缺测插补误差比例Table 4. Proportion of error for interpolated integral point temperature误差范围 比例/% (-∞,-1) 7.8 [-1,-0.8) 5.4 [-0.8,-0.5) 9.5 [-0.5,0) 31.2 [-0.5,0.5] 60.5 [0,0.5) 29.3 [0.5,0.8) 11.6 [0.8,1) 3.6 [1,∞) 1.6 [-0.8,0.8] 81.6 3.3 插补结果
针对唐山站2008年1月1日—2013年12月18日的气温整点值缺测数据及错误数据,利用唐山邻近站的同期数据和线性回归模型,对缺测数据插补完整,并修正错误数据,绘制气温整点值的年变曲线,如图 5所示。从图中可以看出,在长时间序列缺测的部位,插补值与前后正确的数据衔接吻合,没有出现突升或突降变化。插补后完整连续的数据符合夏高冬低的年变规律,气温6年的变化形态一致。
4. 震例应用
据中国地震台网中心测定,2012年5月28日10时22分在河北省唐山市辖区、滦县交界处发生4.8级地震,震源深度8km。
利用插补完整的连续数据,分析2012年3—5月震前气温数据的变化情况。首先,选取无震年份(2008—2011年)同期(3—5月)整点值气温,并计算历年同期气温日均值,以此作为3—5月正常的背景值;其次,将2012年3—5月日均值与历年同期背景值做差值,得到2012年3—5月份数据与历年同期均值的偏移程度;最后,以无震年份同期气温波动范围作为基准,即2008—2011年3—5月所有气温值的标准差作为判断标准。根据以往经验将差值大于2倍标准差视为气温前兆异常。
从2012年3—5月当年日均值与历年同期(2008—2011年)多年日均值的差值及标准差(图 6)中可以看出,从3月27日开始,唐山站数据出现大幅度增温异常现象,4月增温天数也较多,5月1日—11日的气温日均值仍然高于历年同期,特别是5月10日(即震前2天)增幅达到约8℃,且大于2倍标准差;震后差值开始变小,并逐步恢复至平静。
5. 结论
长期连续完整的历史气象资料是震前气温异常判别研究的重要数据基础,但由于观测环境、仪器故障等原因,造成气温观测数据缺测或错误数据,且部分数据缺测的时间较长。为此,本文利用线性回归模型,插补缺测和错误的气温整点值数据,较好地解决了长期连续缺测的情况。
通过对唐山观测站2008年1月1日—2013年12月18日的气温整点值缺测数据及错误数据进行插补,使得数据完整连续,并应用插补完整的气温整点值数据,分析研究了2012年5月28日唐山4.8级地震的气温前兆异常现象,主要得出以下结论:
(1)唐山观测站的插补值与其前后的观测数据衔接吻合,插补后完整连续的数据符合夏高冬低的年变规律。
(2)插补误差在±0.5℃范围内的比例为60.2%,在±0.8℃范围内的比例为80.3%,其误差绝对值大于1.0℃的为9.6%,平均绝对误差为0.84℃。插补值与观测值的相关系数大部分在0.9以上,可见插补结果真实可靠。
(3)从3月27日起,唐山观测站数据出现增温异常,震前2天增温幅度约8℃。
-
-
陈向军, 上官文明, 宋秀青等, 2014.新疆全区和分区地壳速度模型的分析.中国地震, 30(2):178-187. doi: 10.3969/j.issn.1001-4683.2014.02.005 邓起东, 冯先岳, 张培震等, 1999.乌鲁木齐山前坳陷逆断裂-褶皱带及其形成机制.地学前缘, 6(4):191-201. doi: 10.3321/j.issn:1005-2321.1999.04.001 高原, 周蕙兰, 郑斯华等, 1997.测定震源深度的意义的初步讨论.中国地震, 13(4):321-329. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199700698774 黄禄渊, 张贝, 瞿武林等, 2017.2010年智利Maule特大地震的同震效应.地球物理学报, 60(3):972-984. http://www.cnki.com.cn/Article/CJFDTotal-DQWX201703013.htm 孔祥艳, 陈向军, 赵石柱等, 2016.新疆数字测震台网"九五"与"十五"运行期间震源深度对比分析.地震地磁观测与研究, 37(2):67-71. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=dzdcgcyyj201602012 李赫, 张辉, 王熠熙等, 2017.2012年5月28日唐山MS 4.8地震震源机制解和震源深度精确测定.地震, 37(3):107-116. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=hbdzkx201701007 李艳永, 热依木江, 唐明帅等, 2016.利用震相方位角改善地震台网稀疏地区地震定位精度.地震地磁观测与研究, 37(2):57-62. http://d.old.wanfangdata.com.cn/Periodical/dzdcgcyyj201602010 李莹甄, 高国英, 高歌等, 2008.新疆乌鲁木齐地区震源深度分布与断层关系研究.地震地质, 30(2):534-543. doi: 10.3969/j.issn.0253-4967.2008.02.018 李志伟, 黄志斌, 王晓欣等, 2015.USGS地震目录中4-5级震源深度异常地震可靠性初步研究:以南北地震带若干地震为例.地球物理学报, 58(4):1236-1250. http://www.cqvip.com/QK/94718X/201504/664747561.html 刘建明, 王琼, 李金等, 2018.2016年12月8日呼图壁MS 6.2地震序列重定位与发震构造.地震地质, 40(3):566-578. http://d.old.wanfangdata.com.cn/Periodical/dzdz201803005 刘永梅, 刘芳, 刘改梅等, 2017.蒙宁交界区地震震源深度的对比分析.中国地震, 33(2):290-300. doi: 10.3969/j.issn.1001-4683.2017.02.010 罗艳, 曾祥方, 倪四道, 2013.震源深度测定方法研究进展.地球物理学进展, 28(5), 2309-2321. http://d.old.wanfangdata.com.cn/Periodical/dqwlxjz201305013 吕作勇, 杨传成, 房立华, 2014.广东地震台网近震定位精度的经验性GT准则.震灾防御技术, 9(增刊):657-664. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=2014s112&journal_id=zzfyjs 邵学钟, 张家茹, 范会吉等, 1996.天山造山带地壳结构与构造:乌鲁木齐-库尔勒地震转换波测深剖面.地球物理学报, 39(3):336-346. doi: 10.3321/j.issn:0001-5733.1996.03.007 宋秀青, 缪发军, 刘双庆等, 2014.2014年新疆于田MS 7.3地震序列的震源深度测定.中国地震, 30(2), 198-207. http://d.old.wanfangdata.com.cn/Periodical/hndz201404003 宋秀青, 2017.用PTD方法测定2017年精河MS 6.6地震序列的深度.中国地震, 33(4):694-702. 王椿镛, 楼海, 魏修成等, 2001.天山北缘的地壳结构和1906年玛纳斯地震的地震构造.地震学报, 23(5):460-470. doi: 10.3321/j.issn:0253-3782.2001.05.002 王海涛, 李志海, 赵翠萍等, 2007.新疆北天山地区MS ≥ 2.0地震震源参数的重新测定.中国地震, 23(1):47-55. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=zgdz200701006 王晓楠, 唐方头, 邵翠茹, 2018.南迦巴瓦构造结周边地区主要断裂现今运动特征.震灾防御技术, 13(2):267-275. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=20180202&journal_id=zzfyjs 王周元, 1987.测震定位中模型影响的数值实验.西北地震学报, 9(3):63-70. http://www.cnki.com.cn/Article/CJFDTotal-ZBDZ198703009.htm 谢卓娟, 吕悦军, 彭艳菊等, 2008.渤海海域地震震源深度的分布特征.震灾防御技术, 3(3):311-320. doi: 10.3969/j.issn.1673-5722.2008.03.014 杨晓平, 顾梦林, 孙振国等, 2002.1906年新疆玛纳斯大震区的多层次逆冲构造与深部结构.地震地质, 24(3):303-314. doi: 10.3969/j.issn.0253-4967.2002.03.002 张国民, 汪素云, 李丽等, 2002.中国大陆地震震源深度及其构造含义.科学通报, 47(9):663-668. doi: 10.3321/j.issn:0023-074X.2002.09.004 张培震, 邓起东, 徐锡伟等, 1994.盲断裂、褶皱地震与1906年玛纳斯地震.地震地质, 16(3):193-204. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=QK199400054115 张志斌, 朱皓清, 李艳永, 2015.基于地震波反演研究南天山中西段的震源深度.震灾防御技术, 10(增刊):701-711. http://zzfy.eq-j.cn/zzfyjs/ch/reader/view_abstract.aspx?flag=1&file_no=2015s101&journal_id=zzfyjs 张志斌, 朱皓清, 李艳永, 2018.新疆数字地震台网地方性震级量规函数的初步研究.中国地震, 34(4):667-674. doi: 10.3969/j.issn.1001-4683.2018.04.007 朱元清, 石耀林, 李平, 1990.一种确定地震深度的新方法.地震地磁观测与研究, 11(2):4-12. http://www.cqvip.com/qk/96509X/199002/228831.html 胥颐, 刘福田, 刘建华等, 2000.天山地震带的地壳结构与强震构造环境.地球物理学报, 43(2):184-193. doi: 10.3321/j.issn:0001-5733.2000.02.006 Bondár I., Myers S. C., Engdahl E. R., et al., 2004. Epicentre accuracy based on seismic network criteria. Geophysical Journal International, 156(3):483-496. doi: 10.1111/gji.2004.156.issue-3 Lu R. Q., He D. F., Xu X. W., et al., 2017. Seismotectonics of the 2016 M 6.2 Hutubi earthquake:Implications for the 1906 Manas M 7.7 Manas earthquake in the northern Tian Shan belt, China. Seismological Research Letters, 89(1):13-21. http://d.old.wanfangdata.com.cn/Periodical/dqwlxb201810011 Schweitzer J., 2001. HYPOSAT- an enhanced routine to locate seismic events. Pure and Applied Geophysics, 158(1-2):277-289. http://cn.bing.com/academic/profile?id=d5ce922cdc793927cd4dcba2fad33d43&encoded=0&v=paper_preview&mkt=zh-cn Stein S., Wiens D. A., 1986. Depth determination for shallow teleseismic earthquakes:Methods and results. Reviews of Geophysics, 24(4):806-832. doi: 10.1029/RG024i004p00806 Zhao L. S., Helmberger D. V., 1994. Source estimation from broadband regional seismograms. Bulletin of the Seismological Society of America, 84(1):91-104. http://d.old.wanfangdata.com.cn/NSTLQK/NSTL_QKJJ027969581/ Zhu L. P., Helmberger D. V., 1996. Advancement in source estimation techniques using broadband regional seismograms. Bulletin of the Seismological Society of America, 86(5):1634-1641. http://cn.bing.com/academic/profile?id=b303f16dd343f999802dfd2d2577b5ae&encoded=0&v=paper_preview&mkt=zh-cn Zhu L. P., Ben-Zion Y., 2013. Parametrization of general seismic potency and moment tensors for source inversion of seismic waveform data. Geophysical Journal International, 194(2):839-843. doi: 10.1093/gji/ggt137 -