CHARACTERISTICS OF TECTONIC STRESS FIELD OF THE XIAOWAN RESERVOIR BEFORE AND AFTER THE IMPOUNDMENT

Using the seismic waveform data of Xiaowan seismic network and Yunnan seismic network, we determined the focal mechanisms of 36 earthquakes(M_L ≥ 3.0)from Jun. 2005 to Dec. 2008 and 51 earthquakes(M_L ≥ 2.5)from Jan. 2009 to Dec. 2015 by generalized polarity and amplitude technique. We inverted tectonic stress field of the Xiaowan reservoir before impounding, using the focal mechanisms of 36 earthquakes(M_L ≥ 3.0)from Jun. 2005 to Dec. 2008 and CAP solutions of 58 earthquakes(M_L ≥ 4.0)collected and the solutions in the Global Centroid Moment Tensor(GCMT)catalog; We inverted local stress field of the reservoir-triggered earthquake clustering area, using 51 earthquakes(M_L ≥ 2.5)from Jan. 2009 to Dec. 2015. Focal mechanisms statistics show that, the Weixi-Qiaohou Fault is the seismic fault. Focal mechanisms were strike-slip type in initial stage, but normal fault type in later stage. Focal depths statistics of 51 earthquakes(M_L ≥ 2.5)show that, the average value of focal depths in period Ⅰ, period Ⅱ and period Ⅲ are 8.2km, 7.3km and 7.8km respectively and the standard deviations are 4.3km, 3.5km and 6.0km respectively. The average value of focal depths is basically stable in different period, only the standard deviation is slightly different. Therefore, there is not positive connection between focal depth and deviation of focal mechanisms. What's more, there are 2 earthquakes(number 46 and number 47 in Fig.5 and Table 3)with almost the same magnitude, epicenter and focal depth, but they have different faulting types as normal and strike-slip. The focal mechanism of event No.46 is strike:302°, dip:40° and rake:-97° for plane Ⅰ, however, the focal mechanism of event No.47 is strike:292°, dip:82° and rake:140° for plane Ⅰ. Likewise, earthquake of number 3 and number 18 have similar characteristic. Therefore, the obvious focal mechanism difference of similar earthquake pair indicates the complexity of Weixi-Qiaohou Fault. Considering the quiet-active character of reservoir-triggered earthquakes, we discussed the change of local stress field in different period. The σ₁ of tectonic stress field was in the near-south direction, with a dip angle of 14° before the impoundment, however, the direction of σ₁ of local stress field changed continuously, with the dip angle getting larger after the impoundment. The direction of σ₁ of local stress field of reservoir-triggered earthquake clustering area is close to the strike of Weixi-Qiaohou Fault, and reservoir impoundment increased the shear stress in the fault, so the weakening of fault was beneficial to trigger earthquakes. Comprehensive analysis suggests that fluid permeation and pore pressure diffusion caused by the water impounding, and the weakening of fault caused by local stress field are the key factors to trigger earthquake in the Xiaowan reservoir.     

南黄海西部2007年4月的逆温跃层

利用南黄海西部2007-04的温盐实测资料,采用海洋层结谱表达法及自适应识别,得到逆温跃层的"五点三要素",形成强度要素平面分布图.分析表明,逆温跃层的存在与黄海暖流水有直接的关系:1)4月份,黄海暖流水受到的海面冷却仍是产生逆温跃层的普遍原因,在该海区黄海暖流向北延伸和向两侧拓展的区域都有该种类型的逆温跃层存在,位置相对较浅;2)但在偏南的黄海暖流主干区,海面冷却产生的效应被主流区的热量补充所抵消,逆温跃层很弱甚至消失,这是该月份逆温跃层分布区向北退缩并在南部中心附近呈现缺失区的主要原因;3)南下的鲁北沿岸流水的冷水叠加在黄海暖流水的暖水上方,使逆温跃层加强,使得冷暖水的作用区成为强逆温跃层区;4)黄海暖流左侧冷沿岸流水及右侧冷水的前端向黄海暖流楔入,其前端往往覆盖在底层高温高盐的黄海暖流水上方形成下逆温跃层,从而形成双逆温跃层.这些特点,较以前认知更加客观、全面、细致和准确.     

Joint PP and PS AVO inversion based on Zoeppritz equations

Considering Zoeppritz equations, reflections of PP and PS are only the function of ratios of density and velocity. So the inversion results will be the same if the ratios are the same but values of density, velocities of P-wave and S-wave are different without strict constraint. This paper makes efforts to explore nonlinear simultaneous PP and PS inversion with expectation to reduce the ambiguity of AVO analysis by utilizing the redundancy of multi-component AVO measurements. Accurate estimation of ratio parameters depends on independence of input data. There are only two independent AVO attributes for PP reflectivity (i.e. intercept and gradient) and two for PS reflectivity (i.e. pseudo-intercept and pseudo-gradient or extreme amplitude), respectively. For individual PP and PS inversion, the values of least-squares objective function do not converge around a large neighborhood of chosen true model parameters. Fortunately for joint PP and PS inversion the values of the least-squares objective function show closed contours with single minima. Finally the power function fitting is used to provide a higher precision AVO attributes than traditional polynomial fitting. By using the four independent fitting attributes (two independent attributes for PP and PS respectively), the inversion of four ratio parameters (velocities and densities) would be estimated with less errors than that in traditional method.     

