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

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

自适应卡尔曼滤波在航空重力异常解算的应用研究

依据航空重力测量基本原理,构建了航空重力异常解算的卡尔曼滤波模型,将新息自适应卡尔曼滤波器(IAE,Innovation based Adaptive Estimation)应用于量测噪声未知的航空重力异常解算.针对IAE滤波器滑动窗口宽度难以准确确定的问题,通过对多个不同滑动窗口新息协方差估计的加权平均,获得改进的IAE滤波器,该IAE滤波器不仅具有量测噪声自适应估计能力,还能实现滑动采样窗口的优化选取.试验结果表明,IAE滤波器可以降低因量测噪声统计信息不明引起的解算误差,改进IAE解算的重力异常误差约为1 mGal.     

GPS应变率场解算方法的对比分析及其可靠性研究

利用GPS速度场资料进行应变率场计算存在多种不同的计算方法,虽然理论上速度场与应变率为简单的偏导关系,但由于GPS观测为点测量方式,如何构建逼近真实的速度场进而计算可靠的应变率场就有其复杂的一面,因此实际计算中不同研究者即使利用相同的数据也很难得到相同的结果,究竟如何评价不同计算方法的可靠性和适用性就显得尤为重要.     

基于多视线向D-InSAR技术的三维同震形变场解算方法研究及应用

差分合成孔径雷达干涉测量技术(differential interferometry synthetic aperture radar,D-InSAR)是地震形变监测的一种重要方法,但是其存在视线向(line of sight,LOS)模糊的问题,即D-InSAR测量的形变场并不是地表真实的垂直向、东西向、南北向形变场,而是地表     

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

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

基于PSInSAR的2008年汶川MS8.0地震震后形变场特征

2008年5月12日汶川M_S8.0地震造成了巨大的同震形变, 地表破裂带分布广泛。 为了研究汶川地震震后形变的分布特征, 收集了震后2008年5月28日至2010年9月15日的多时相ENVISAT ASAR数据, 该数据包括两个相邻条带, 能够覆盖震后形变区域。 使用基于单一主图像的PSInSAR处理技术分别对两个条带进行处理并拼接, 获得了震后时序形变图和年平均形变速率图。 研究结果显示, 龙门山断裂带两侧的震后形变分布特征明显不同, 断裂带西北盘形变量表现为LOS向隆升, 最大累积形变约为26 mm, 最大形变速率约为17±7 mm/a; 断裂带东南盘形变表现为LOS向沉降, 形变量较小, 最大累积形变约为-10 mm, 最大形变速率约为-5 mm/a。 汶川地震后的形变特征反映了龙门山断裂在此段具有逆冲为主兼有走滑分量的活动性质, 本文研究结果与震后GPS观测结果在形变趋势和量级方面具有较好的一致性。     

利用PSInSAR时间序列研究拉奎拉地震位移场变化特征

对2009年4月6日意大利拉奎拉(LAquila)Mw6.3级地震9景ENVISAT重轨单视复数据采用双轨模式进行永久散射体(Persistent Scatterer,PS)处理,获得了此次地震时序位移场,结合PS位移场Delaunay算法数值分析,结果表明:1)Envisat干涉雷达完整而清晰地探测到了地震前后位移场变化过程及其在不同阶段与震源断裂相关的不同形变特征:震前蠕动位移-明显变形-震期快速突变-震后量级明显减缓的持续变形;2)在此次分析的约54×59km2区域内,震中以西为大面积上升区,视线向最大上升量为130mm;震中以东为下降区域,但下降集中在破裂区,主要在发震过程及震后形成,视线向最大下沉量为210mm,与GPS监测结果一致;3)地震产生的破裂主要分布在约22×14km2范围内,沿NW走向、SW倾斜的Paganica-S.Demetrio正断层展布,方向约135°;4)本文为探索利用PSInSAR时间序列进行地震趋势预测提供了一个较好的案例.     

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

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

井间电磁场的一维、二维联合反演方法

提出了针对轴对称二维井间模型的一维、二维联合反演方法.该方法基于纵向成层背景地层的Green函数,以一维反演结果作为联合反演的迭代初始值,同时反演一维背景地层和二维异常剖面的电导率.采用递推技术计算Green函数的系数,可以很方便地同时得到Green函数对空间坐标的积分和对背景层电导率的微分,大大加快了计算雅可比矩阵〖WTHX〗M〖WTBZ〗的速度并使计算方便、准确.反演实例说明了该方法的有效性.反演实例还显示,背景层电导率较二维成像剖面的电导率收敛速度要快,测量数据的精度对二维剖面电导率成像质量的影响大于对背景层电导率分布的影响.     

