四川省CORS网地壳运动速度场模型研究

使用GAMIT/GLOBK解算2006~2015年覆盖四川地区的50个CORS站的观测数据，并提取坐标时间序列，获得速度场模型。对四川省整体地壳运动以及区域应变场进行分析，结果表明，四川省CORS站水平方向平均速度为38.72 mm/a，速度场的优势方向为S73.9°E；欧亚框架下水平方向平均速度为14.40 mm/a，优势方向为S77.78°E，较ITRF2008框架下的速度降低了24.32 mm/a；垂直方向平均速度为3.43 mm/a，整体上表现为隆升状态，位于东部平原的少部分测站处于下降趋势，可能是由于城市建设、地下水或者煤矿、石油等自然资源的过度采伐导致；最大主应变率达到7.31×10-11/a，最小主应变率为-5.4×10-11/a，其中区域地壳面膨胀率高值区域集中在川西，低值区域和最大剪切应变率高值区域分布在三大断裂带附近。     

环渤海区域及邻区现今地壳构造运动形变特征分析

将环渤海区域划分为5个次级块体，利用整体旋转与线性应变模型分析2009~2014年中国大陆构造环境监测网络的GPS对地观测速度场数据，得到相对于欧亚板块0.5°×0.5°的水平运动速度场和应变场。分析形变速度场和应变场的空间变化发现，环渤海区域整体上呈现0.25×10^-9/a的NW-SE 111.3°的双向趋势性扩张运动，太行块体、山西块体分别呈NW-SE 116.3°、NW-SE 130°的微弱扩张运动，胶辽块体、冀鲁块体和阴山-燕山块体分别存在NW-SE 144.3°、NE-SW 39.5°和NW-SE 155.6°的微弱压缩运动，其中太行块体的扩张运动量级相对于其他次级块体较大。尽管环渤海区域内部各次级块体的形变在空间上存在以上差异，但整体上环渤海区域主压应变轴方向基本为NEE-SWW，主张应变轴方向为NNW-SSE，且主压应变轴走向为NE 47.70°～89.74°，与利用地球物理方法得到的主压应变轴优势方向大体一致，表明环渤海区域现今地壳运动相对稳定。     

青藏高原西部及邻区构造应力场反演和断层滑动趋势分析

基于Global CMT提供的1976-01-01~2015-12-31期间1 255个M_W≥3.5地震的震源机制解,利用FMSI方法对青藏高原西部及邻区（新疆西部地区）应力场进行反演,得到区域应力场特征,并分析西昆仑断裂带及阿尔金断裂带的断层滑动趋势。结果表明,青藏高原西部及邻区的最大主压应力轴σ1整体呈NNE-SSW向,自西向东顺时针旋转;最小主压应力轴σ3整体呈NWW-SEE向。最大主压应力轴σ1倾角大小约为5°~9°,说明该地区的构造应力场以水平运动为主。西昆仑断裂带及阿尔金断裂带具有较强的滑动趋势（T_s≥0.5）,说明西昆仑断裂带及阿尔金断裂带未来失稳滑动的可能性较大。     

2009年陕西高陵ML4.8地震震源参数及构造特征分析

利用双差定位法对2009-11-05陕西高陵M_L4.8地震及其余震序列进行重定位研究,并采用CAP方法求解主震震源机制解,根据所得结果综合分析地震构造特征。结果表明,高陵地震序列主震深度约为6 km,余震主要分布在2~12 km深度范围内,地震序列在平面上呈簇状分布,在深度上大致呈垂直分布;主震震源机制解节面Ⅰ走向103°,倾角61°,滑动角-74°;节面Ⅱ走向252.40°,倾角32.78°,滑动角-116.44°;最佳矩心深度为7 km,矩震级近似为4.2级,P轴方向为北东方向;主震为正断兼少量走滑特征。综合分析认为,渭河断裂与渭南-泾阳断裂交会处为重要的地震潜在地段。     

山西地区地壳S波速度结构

利用山西及周边共44个宽频地震台2014-01~07地震噪声数据,反演瑞利波群速度和相速度频散,得到研究区域的三维剪切波速度。上地壳中,山西断陷带表现为较低的速度异常,而两侧的隆起区则显示高速异常;中下地壳主要表现为以北纬38.5°左右为界的南高北低的速度特征。大同盆地从浅部到中下地壳整体呈现为低速特征,与该地区广泛存在的火山群有关。     

青藏高原东北缘及其周边地区现今构造应力场及地震活动性

收集青藏高原东北缘及其周边地区371个震源机制解数据,对该区域进行1°×1°经纬度网格的构造应力场反演。结果表明,青藏高原东北缘内缘最大主压应力以NE向为主,外缘最大主压应力由西向东呈现出顺时针旋转的特征,即最大主压应力轴由西段近NNE向转为中段NE向,再转为东段NWW向。从不同深度对青藏高原东北缘及其周边地区进行2°×2°经纬度网格的构造应力场反演。结果显示,地壳上下部分断层体系存在巨大的变化差异,产生这种差异的原因可能与地壳的挤压增厚有关。     

