基于多卫星平台永久散射体雷达干涉提取三维地表形变速度场

永久散射体雷达干涉(PSI)技术及其应用于区域地表形变监测已成为雷达遥感领域的研究热点之一.使用单一卫星平台所获取的单侧视SAR影像时间序列进行PSI分析,仅能获取沿雷达视线(LOS)方向的一维地表位移信息.本文提出了基于多平台永久散射体雷达干涉提取三维地表形变速度场的模型与算法,其基本策略是:首先针对每一卫星平台的SAR影像时间序列进行PSI分析,并计算各地面目标沿LOS向的位移速度值,然后联合各平台所对应的LOS向位移速度值进行建模,并基于最小二乘方法解算各地面目标的三维位移速度分量.实验选取天津市西北部作为测试区,使用2007—2010年所获取的39幅TerraSAR-X影像、23幅ENVISAT ASAR影像和16幅ALOS PALSAR影像进行分析,经联合解算得到了该测试区域的垂直位移速度场以及南北向和东西向水平位移速度分量.与地面水准和已有GPS观测结果对比分析表明:多平台PSI的垂直位移速度场精度可达毫米级,而其水平位移速度分量与已有GPS结果基本一致.多平台PSI分析无需引入任何外部形变参考信息,便可以实现形变场的偏差校准和三维形变场的恢复.     

基于多卫星平台永久散射体雷达干涉提取三维地表形变速度场

永久散射体雷达干涉(PSI)技术及其应用于区域地表形变监测已成为雷达遥感领域的研究热点之一.使用单一卫星平台所获取的单侧视SAR影像时间序列进行PSI分析,仅能获取沿雷达视线(LOS)方向的一维地表位移信息.本文提出了基于多平台永久散射体雷达干涉提取三维地表形变速度场的模型与算法,其基本策略是:首先针对每一卫星平台的SAR影像时间序列进行PSI分析,并计算各地面目标沿LOS向的位移速度值,然后联合各平台所对应的LOS向位移速度值进行建模,并基于最小二乘方法解算各地面目标的三维位移速度分量.实验选取天津市西北部作为测试区,使用2007—2010年所获取的39幅TerraSAR-X影像、23幅ENVISAT ASAR影像和16幅ALOS PALSAR影像进行分析,经联合解算得到了该测试区域的垂直位移速度场以及南北向和东西向水平位移速度分量.与地面水准和已有GPS观测结果对比分析表明:多平台PSI的垂直位移速度场精度可达毫米级,而其水平位移速度分量与已有GPS结果基本一致.多平台PSI分析无需引入任何外部形变参考信息,便可以实现形变场的偏差校准和三维形变场的恢复.     

雷达干涉PS网络的基线识别与解算方法

时序雷达干涉图中的永久散射体(PS)可看作“天然GPS点”, 以构成网络用于监测长期的地表形变. 本文提出采用邻接矩阵拓扑模型对基于Delaunay剖分算法生成的PS网络进行基线识别, 并采用时序相干最大化算法求解PS基线的线性形变速度增量和高程误差增量. 该数据模型和计算方法被应用于探测香港地区2006~2007年间的区域地表沉降. 实验研究采用由Envisat卫星ASAR传感器对该地区成像所获取的时序SAR影像作为数据源, 并联合该地区12个GPS连续运行参考站的观测数据予以大气修正和地面控制. 实验结果表明, 该模型和方法应用于地表形变测量是有效的和可靠的, PS网络方法探测地面沉降的精度约为±2.0 mm/a.     

雷达差分干涉测量技术及其应用研究--以1997年西藏玛尼7.9级地震区域形变测量为例

雷达差分干涉测量是一种最新的大地形变测量遥感技术方法,本文对差分干涉技术的原理及具体实现进行了深入研究,并讨论了差分干涉测量结果的误差.利用欧洲空间局1996年4月15日获取的ERS-1数据和4月16日获取的ERS-2数据,以及1997年12月2日获取的ERS-2数据,应用差分干涉测量技术对发生于1997年11月8日的西藏玛尼地震进行了提取区域形变场的应用研究,从得到的变化检测条纹图中可识别出地表破裂带,还可定量推算震中周围和两条断裂带附近的变形情况,差分干涉测量结果与地面调查符合得很好.     

雷达差分干涉测量技术及其应用研究——以1997年西藏玛尼7.9级地震区域形变测量为例

雷达差分干涉测量是一种最新的大地形变测量遥感技术方法，本文对差分干涉技术的原理及具体实现进行了深入研究，并讨论了差分干涉测量结果的误差。利用欧洲空间局1996年4月15日获取的ERS-1数据和4月16日获取的ERS-2数据，以及1997年12月2日获取的ERS-2数据，应用差分干涉测量技术对发生于1997年11月8日的西藏玛尼地震进行了提取区域形变场的应用研究，从得到的变化检测条纹图中可识别出地表破裂带，还可定量推算震中周围和两条断裂带附近的变形情况，差分干涉测量结果与地面调查符合得很好。     

