Matlab平台下遥感影像的北方岸线提取研究——以大连长兴岛为例

在Matlab平台上, 利用卫星遥感影像数据, 结合RS与GIS技术, 选取大连长兴岛作为研究区域, 使用阈值分割法、边缘检测法及小波变换法进行了海岸线的提取, 并对结果进行对比分析。试验结果表明, 边缘检测法中Canny算子与小波变换法对海岸线的提取结果与真实地物最为接近, 其中Canny算子可以较为准确地提取围填海区域的海岸线。     

爆炸荷载作用下深埋圆形隧洞的饱和黏弹性土-衬砌系统动力响应

假定土骨架服从标准线性固体黏弹性本构关系,研究了深埋圆形隧洞的饱和黏弹性土-弹性衬砌耦合系统在轴对称爆炸作用下的瞬态动力响应。首先,基于饱和土的Biot模型和衬砌的弹性理论,通过引入势函数和Laplace变换,利用弹性衬砌和饱和黏弹性土界面处的连续性条件以及边界条件,得到饱和黏弹性土体和弹性衬砌位移、应力和孔隙水压力等在Laplace变换域中的解析解。其次,利用Laplace数值Crump逆变换得到耦合系统在时间域的动力响应,数值分析了不同土体模型下土体-衬砌耦合系统的径向位移和环向应力以及土体孔隙水压力等。结果表明：对不同土体模型的土体-衬砌耦合系统,其在爆炸载荷作用下的动力响应性态基本一致,但动力响应的振动周期和幅值等具有明显的差异。同时,对于饱和黏弹性土-弹性衬砌系统,土体黏性参数对土体径向位移和孔隙水压力有明显的影响,但对土体环向应力影响较小。     

关于频率电磁测深几个问题的探讨(八)——频率电磁测深与瞬变电磁测深的关系

从探测原理出发,说明人工源频率电磁测深与瞬变电磁测深的关系密切,这是一个问题的两个方面,可以互相转换。两种探测方法实质是一样的,具有共性;但毕竟从频率域和时间域两个侧面反映问题,又具有各自特性。两种方法的原理所依据的麦克斯韦方程,其正演公式可通过傅立叶变换相互转换;从实际探测所发收的信号特征,场区划分与探测深度等主要方面阐明两种探测方法异同,以加深对探测理论与方法的认识,结合实际择优应用。      

流变边界层及其对华北克拉通热/地震岩石圈底界差异的意义

二维热传导/对流数值模型显示,纯传导的固体岩石圈与纯对流的流体软流圈之间存在一过渡层,即流变边界层,其间传导与对流共同作用来传递热量.流变边界层厚度主要由软流圈黏性系数(η)控制,而受固体岩石圈厚度及热状态影响很小.随着η从1×1021Pa·s降低至1×1019Pa·s,流变边界层也随之减薄,流变边界层的厚度与lg(η)成正比.流变边界层的存在是造成热岩石圈与地震岩石圈厚度差异的重要因素.全球典型克拉通岩石圈的对比结果表明,地震岩石圈厚度普遍大于热岩石圈厚度,二者的差异多数在70~90km,很好地验证了流变边界层的存在.研究发现二者的差异在华北克拉通自西向东逐渐减小:由西部鄂尔多斯的约80km减少至渤海湾盆地的约20km.反映出华北克拉通岩石圈下部流变边界层厚度自西向东减薄,意味着软流圈黏性系数自西向东逐渐降低.这可能与中生代太平洋俯冲脱水形成的低黏大地幔楔有关,从一侧面印证了太平洋俯冲对华北克拉通破坏的影响.     

基于压缩感知的Curvelet域联合迭代地震数据重建

由于野外采集环境的限制,常常无法采集得到完整规则的野外地震数据,为了后续地震处理、解释工作的顺利进行,地震数据重建工作被广泛的研究.自压缩感知理论的提出,相继出现了基于该理论的多种迭代阈值方法,如CRSI方法（Curvelet Recovery by Sparsity-promoting Inversion method）、Bregman迭代阈值算法（the linearized Bregman method）等.CSRI方法利用地震波形在Curvelet的稀疏特性,通过一种基于最速下降的迭代算法在Curvelet变换域恢复出高信噪比地震数据,该迭代算法稳定,收敛,但其收敛速度慢.Bregman迭代阈值法与CRSI最大区别在于每次迭代时把上一次恢复结果中的阈值前所有能量都保留到本次恢复结果中,从而加快了收敛速度,但随着迭代的进行重构数据中噪声干扰越来越严重,导致最终恢复出的数据信噪比低.综合两种经典方法的优缺点,本文构造了一种新的联合迭代算法框架,在每次迭代中将CRSI和Bregman的恢复量加权并同时加回本次迭代结果中,从而加快了迭代初期的收敛速度,又避免了迭代后期噪声干扰的影响.合成数据和实际数据试算结果表明,我们提出的新方法不仅迭代快速收敛稳定,且能得到高信噪比的重建结果.     

InSAR约束下的2008年汶川地震同震和震后形变分析