安徽省旅游发展空间错位的模型分析

旅游资源与旅游区位的地域差异是客观存在的,对区域旅游发展空间错位及其原因的分析可为促进区域旅游均衡发展提供科学依据。本文以安徽省16个地级市为研究单元,运用重力模型和二维组合矩阵,选取旅游收入、资源丰度和旅游区位作为评价指标,定量分析旅游发展空间错位关系,并利用ArcGIS软件将结果进行可视化表达。研究发现:安徽省旅游收入、资源丰度和旅游区位的重心分别处于(117.63°E,31.18°N)、(117.51°E,31.12°N)和(117.20°E,32.00°N),经度方向最大错位0.43°,最小0.12°,纬度方向最大错位0.88°,最小0.06°;与区域几何中心相比,旅游收入重心和资源丰度重心均向南偏移;从旅游收入-资源丰度、旅游收入-旅游区位的组合矩阵来看,呈错位和同步发展的城市均为8个。整体而言,安徽省旅游空间错位主要是因旅游区位和旅游资源的条件差异,四大旅游板块要结合自身的特点,发挥优势,挖掘潜力,通过采取不同的空间矫正策略达到区域旅游协同发展的目的。     

使用震源机制解反演内蒙古东部地区构造应力场

使用内蒙古地震台网2008-01~2017-12记录的数字地震波形观测资料,采用基于P波初动和振幅比联合求解（FOCMEC）方法计算震源机制,获得内蒙古东部地区发生的138个ML≥2.5地震的震源机制解。对震源机制解特征进行分析显示,研究区地震震源类型以走滑型为主,兼有一定比例的逆冲型,而正断型相对较少,P轴和T轴走向的优势方向分别为NEE-SWW向和NNW-SSE向。将研究区划分为1°×1°的网格,采用阻尼区域应力反演（MSATSI）方法获得研究区的应力张量和应力场图像。应力场反演结果显示,内蒙古东部地区应力状态主要表现为走滑,局部区域存在剪切、挤压,最大主应力方向主要分布在NEE-SWW向,由北向南顺时针旋转。     

2019年四川长宁MS6.0地震孕震机理分析

基于四川省、贵州省、重庆市数字地震台网的三分量宽频带波形资料和震相报告数据,使用Geiger定位法和ISOLA反演方法,对长宁M_S6.0地震序列进行重定位,并反演序列中5次M_S≥5.0地震震源机制解。重定位结果显示,地震序列主要分布在白象岩-狮子滩背斜轴部,总体沿NW65°方向展布,沿该方向震源深度特征为北西深、南东浅。震源机制解反演结果显示,主震破裂表现为高倾角左旋走滑特征,节面Ⅰ走向/倾角/滑动角为20°/72°/151°,节面Ⅱ走向/倾角/滑动角为120°/62°/21°,与白象岩-狮子滩背斜轴部走向基本一致。5次地震纯双力偶（double-couple,DC）成分含量均小于70%,矩心深度均小于10 km。综合分析认为,此次长宁M_S6.0地震发震断层为发育在白象岩-狮子滩背斜的左旋走滑断层,震因为流体注入和复杂的地质构造活动。     

近震波形反演震源机制解方法对比研究

分别利用近震波形反演求解震源机制解的gCAP方法和全波形反演方法反演了鲁甸MS6.5地震的震源机制解,并比较两种方法的优点与特性。结果发现,数据方位角的覆盖范围会影响解的稳定性,在使用gCAP方法时尤其需要注意;两种方法获取的震源机制解结果较为一致,走向为73°~77°,倾角为72°~85°,滑动角为157°~180°。     

Sentinel-1 SAR数据约束的2018年谢通门MW5.8地震同震破裂模型

通过Sentinel-1卫星升降轨数据获取谢通门地震的同震形变场,并基于均匀弹性半无限位错模型反演地震的同震滑动分布模型。InSAR同震形变场表明,升降轨视线向最大形变量分别为0.049 m和0.051 m,形变场长轴大致呈南北方向,位于甲岗-定结断裂西侧。通过对倾角和倾向进行格网搜索发现,西倾节面更可能为该地震的发震节面。反演结果表明,滑动分布主要位于2~10 km深度范围内,平均滑动量为0.02 m,最大滑动量为0.10 m,发震断层倾角为47°,平均滑动角为-81.60°,显示该地震以正倾滑动为主。大地测量数据约束的该地震震中为30.27°N、87.75°E,震源深度为6.58 km,释放地震矩为5.056×10¹⁷ Nm,对应矩震级为M_W5.7,与GCMT、USGS公布的震级基本一致。综合分析震中位置和滑动机制认为,甲岗-定结断裂的分支断层为本次谢通门地震的发震断层。     