CRInSAR与PSInSAR联合探测区域线性沉降研究

本文提出一种CRInSAR与PSInSAR联合解算算法,并将该算法应用于区域线性沉降探测中.该算法将CR点上计算得出的形变速率值及高程改正值作为研究区域PS基线网络的约束,进而通过间接观测平差法估计出PS网沉降速率和高程改正值的全局最优解.算法实现了人工角反射器雷达干涉测量技术(CRInSAR)与永久散射体雷达干涉测量...     

基于PSInSAR技术的海原断裂带地壳形变初步研究

常规差分干涉测量(DInSAR)受时间、空间失相干的严重制约和和大气延迟等相位误差的影响,难以实现对长期累积微小地壳形变场的有效探测.PSInSAR技术克服了常规DInSAR的局限性,能够高精度监测微小地壳形变.本文首先介绍了PSInSAR技术的算法模型和处理方法.该方法通过二维线性相位模型,对时序干涉图象上相干点目标的差分干涉相位进行回归分析,逐次消除大气延迟、轨道残余和地形残余等相位误差,提取出准确可靠的形变相位,进而得到相干点目标上的累积形变量和形变速率.在此基础上,以祁连山海原断裂带为实验研究区,利用2003~2009年的21景ENVISAT ASAR数据,采用上述相干点目标处理方法进行了海原断裂带地壳微小形变的探索性研究,得到海原断裂东段绝大多数高相干点的形变速率在6~7 mm/a,2003~2009年共6.3年累积位移平均值约为4.2 cm,运动性质为左旋走滑.这一结果与GPS,库仑应力反演和地质学方法得到的结果基本一致,说明PSInSAR技术在长期累积微小地壳形变探测中具有广阔的应用前景和巨大发展潜力,有望成为探测震间微小形变,获取地震形变异常的有效途径之一.     

基于相邻虚拟道叠加的超虚折射干涉法及其在广角OBS折射波增强中的应用

广角地震的远偏移距折射初至含有地下深层的信息,提高远偏移距折射波信噪比能够有效提高初至拾取精度,对于深部结构的层析速度建模十分有利.超虚折射干涉法基于干涉原理对远道折射波进行增强,它在炮点和检波点都非常密集的情况下效果较好.然而,在应用于台站间距较大且环境噪声较强的广角海底地震仪（OBS）观测时,该方法对折射波的增强能力不足,而且容易产生虚假波形,造成增强后的折射波信噪比仍然较低.针对这种情况,本文提出基于相邻虚拟道叠加的超虚折射干涉法,通过叠加相邻虚拟道来提高远道与近道互相关的准确度,以达到稳定增强远道折射波信噪比的目的.理论实验和实际资料测试均显示,基于相邻虚拟道叠加的超虚折射干涉法在台站间距较大和信噪比较低的情况下能够准确构建虚拟道,且增强后的波形同相轴连续程度和信噪比均较高,有利于拾取高精度的折射初至到时.本文方法也可用于增强陆上炮间距较大的广角地震数据折射波.

     

探地雷达在海塘堤脚淘空损伤检测中的应用试验及分析

本文在简要介绍探地雷达方法原理的基础上,详细分析了浙江某海塘堤脚沉井防冲结构及其相关介质的物理性质和空间分布特点,从而指出了堤脚淘空损伤雷达检测的基本可行性及复杂性.最后通过笔者实际参与的雷达检测试验和综合分析表明,采用雷达检测可指出淘空损伤现象存在的可能性及大致位置和范围,而难以断定损伤的性质及尺寸.     

用接收函数建立区域模型的震源机制反演及其在芦山地震序列研究中的应用

本文提出并试验了一种基于接收函数建立区域模型进行震源机制反演的方法.选取四川地震台网记录的M≥3且信噪比高的近震波形资料,反演得到了芦山地震序列中74个地震的震源机制.通过对震源深度和震源机制的综合分析,探讨了芦山地震的发震构造和区域应力场状态.采用接收函数方法反演获取了26个台站下方的S波速度结构,对不同区域的台站反演结果进行叠加平均,以此区域平均S波速度作为本文震源机制反演使用的区域模型的S波速度;区域模型的P波速度由经验公式给出.反演稳定性测试表明,使用不同模型或对原始波形记录加入随机噪声的反演结果与原始反演相比,震源深度最大误差为1 km,断层面各参数误差水平也很低,且显示的发震类型是一致的,其中随机噪声带来的误差小于模型带来的误差.主震反演得到的震源机制解为:震源深度17 km,矩震级6.47;节面Ⅰ走向213°,倾角51°,滑动角98°;节面Ⅱ走向20°,倾角40°,滑动角80°;显示芦山主震可视为纯逆冲型地震,发震构造可能是某个具有较大倾角的逆冲断层,而不是低缓的推覆构造的基底滑脱面.同时本文反演获取的73个M≥3余震的震源机制绝大多数也显示了类似的发震类型,逆冲型地震为67个,占92%,具有绝对优势;走滑型地震为5个,正断型地震为1个.其中5个走滑型地震中的4个均分布在震源区的东北端.整个芦山地震序列深度集中在12~20 km,且沿震源区短轴的余震深度剖面有自西向东呈逐步变浅的趋势,呈现清晰的铲形断面结构,结合本地地质构造,可以推断芦山地震序列主要发生在龙门山前山断裂以东的逆冲推覆体内的一个隐伏断裂上.P轴方位角优势方位与区域应力场及汶川震源区南段的相一致,表明芦山序列地震活动主要受区域应力场控制,且汶川震后该区应该不存在应力场变化.P轴仰角随深度分布则显示了孕震层在浅部为脆性上地壳,而深部已经进入了中地壳低速层.断层面的几何形态简单,倾角均值在不同深度保持稳定在55°左右,与主震倾角接近,这与汶川震源区南段的研究结果明显不同,揭示了龙门山断裂带南段与此次芦山发震断裂在断层面几何形态上的明显差异.     

