基于分区解缠方式研究2010年玉树M_W6.9地震同震形变

2010年玉树MW6.9地震发生后,已有学者利用DInSAR技术得到了该次地震的同震形变场,并在此基础上反演了其震源破裂滑动分布。本文以提高玉树地震同震形变场准确度作为出发点,设计了一种新的解缠方式获取相应的同震形变场:首先沿地表同震破裂迹线将InSAR干涉图像分割成上下两部分,然后利用网络流解缠算法对其分别进行相位解缠,最后通过统计重叠部分的相位一致性信息将两部分结果进行拼接。利用震中区GPS同震位移数据对分区解缠结果进行验证。结果表明:相对于整体解缠结果,分区解缠方式得到的最大视线向沉降值由42.6cm增加到48.1cm,所得结果与GPS观测结果更加接近,采用分区解缠方式提高了断层附近形变场的准确性。     

相位空间内的解缠绕相位反演

全波形反演方法是一种数据域高精度反演方法, 该方法通过匹配观测数据与模拟数据的地震波形, 利用梯度法准确反演地下介质参数的分布情况.由于观测数据普遍缺少低频信息, 该方法易受周期跳跃现象影响.特别是当地下存在大尺度强反射界面的构造时, 地下介质的反演转化为强非线性问题求解.该情形下, 即使观测数据包含充足的低频信息, 全波形反演也难以给出准确的反演结果.一般可以通过减弱反演对初始模型参数的依赖性来克服上述问题, 具体表现为使用新变量（例如瞬时相位、包络等）代替目标函数中的采样后波场, 以增强新目标函数的凸性.但是, 对该新目标函数进行反演时, 伴随状态方程中存在关于新变量和波场的一个链式微分项, 该项保留了反演问题的非线性, 导致新的反演方法难以处理包含大尺度构造的强非线性反演问题.此外, 基于新变量的反演问题依然在波场空间中计算模型梯度, 难以充分利用新变量与模型参数之间的弱非线性关系.因此, 本文提出用频率域波动方程的相位形式代替传统的波动方程来消除伴随状态方程中的链式微分项, 用解缠绕的相位代替目标函数中采样前波场并在相位空间进行反演.该方法可以最大程度地利用地下介质参数和解缠绕相位之间的弱非线性关系, 从而削弱反演的非线性性.由于基于频率域波场计算得到相位有严重的缠绕问题, 本文采用基于振幅排序的多聚类算法来对相位进行解缠绕.虽然将介质参数到波场的映射替换为介质参数与解缠绕相位的映射, 会导致反演结果的分辨率有所下降, 但该方法可以在相位空间恢复介质参数的大尺度低波数分量.Marmousi模型测试证明了该方法的有效性和准确性, 针对部分BP模型的测试也证明了该方法处理强非线性问题的能力.

     

形变观测中小数据集PSInSAR的误差分析与初步应用

PS InSAR技术在微小形变的监测中有着独特的优势, 但同时受到数据量的严重制约。 本文以形变速率稳定的怀来县城为实验区, 利用15景ASAR数据, 依据主影像选择和高程改正选取了三种略有差异的模型分别进行PS处理, 并对形变提取中的相位解缠、 基线改正和高程改正进行误差分析。 结果表明, 小数据集的相位解缠和高程改正在部分PS点上容易发生错误, 但形变估计整体上能够正确反映实际情况。 三个模型中MS(Multi-Reference Single Regression) 模型得出结果最为可靠。 小数据集解算的整体PS点质量不高, 可以在满足PS点密度要求的前提下, 降低相位标准偏差阈值以提高解算的可靠性。     

基于PSInSAR技术的长白山天池火山形变监测