东北地区近期水平形变应变场研究

根据2011年和2013年东北地区GPS观测资料计算得到的地壳水平运动场和连续主应变场,研究了东北地区地壳运动特征和应变应力状态, 求解块体运动参数,计算东北主要断裂带的错动速率和断层应变应力状态,分析断裂带活动性。研究表明,2011～2013年期间,东北地区东部运动规律性相对较好,且东部运动场收敛于珲春-龙井一带,向东位移;中部和西部运动场一致性较弱,位移方向相对分散。主应变反映出,东北地区2011～2013年主压应力轴优势方向为NE55°。两个水平主应变轴都缩短的挤压区分布面积较小,两个主应变轴都伸长的拉张区次之,最小主应变轴缩短的区域和最大主应变轴伸长的区域所占比重最大。该时段Ms≥5.0地震集中在2013年,地震多发生在运动方向转折带附近及挤压-剪切活动区,最大的5.8级地震发生在不同运动区交界带和挤压作用相对强烈的构造带。     

2021年滦州MS4.3地震和2015年昌黎MS4.2地震震源参数及构造意义研究

基于京津冀地震台网观测资料,利用CAP方法反演滦州M_S4.3地震、昌黎M_S4.2地震及其M_S3.3余震的震源机制,并利用近震深度震相获得更为准确的震源深度,结合双差定位法获得2个地震序列的震源分布结果,对发震构造及成因关联开展分析研究。结果显示：1）滦州M_S4.3地震的节面Ⅰ走向、倾角、滑动角分别为211°、85°、168°,节面Ⅱ分别为302°、78°、5°,震源错动类型为走滑型,震源深度为8 km,地震序列的震源分布呈NNE向,短轴剖面显示断层面倾角近垂直,认为其发震断层面为节面Ⅰ;昌黎M_S4.2地震及M_S3.3余震的节面Ⅰ分别为189°、68°、161°及190°、61°、170°,节面Ⅱ分别为286°、72°、23°及285°、81°、29°,震源错动类型同为走滑型,震源深度为10.5 km,地震序列的震源分布呈NNE向,短轴剖面显示断层面倾角近垂直,认为其发震断层面为节面Ⅰ;2）基于滦州M_S4.3地震、昌黎M_S4.2地震及滦州M_S7.1地震的震源参数结果,结合区域地质构造等资料分析认为,3次地震的发震构造不是上地壳先存断裂,而可能与震源区的深部构造背景密切相关,即壳内包体现象是孕育这些地震的共同基础。     

滇西北三江并流地区地壳S波速度结构北大核心CSCD

对滇西北三江并流地区36个台站记录的远震波形数据进行处理,提取径向P波接收函数,基于中国地震科学实验场地区地壳剪切波速模型获取不同区域的初始模型,采用时间域线性反演方法和bootstrap重采样技术,反演各台站下方的S波速度结构。结果表明:1)研究区地壳S波速度结构整体以澜沧江断裂带和27.5°N为界,表现出明显的区域差异。在上地壳,S波速度在0~25 km深度范围内存在由地表向壳内延伸的低速层,但不同区域低速层厚度存在较大差异,27.5°N以北低速层较为发育,厚度在20~25 km之间;27.5°N以南,澜沧江断裂带东西两侧低速层厚度差异明显,断裂带以西在10~25 km之间,断裂带以东约4 km。2)在下地壳,S波速度在横向上以澜沧江断裂带和27.5°N为界,表现出NE与SW区域以相对低速为主,NW和SE区域以相对高速为主的分布特征;垂向上,怒江断裂带中南部和金沙江-红河断裂带北端附近局部区域25~45 km深度范围内存在低速体,与上地壳低速层相连。3)澜沧江断裂带为地壳厚度和S波速度变化的重要分界面,对滇西北三江并流地区的构造演化过程具有控制作用。     

走滑断裂活动导致地应力解耦的机理研究——以龙门山断裂带东北段为例