An automated cross‐correlation based event detection technique and its application to a surface passive data set

In studies on heavy oil, shale reservoirs, tight gas and enhanced geothermal systems, the use of surface passive seismic data to monitor induced microseismicity due to the fluid flow in the subsurface is becoming more common. However, in most studies passive seismic records contain days and months of data and manually analysing the data can be expensive and inaccurate. Moreover, in the presence of noise, detecting the arrival of weak microseismic events becomes challenging. Hence, the use of an automated, accurate and computationally fast technique for event detection in passive seismic data is essential. The conventional automatic event identification algorithm computes a running‐window energy ratio of the short‐term average to the long‐term average of the passive seismic data for each trace. We show that for the common case of a low signal‐to‐noise ratio in surface passive records, the conventional method is not sufficiently effective at event identification. Here, we extend the conventional algorithm by introducing a technique that is based on the cross‐correlation of the energy ratios computed by the conventional method. With our technique we can measure the similarities amongst the computed energy ratios at different traces. Our approach is successful at improving the detectability of events with a low signal‐to‐noise ratio that are not detectable with the conventional algorithm. Also, our algorithm has the advantage to identify if an event is common to all stations (a regional event) or to a limited number of stations (a local event). We provide examples of applying our technique to synthetic data and a field surface passive data set recorded at a geothermal site.     

高分辨率非线性储层物性参数反演方法和应用

对于陆相沉积环境下的复杂隐蔽岩性储层,由于观测信息不准确,如信息重叠、信息缺失和噪音污染,以及岩石物理关系模糊等原因,储层横向预测存在不惟一性、不稳定性和不确定性.基于线性假定的常规储层横向预测技术已不适用于复杂隐蔽岩性储层的勘探.本文采用一种非线性储层岩性物性褶积模型,建立波阻抗与孔隙度/泥质含量的函数关系;通过多级结构分解和双向边沿子波检测来刻画复杂岩石物理关系;通过Caianiello褶积神经网络实现确定性反演、统计反演和非线性理论三者有机结合;最后联合应用基于逆算子的反演方法和基于正算子的重建算法实现了综合地质、测井和地震波阻抗信息进行高分辨率储层物性参数反演.非线性储层物性参数反演采用多井约束机制和分频反演方式,在陆相和近海油气勘探资料的实际应用中,取得了明显应用效果.     

基于ANSYS的灵敏度分析及其在单层网壳中的应用

基于实验数据的灵敏度分析是一种可靠有效的分析方法,详细阐述了ANSYS中基于实验数据的灵敏度分析基本原理,为基于ANSYS的结构灵敏度分析提供参考;应用ANSYS对一个Schwedler型单层球面网壳进行了灵敏度分析,得出了几种结构响应对于杆件的灵敏度值,将不同位置的杆件对于几种结构响应的影响程度进行了量化,可为该类型单层网壳的设计及响应控制提供参考,通过对计算结果的分析验证了基于实验数据灵敏度分析的高效性和可靠性.     

Petrophysical evaluation and its application to AVO based on conventional and CMR-MDT logs

常规测井为 AVO 分析提供了基础的资料, 成为联系岩石物理与地震资料的桥梁。然而如果储层存在有复杂的流体系统, 如地层被严重地层流体侵入、电阻率响应低及盐水矿化度复杂等的现象, 则常规测井无法提供高质量的测井资料, 导致得出错误的弹性计算结果, 使 AVO 结果与地震资料不吻合。中国渤海湾地区第三系裂缝性储层复杂, 我们利用常规测井和核磁共振测井与模块地层动态测试相结合的组合仪完成了地层评价和储层描述。研究结果表明岩石物理学家利用上述方法技术可以获得诸如空隙度、渗透率、含水饱和度、束缚流体以及空隙压力等重要的储层参数并进一步综合应用这些结果和以实验室测量数据为基础的岩性分析结果进行在地震域岩石物理研究和 AVO 分析。