将PSInSAR技术引入长白山天池火山的形变监测,获取了1992—1998年和2007—2010年2个时段的火山形变信息。结果显示:天池火山在这2个时段内整体抬升,1992—1998年火山较为活跃,雷达视线向平均形变速率为6mm/a,2007—2010年火山活动趋于平缓,雷达视线向平均形变速率为3mm/a;结合水准和GPS数据分析,发现火山口区域地表抬升明显,远离火山口处较为稳定。文中PSInSAR结果与水准数据能较好地吻合,且在空间上有较大覆盖范围,能更直观地反映火山地表的形变特征。     

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

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

基于广义S变换的地震资料振幅谱补偿和相位谱校正方法研究（英文）

由于地下介质对地震波振幅的影响和地震波频散因素,地震波振幅和相位随时间、空间及频率的变化而发生改变,本文提出一种基于广义S变换的振幅谱补偿和相位谱校正新方法。该方法在S域中分为振幅谱补偿和相位谱校正两个步骤进行处理:振幅谱补偿是在地震记录可靠频带范围内恢复反射系数的振幅谱,其具体实现是在S域中利用谱模拟技术来拟合时变子波振幅谱,从而补偿由地层吸收所引起的振幅衰减;相位谱校正是消除子波剩余相位的影响,其具体实现是在S域中利用相位扫描来拾取随时间、空间和频率而变化的相位校正量,并由Parsimony准则来进行最佳相位判别。本文方法不需要直接求取Q值,能够适用于变Q值情况。理论模型和实际资料处理表明,该方法不仅能恢复地层反射系数的振幅谱,还可以有效消除子波剩余相位的影响,使子波接近或达到零相位,从而提高地震资料分辨率。     

改进的阻尼法低纬度磁异常化极方法

中低纬度地区的磁异常资料进行化极处理能减少或消除斜磁化影响,对提高磁测资料的利用程度有重要意义.本文在总结分析原有的低纬度磁异常化极过程的基础上,改进了原有化极方法的滤波过程.改进后的阻尼化极方法,采用将振幅谱和相位谱分开进行滤波的过程,过程中只需要考虑对化极因子的振幅谱进行压制即可.通过公式分析表明该方法更加合理,减少了人为因素的影响;经单个模型和组合模型试验表明,此化极方法准确、稳定、可靠,可以有效地突出低纬度地区的异常体信息.

     

求解线性方程组的一种迭代方法——ART算法的一种改进

本文用Ｓｃｈｍｉｄｔ正交化方法对阻尼ＡＲＴ算法进行了改进，并以此为基础，通过与一种投影迭代法的结合给出一种求解线性方程组的迭代方法，同时指出即使线方程组不相容时，由此法所构造的向量序列将收敛于其Ｍｏｏｒｅ－Ｐｅｎｒｏｓｅ广义逆解。     

强地震动场相位谱的一种工程预测方法

本文在强震地震波频散研究的基础上，进行了强地震动场相位预测方法和多点地面运动时程合成方法的研究，作者首先运用ＳＭＡＲＴ－１台阵资料对强地震波等效频散模型进行了深入研究，发现表示地震动相位随机性的随机过程η（ｆ）的相位服从正态分布，并用此分布取代了原模型中的均匀分布的假定，然后将这一改进的等效频散模型推广应用于强地震动场相位的预测论文中提出的方法既考虑了相邻点相位的行波相关性，亦考虑了多种复杂因素产     

基于保角变换的双谱域相位谱估计方法及其在地震子波估计中的应用(英文)

地震子波估计是地震资料处理与解释中的重要环节,它的准确与否直接关系到反褶积及反演等结果的好坏。高阶谱(双谱和三谱)地震子波估计方法是一类重要的、新兴的子波估计方法,然而基于高阶谱的地震子波估计往往因为高阶相位谱卷绕的原因,导致子波相位谱求解产生偏差,进而影响了混合相位子波估计的效果。针对这一问题,本文在双谱域提出了一种基于保角变换的相位谱求解方法。通过缩小傅里叶相位谱的取值范围,有效避免了双谱相位发生卷绕的情况,从而消除了原相位谱估计中双谱相位卷绕的影响。该方法与最小二乘法相位谱估计相结合,构成了基于保角变换的最小二乘地震子波相位谱估计方法,并与最小二乘地震子波振幅谱估计方法一起,应用到了地震资料混合相位子波估计中。理论模型和实际资料验证了该方法的有效性。同时本文将双谱域地震子波相位谱估计中保角变换的思想推广到三谱域地震子波相位谱估计中。     