基于逆算子估计的AVO反演方法研究

传统反演算法以优化算法为主,而基于逆算子估计的AVO反演算法则利用了直接求逆的思路.算法的关键在于寻找存在逆函数的子域,进而可以在子域内直接求逆,这种解决反问题的思路不同于一般的优化类算法所采用的直接搜索解的方式,具有更高的效率.AVO反演利用了振幅随着偏移距的变化特征,反演的精度受到地震资料质量的影响,通过加入L1范数约束以及合理的初始模型有助于提高反演的稳定性以及准确度.模型测算和实际应用表明,基于逆算子估计的AVO反演方法具有较高的精确程度和可靠性.     

A method of obtaining a velocity-depth envelope from wide-angle seismic data

Summary. Due to the non-uniqueness of traveltime inversion of seismic data, it is more appropriate to determine a velocity-depth ( v-z ) envelope, rather than just a v-z function. Several methods of obtaining a v-z envelope by extremal inversion have been proposed, all of which invert the data primarily from either x-p , or T-p , or both domains. These extremal inversion methods may be divided into two groups: linear extremal and non-linear extremal. There is some debate whether the linearized perturbation techniques should be applied to the inherently non-linear problem of traveltime inversion. We have obtained a v-z envelope by extremal inversion in T-p with the constraint that the inversion paths also satisfy x-p observations. Thus we use data jointly in r-p and x-p , and yet avoid the linearity assumptions.
This joint, non-linear extremal inversion method has been applied to obtain a v-z envelope down to a depth of about 30 km in the Baltimore Canyon trough using x-t data from an Expanding Spread Profile acquired during the LASE project. We have found that the area enclosed by the v-z envelope is reduced by about 15 per cent using x-p control on the T-p inversion paths, compared to the inversion without x-p control.     

重磁解释中健全的极大似然法及三维反演

本文介绍了JOAO.B.C.SILVA提出的重磁解释中健全的(Robust)极大似然反演方法的基本原理,讨论了重力二度体非线性几何参数的反演情况,并将该方法扩展到三维重、磁非线性几何参数反演中。重点讨论了场值中有较大随机干扰、有不同类型的地质噪音源影响下的反演结果。理论反演结果表明,该方法有比最小二乘反演方法更接近实际的解,从而提高了解的可靠性,对于重磁局部异常或孤立异常有较好的反演效果。     

Crust and Upper Structure of Qinghai-Tibet Plateau and Its Adjacent Regions From Surface Waveform Inversion

CRUST AND UPPER STRUCTURE OF QINGHAI-TIBET PLATEAU AND ITS ADJACENT REGIONS FROM SURFACE WAVEFORM INVERSION     

基于二维NLCG反演的水文地质结构辨识研究

地下水调查是草原水资源合理开发利用的前提,而探明水文地质结构对地下水调查尤为重要。笔者充分利用大地电磁法（MT）和地面核磁共振（SNMR）各自的优点,尝试采用MT和SNMR对水文地质结构进行探测,并通过二维非线性共轭梯度（NLCG）反演解译了电性与地层两者的对应关系;实现了典型层状含水层模型和局部含水层模型正反演试验;通过实测数据反演与钻孔资料对比验证了该反演算法对低阻含水层地质体解译的有效性。试验结果对水文地质调查具有一定参考意义,也为同类数据反演提供了方法借鉴。     

一类双曲-抛物型方程的广义解

讨论了一类奇摄动双曲-抛物型方程广义初边值问题,在适当的条件下,用Galerkin方法研究了广义解的存在性、唯一性,同时得到了解的渐近估计式。     

控制点与检查点数量和比例对水深遥感反演精度的影响分析

在水深遥感反演半经验模型中,水深控制点和检查点选取是一项非常重要的工作,前者用于建立实测水深值与遥感影像灰度值之间的数量关系,后者用于评价水深反演精度。但前人在相关研究中并未就控制点和检查点的选取数量和比例给出有实验依据的建议,亦未开展其对水深反演精度的影响分析。针对这个问题,作者使用6020个实测水深点和World View-2、GF-1 WFV、Landsat8 OLI 3种多光谱影像,基于三波段水深反演模型开展了81组实验,分析比较了不同数量和比例的控制点与检查点对反演结果的影响,评价过程利用平均绝对误差(Mean Absolute Error,MAE)、平均相对误差(Mean Relative Error,MRE)和反映实测水深与反演水深相关性的决定系数R2等3种参数。结果表明:(1)当控制点数量达到31个时,水深反演精度即趋于稳定;(2)检查点数量在30个时,其评价指标已可以代表模型反演精度;(3)控制点和检查点的数量比例对反演精度并无影响,但控制点选取过多则会产生反演精度降低的现象,采用人工选取检查点时剔除浅水区的异常点即可有效避免这一问题。