2008年5月12日,青藏高原东缘的龙门山断裂带上发生了M_w7.9级汶川地震.本文通过分析覆盖汶川地震震中区域的ALOS/PALSAR像对的方位向偏移量来选择无明显电离层扰动影响的像对进行干涉处理,获取了高精度、连续的InSAR地表形变场.在此基础上,结合高精度GPS同震形变数据,采用同震、黏弹性松弛震后形变联合反演模型同时确定了汶川地震的同震滑动分布和龙门山地区的流变结构参数.研究结果表明,汶川地震是一个断层破裂非常复杂的地震事件,其中,北川段、岳家山段、虹口段和汉旺段的滑动以逆冲为主,而青川段以右旋走滑为主.滑动主要发生在10 km深度以上的区域,最大滑动量位于虹口段的东北端,达10.7 m.地震释放的总能量为9.28×10²⁰ N·m（M_w7.91）,与地震学的结果一致.联合反演模型确定的龙门山地区中下地壳的黏性系数为2×10¹⁸ Pa·s,为青藏高原东部地区的黏性系数提供了一个可靠的下限值.如果有更长时间的震后形变观测时间序列,将为该区域提供更为可靠的流变结构.     

流变边界层及其对华北克拉通热/地震岩石圈底界差异的意义

二维热传导/对流数值模型显示,纯传导的固体岩石圈与纯对流的流体软流圈之间存在一过渡层,即流变边界层,其间传导与对流共同作用来传递热量.流变边界层厚度主要由软流圈黏性系数（η）控制,而受固体岩石圈厚度及热状态影响很小.随着η从1×10²¹Pa·s降低至1×10¹⁹ Pa·s,流变边界层也随之减薄,流变边界层的厚度与lg（η）成正比. 流变边界层的存在是造成热岩石圈与地震岩石圈厚度差异的重要因素. 全球典型克拉通岩石圈的对比结果表明,地震岩石圈厚度普遍大于热岩石圈厚度,二者的差异多数在70~90 km,很好地验证了流变边界层的存在. 研究发现二者的差异在华北克拉通自西向东逐渐减小：由西部鄂尔多斯的约80 km减少至渤海湾盆地的约20 km. 反映出华北克拉通岩石圈下部流变边界层厚度自西向东减薄,意味着软流圈黏性系数自西向东逐渐降低.这可能与中生代太平洋俯冲脱水形成的低黏大地幔楔有关,从一侧面印证了太平洋俯冲对华北克拉通破坏的影响.     

太阳风压力系数的研究

利用全球磁流体力学(MHD)的模拟结果,研究了太阳风压力系数与上游太阳风参数和日下点磁层顶张角的相关性.在识别出日下点附近磁层顶位置后,通过拟合得到日下点附近的磁层顶张角.在考虑上游太阳风中的磁压和热压以及磁层顶外侧的太阳风动压的情况下,计算了太阳风压力系数.通过分析行星际磁场不同方向时太阳风动压在日地连线上与磁压和热压的转化关系,详细研究了太阳风参数和日下点磁层顶张角对太阳风压力系数的影响,得到以下相关结论:(1) 在北向行星际磁场较大(B_z≥5 nT)时,磁层顶外侧磁压占主导,南向行星际磁场时磁层顶外侧热压占主导;(2) 太阳风压力系数随着行星际磁场的增大而增大,随着行星际磁场时钟角的增大而减小;并且在行星际磁场大小和其他太阳风条件相同时,北向行星际磁场时的太阳风压力系数要大于南向行星际磁场时的;北向行星际磁场时,太阳风压力系数随着太阳风动压的增大而减小,南向行星际磁场时,太阳风压力系数随着太阳风动压的增大而增大;以上结论是对观测结果的扩展;(3) 最后,我们还发现太阳风压力系数随着日下点磁层顶张角的增大而增大.     

海原断裂带中东段地貌差异及其成因探讨

以定量化地形因子为切入点的构造地貌学方法已成为活动构造研究的有效手段,被广泛用于定性或半定量解析地貌对新构造运动的响应及其演化过程。针对海原断裂带中东段现今地貌差异,以SRTM 90m分辨率DEM为基础,利用ArcGIS软件和Matlab程序脚本,提取了海原断裂带中东段高程、坡度、地形起伏、地形侵蚀以及河流陡峭系数等地形因子。从空间分布上看,上述各项地形因子沿断裂走向均呈现 “西高东低”的整体分布特征。西段海拔高、坡度陡、起伏大、侵蚀强、抬升快,中段和东段海拔低、坡度缓、起伏小、侵蚀弱、抬升慢,此外,在断裂带的东南尾端呈略微增加趋势,达到小范围内的峰值。在此基础上,通过对比分析地形因子与年降水量、基岩岩性,初步探讨了构造与降水、岩性等因素对地形地貌的控制作用,认为不同降水条件对地貌后期改造起显著作用,基岩岩性与现今地貌之间并无显著关系,该区域地貌类型主要受构造抬升差异所控制。沿断裂带走向上的现今地貌差异表明,西段处于相对快速的构造隆升和强挤压造山构造背景,中段由于受到黄河下切及河流冲积作用影响,地貌参数记录的抬升特征并不显著,而东段则反映出大型断裂带尾端挤压调整效应。     

高采样率四分量钻孔应变同震响应分析

In this paper, we adopt the observation data of two 10Hz YRY 4-component borehole strainmeters and one 100Hz FBS-3B broadband seismometer to study the responses of their observations to the April 11, 2012 Mw 8.6 earthquake off the west coast of northern Sumatra, and compare the response characteristics during the coseismic stage after preprocessing of the data. We analyze the spectrum＇s dynamic process during the coseismic stage with the S transformation method, and discuss the spectrum＇s details in different seismic phases. The S transformation results of the data of seismometer and borehole strainmeter at the same station show that the coseismic spectrum evolution of the two kinds of observation are similar to each other and can be cross-referenced in indentifying the wave signals, which proves that the borehole strainmeter measurements at high frequencies are reliable.