An improved implicit-explicit time integration method for structural dynamics

An explicit predictor-corrector algorithm is derived from the implicit α-method. This explicit algorithm is shown to have better stability and accuracy properties than its Newmark-based predecessor. This algorithm is then combined with the implicit α-method, resulting in an implicit-explicit α-method which can be effectively utilized for linear and non-linear structural dynamics calculations.     

改良的数字化测图新方法

阐述了在地形复杂、地貌破碎、作业组较多时一种改良的数字化测图新方法的作业流程。该方法发挥了数字化测图的高精度高效率和模拟法测图善于处理勾绘复杂地貌等高线的长处，为针对不同情况制定合适的测图作业方案增加了一种选择。     

An improved capacity spectrum method for ATC‐40

In order to account for the non‐linear behavior of structures via non‐linear static procedure, the capacity spectrum method has been adopted by ATC‐40 for evaluation and retrofit of reinforced concrete buildings. For elastic‐perfectly‐plastic SDOF systems, the accuracy of the capacity spectrum method depends only on the acceleration response spectrum chosen to form the demand spectrum and the adopted model for calculating the equivalent viscous damping ratios. According to this method, the pseudo‐acceleration response spectrum (PS_a) is used to create the demand diagram. It is found that the ATC‐40 procedure, using its Type A hysteretic model, may be inaccurate especially for systems with damping ratios greater than 10% and periods longer than 0.15sec. In order to improve the accuracy of the capacity spectrum method, this study proposes to use the real absolute acceleration response spectrum (S0._a) instead of the PS_a to establish the demand diagram. The step‐by‐step procedure of the improved method and examples are implemented in this paper to illustrate the calculations of earthquake‐induced deformations. In addition, three selected models of equivalent viscous damping are also compared in this paper to assess the accuracy of the model used in the ATC‐40 procedure. Results show that the WJE damping model may be used by the capacity spectrum method to reasonably predict the inelastic displacements when the ductility demand (μ) of the structures is less than 4, whereas the damping model proposed by Kowalsky can be implemented when μ>4.0. Alternatively, the damping model proposed by Kowalsky may be used to calculate the equivalent viscous damping for the entire range of ductility. Copyright © 2003 John Wiley & Sons, Ltd.     

An improved oxygen method for measuring primary production in lakes

The accuracy of the Winkler technique for measuring dissolved oxygen was increased to 0.002–0.005 mg O₂/l with a photometric endpoint detection method. This precision is valuable for accurate measurements of primary production only if an oxygen release on illuminating dark adapted samples and the instability of the oxygen concentration of undersaturated water are taken into account. This work was supported by the ‘Swiss National Science Foundation’.     

An improved optimal elemental method for updating finite element models

The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results. This situation occurs when the test modal model is incomplete, as is often the case in practice. An improved optimal elemental method is presented that defines a new objective function, and as a byproduct, circumvents the need for mass normalized modal shapes, which are also not readily available in practice. To solve the group of nonlinear equations created by the improved optimal method, the Lagrange multiplier method and Matlab function fmincon are employed. To deal with actual complex structures, the float-encoding genetic algorithm (FGA) is introduced to enhance the capability of the improved method. Two examples, a 7degree of freedom (DOF) mass-spring system and a 53-DOF planar frame, respectively, are updated using the improved method.The example results demonstrate the advantages of the improved method over existing optimal methods, and show that the genetic algorithm is an effective way to update the models used for actual complex structures.     

An improved method of aggregate stability measurement