2008年汶川MS8.0地震后,沿龙门山断裂带系列水压致裂地应力测量结果表明,断裂带东北段(茂县-绵竹连线东北侧)北川、平武及广元等地浅表层最大主应力方向与区域应力场相比发生了25°~69°不等的偏转,亦称地应力解耦。初步分析认为,龙门山断裂带东北段的右旋走滑活动直接导致了区域内最大主应力方向的解耦,且最大主应力方位向偏转角度与断裂右旋走滑活动量正相关。为合理解释这一现象,探讨走滑断裂活动导致地应力解耦的机理,本文采用二维有限元数值模拟方法,选取汶川MS8.0地震断裂(龙门山中央断裂)上的同震右旋走滑位移作为边界条件,首先模拟了龙门山断裂带东北段右旋走滑活动产生的扰动应力特征,之后结合断裂带附近典型钻孔地应力测量结果,叠加该扰动应力作用,从力学机理上深入剖析断裂右旋走滑活动产生的扰动应力对区域应力场(或初始应力场)的影响,并对其附近不同构造部位最大主应力方向偏转现象进行定量解释。初步结果表明：1) 走滑断裂活动导致地应力解耦主要是通过断裂滑动产生的扰动应力对区域应力场的附加影响实现的;2) 在固定的断裂活动长度和摩擦强度情况下,走滑断裂活动导致区域最大主应力方向偏转角度基本与活动断裂走滑位移量呈正比;3) 受汶川MS8.0地震影响,龙门山断裂带大致以茂县-绵竹为界,西南侧的宝兴-映秀-汶川一带最大主应力方向偏转角度很小,基本与区域应力场一致,而东北侧的北川-平武-广元一带则发生角度不等的偏转,3个地区最大主应力方向平均偏转角度分别为65°、56°和29°,该角度与龙门山断裂带东北段同震右旋走滑位移量呈正相关。     

2020-01-23四川石渠MS4.3地震震源机制及发震构造研究

基于数字地震台网波形数据,采用gCAP反演方法解算四川石渠M_S4.3地震的震源机制。结果表明,最佳双力偶解为节面Ⅰ走向134°、倾角82°、滑动角11°,节面Ⅱ走向42°、倾角79°、滑动角171°;最佳质心深度为9 km,矩震级为M_W4.53。为测试震源机制解的稳定性和可靠性,从地壳速度结构、定位误差和数据质量影响等方面进行分析,研究结果表明,这些因素对反演结果影响较小。余震序列展布方向为NW向,主要分布在主震偏北东侧附近,发震断层具有向NE倾的趋势。综合震源区地质构造特征、主震震源机制解和余震序列的空间分布特征,初步判断石渠地震的发震断层面为节面Ⅰ,即长沙贡玛断裂,是一次左旋走滑型地震事件。     

基于地壳介质各向异性分析乳山地震近场区构造应力特征

利用乳山台阵和乳山台记录的资料,采用剪切波分裂方法分析乳山地区地壳各向异性参数,获得该区域应力环境特征。研究结果显示,乳山地区地壳介质各向异性受到区域构造应力和局部构造活动的双重影响,台站DSDC的快剪切波偏振优势方向近NW方向。结合地震序列展布方向及较大地震震源机制推测,该台附近存在近NW方向的活动断裂,该断裂可能为乳山地震序列的发震断层;台站BSTZ和CXZX的快剪切波偏振优势方向近EW方向,附近可能存在近EW方向的活动断裂。     

北天山地区近期重力场及地壳密度变化特征研究

利用2014~2019年北天山地区流动重力资料,采用绝对重力控制进行经典平差计算,获取该地区重力场时空演化特征图像,并结合重力段差可视化方法,得到该地区重点构造部位物质迁移规律。采用紧凑反演方法进行北天山中段不同深度地壳密度模拟计算,结果表明,北天山中段地壳物质密度变化最佳深度为30~40 km,是地壳密度变化的正值区域,基本位于北天山北坡和准噶尔盆地南缘,表明该地区物质不断汇聚,或者表明地壳呈现下沉的趋势;而在北天山山体附近,物质密度变化为负值,表明山体物质处于亏损态势,或者山体受到构造应力的作用,使地壳处于不断隆升过程中。     

寻乌及邻区震源机制和视应力特征

利用寻乌及邻区2009年以来的地震波形资料,采用Snoke方法反演19次ML2.5以上地震的震源机制,计算89次ML1.5以上地震的视应力。结果显示,寻乌及邻区的走滑型地震占绝对优势,震源机制解P轴优势方向分布在310°~320°和110°~130°之间,倾角绝大多数在40°以内,与该区以NWW-SEE向水平挤压为特征的构造应力场吻合;寻乌震群集中区内8次ML2.5以上地震的震源机制结果一致性较好。寻乌及邻区视应力与震级呈正相关,2017-11-06寻乌ML3.8震群前的差视应力值出现高值波动。     

利用CAP方法分析2019-12-05拜城MS4.9地震震源机制解

选取新疆测震台网17个宽频带数字台站记录仪的地震波形数据,采用CAP方法反演2019-12-05拜城MS4.9地震震源机制解。结果显示,节面Ⅰ走向302°、倾角36°、滑动角124°,节面Ⅱ走向83°、倾角61°、滑动角68°,最佳深度解为18 km。同时,采用P波初动方法计算2018-11以来与该事件震中相近的5次MS≥3.0地震,结合以往地质资料分析认为,节面Ⅰ为断层错动面,断层类型为逆冲断层。