三维复杂黄土塬区变加密不等距网格初至走时计算与特征分析

三维复杂黄土塬区初至走时计算与特征分析,对于分析复杂黄土塬区地震波运动学特征、观测系统设计、初至拾取及层析反演等均具有重要意义.针对该问题：首先,采用三维分区局部变加密不等距网格剖分复杂黄土塬速度模型;然后,建立基于费马原理和迎风格式的双线性插值局部走时算法,采用修正后的群推进法作为波前扩展方式;最后,以理论模型考察了算法的精度与效率、分析了三维复杂黄土塬区地震初至波的走时特征.结果表明：变加密网格策略能在增加很少计算量的情况下,大幅降低因震源奇异性引起的计算误差问题,同时还能无条件稳定且灵活的适应三维复杂黄土塬区的复杂地震地质条件;和层状均匀模型相比,黄土塬区地震初至波的走时场和时距曲线较为复杂.

     

非弹性反应谱计算的一个改进算法

本文就非弹性地震反应谱的计算问题（运动方程的形式，谱参量的选取，方程解算中恢复力模型上拐折点的处理和时间积分方法），给出了几点理论考虑，介绍了一个改进的算法。这一算法具有精度好和效率高的优点，这对做好结构抗震设计和其它地震工程问题中需要处理大量强震动数据的分析研究工作是有利的。文中还就双线型恢复力模型具体计算了2个例子（一个为屈服强化情形，另一个为理想弹塑性情形）。最后联系计算工作和计算结果，对非弹性反应谱计算中的其它有关事项做了补充讨论和说明。     

2.5维复电阻率反演及其应用试验

自复电阻率法被提出以来,其区分矿与非矿方面的潜力就引起了人们的广泛关注,但反演问题一直没有得到较好的解决,已经严重影响和制约了该方法的应用和发展.对此,本文提出并实现了利用多个排列视电阻率和视相位数据的2.5维SIP联合反演方法.算法利用最小二乘原理构建了反演目标函数,为提高反演的稳定性,在反演方程中加入了Occam法的光滑模型约束.通过借助电场的偏导数形式,推导出了灵敏度矩阵的解析表达式,并应用互换定理对其进行精确求取.反演算法充分利用了电磁感应和激电效应中的异常信息,能够同时反演出二维地质断面上所有单元的四种复电阻率参数.最后,利用该反演程序对安徽某地区的SIP实测数据进行了反演成像,通过与已知钻井资料及CSAMT反演结果的对比分析表明,该反演算法具有良好的应用效果.     

Fine sediment carrying capacity of combined wave and current flows

The so-called fine sediment in many coastal areas and estuaries in China is mostly referred to the mixture of cohesive sediment and non-cohesive sediment. To predict the mixed type time sediment transport, sediment carrying capacity formulae combined with the 2-D suspended sediment transport equation and morphologic equation have been widely used in China. In the present study, the sediment carrying capacity formula suggested by Dou et al. （1995） for wave conditions has been improved and implemented for the prediction of sediment transport in nearshore regions where wave activities are significant. The improvement is based on the wave energy dissipation principle inside and outside the surf zone. In the improved formula, sediment in suspension increases with the magnitude of the wave period and this feature complies with general observations. More than 300 laboratory and field measured data sets have been reviewed and 12 of them have been used to verify and determine the major coefficients in the improved formula. The application of the sediment carrying capacity model in combined wave and current situations shows that the model can faithfully reproduce the cross-shore sediment concentration distributions at the southwest coast of Bohai Bay.     

基于DOA的相干算法研究与应用

波达方向估计(Direction of Arrival estimate)是一种阵列信号处理技术,这种技术普遍应用于雷达和移动通信系统中用来确定信号的到达方向.和其他方法相比,这种方法具有高的角度分辨能力.把DOA引入到地震数据的相干体算法中,通过频率域的聚焦,可以将频率域原本分散的相位聚焦到同一个频率上.对聚焦后的振幅和相位进行多项式拟合,可以得到反映振幅和相位空间变化的地震属性.这些属性具有明确的物理意义,可以应用这些属性识别小尺度断层,检测地质体的轮廓.对实际资料的计算结果表明本方法具有很大的应用前景.