The stability of individual soil aggregates to impacting water drops is now regarded as a fundamental parameter required in any soil erodibility assessment. This stability has been measured by the use of single drop rainfall simulators, the characteristics and operation of which have, in the past, shown considerable variation. The drop sizes, impact frequency, drop fall height, chemistry of drop-forming liquid, aggregate pretreatment, and definition of breakdown affect the way in which the apparatus can be used to define stability. Within this paper each one of these characteristics is considered and a piece of apparatus proposed which is believed to give results with increased consistency and replicability. The described apparatus can be readily reproduced in any laboratory situation, and has the added advantage of being not only appropriate for the measurement of indices of stability, but also of providing an environment in which the mechanisms of aggregate breakdown can be investigated.     

An improved pseudo-static method for seismic resistant design of underground structures

This paper describes a commonly used pseudo-static method in seismic resistant design of the cross section of underground structures. Based on dynamic theory and the vibration characteristics of underground structures, the sources of errors when using this method are analyzed. The traditional seismic motion loading approach is replaced by a method in which a one-dimensional soil layer response stress is differentiated and then converted into seismic live loads. To validate the improved method, a comparison of analytical results is conducted for internal forces under earthquake shaking of a typical shallow embedded box-shaped subway station structure using four methods： the response displacement method, finite element response acceleration method, the finite element dynamic analysis method and the improved pseudo-static calculation method. It is shown that the improved finite element pseudo-static method proposed in this paper provides an effective tool for the seismic design of underground structures. The evaluation yields results close to those obtained by the finite element dynamic analysis method, and shows that the improved finite element pseudo-static method provides a higher degree of precision.     

An improved straight-line fitting method for analyzing pumping test recovery data

Theis (1935) derived an exact solution for the residual drawdown in a well after the cessation of a pumping test by summing two drawdowns: one (s1), caused by imaginary continuation of the original pumping and the other (s2), due to an imaginary injection at the same constant rate. We approximated the Theis solution to obtain a simple linear relation for determining the transmissivity and storage coefficient from recovery data. Unlike other existing straight-line fitting methods, in our method, we applied different approximations to the well functions in the solutions of s1 and s2. We used the well-known Cooper-Jacob approximation for s1, truncating the expansion of the well function in s2 to its first three terms. For the same level of truncation errors, while the Cooper-Jacob approximation requires the argument u1相似文献   

改进的基础隔震结构地震作用简化计算方法

在《建筑抗震设计规范》(GB50011-2001)关于隔震结构的简化计算方法中,水平向减震系数的表达式和定义有些不符,假定的隔震结构地震作用分布规律也与实际情况略有出入。本文基于水平向减震系数的定义和实际隔震结构的地震作用分布规律提出了一种改进的隔震结构水平向减震系数、隔震结构总地震作用、隔震结构上部地震作用分布的计算方法,并提出了总水平地震作用减震系数的新概念。本文提出的改进方法具有表达准确、物理意义明确的特点。将本文提出的改进算法计算结果与时程分析计算结果比较,结果显示,改进方法的计算结果与时程分析结果接近,且分布规律一致。     

一种改进的基于网格搜索的微地震震源定位方法

震源定位是微地震监测技术要解决的主要问题.目前,井下微地震监测多采用走时拟合法计算震源位置.常规方法受到环境噪声、初至拾取误差、速度模型误差等因素的影响,定位结果存在一定误差.为了提高定位精度,本文提出了一种改进的基于网格搜索的微地震震源定位方法.本文方法根据P波的偏振特征参数计算概率密度函数求取震源方位角,并采用改进的目标函数和搜索算法计算震源的径向距离和深度.模型数据和实际资料的处理结果表明,本文方法具有较强的抗噪性,计算得到的震源方位角更加接近真实值;与常规目标函数相比,本文方法采用的目标函数具有更好的收敛性,其定位结果受初至拾取误差和速度模型误差的影响更小;本文提出的搜索算法能够消除由于错误拾取造成的观测到时中的异常值对定位结果的影响.