基于背景噪声和地震面波联合反演华北克拉通中部岩石圈结构

基于ChinArray三期项目布设于华北克拉通中部的流动台阵观测数据,利用背景噪声互相关和地震面波层析成像获取了研究区内6—140 s周期的瑞雷面波频散,使用蒙特卡罗非线性反演方法获得了华北克拉通中部岩石圈的高分辨率三维S波速度结构。结果显示华北克拉通不同地块的岩石圈速度结构存在显著的横向差异：其中鄂尔多斯盆地腹地整体表现为高速特征,延伸至200 km以下,但其东南缘存在小范围的低速异常;东部的华北盆地整体表现为低速特征,具有较薄的地壳和岩石圈厚度;中部造山带南北两端以及南北重力梯度线下方存在相连接的低速区域,在深处延伸至华北盆地下方;在下地壳和上地幔顶部,大同火山群区域的低速体逐渐向西偏移至鄂尔多斯盆地东北角下方;而在上地幔中,该区域的低速异常随深度增加而逐渐减弱,低速体延伸至东南方向的华北盆地下方。基于本研究获得的S波速度模型,我们认为：鄂尔多斯盆地腹地保持了克拉通特性,但其东南缘存在局部的岩石圈改造作用;华北盆地发生了强烈的岩石圈破坏减薄和地壳伸展变形;中部造山带南北端以及南北重力梯度线下方的岩石圈发生了局部的改造减薄,其机制可能都来源于华北盆地下方地幔热物质的上涌;大同火山群下...     

井中垂直双极源体积分方程法三维模拟研究

类比水平方向电偶极源和垂直方向磁偶极源电磁场求解的原理,又利用偶极沿电流方向积分是双极原理,推导出地下均匀半空间中垂直双极源电场和磁场解析计算公式.采用体积分方程法对井中垂直双极源电磁响应进行了三维模拟.它对提高井中垂直双极源圈定储层范围的精度有着非常重要的意义,同时也为井中垂直双极源三维反演问题奠定了基础.     

波动方程角度域共成像道集

针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.     

共偏移距可分表示法叠前深度偏移

基于双平方根方程的共偏移距可分表示法叠前深度偏移用于复杂介质成像, 该方法在中点-偏移距坐标中同时向下延拓炮点和检波点波场, 实现方式采用正反傅立叶变换.构造的双平方根方程波场延拓算子能够使波数域变量与空间(速度) 域变量分离, 波数域内进行相移计算, 在空间域对因介质横向变速引起的时移作修正.地震数据偏移不需要逐炮计算, 具有较高的效率.对Marmousi模型数据的偏移成像结果显示, 该方法较好地成像强横向变速介质中的复杂构造.      

Estimation of non-point source pollution loads under uncertain information

Many kinds of uncertainties are involved, such as random, fuzzy, grey, unascertained property and so on, in soil erosion process. To exactly predict the non-point source pollution loads, some uncertainties should be taken into consideration. Aiming at the deficiency of present blind number theory being helpless for fuzziness, a novel blind number, i.e. extended-blind number, was introduced by substituting a set of triangular fuzzy numbers （TFNs）, expressed as a-cuts, for interval values in present blind number, and the expected value of extended-blind number was also brought forward by referring to the current blind number theory. On the basis of denoting the parameters of Uni- versal Soil Loss Equation （USLE） as extended-blind parameters, a novel USLE model was established for quantitatively evaluating soil erosion loss and non-point source pollution loads. As a case, the uncertain USLE was employed for predicting the soil erosion loss and non-point source pollution loads of absorbed nitrogen and phosphorus in a dis- trict in the Hangbu-Fengle River basin, in the upstream of Chaohu Lake watershed. The results show that it is feasible in theory to extend blind number into fuzzy environment and reliable on conclusion to apply extended-blind number theory for predicting non-point source pollution loads.     

Two dimensional tensor product B-spline wavelet scaling functions for the solution of two-dimensional unsteady diffusion equations

The fourth-order B spline wavelet scaling functions are used to solve the two-dimensional unsteady diffusion equation. The calculations from a case history indicate that the method provides high accuracy and the computational efficiency is enhanced due to the small matrix derived from this method.The respective features of 3-spline wavelet scaling functions, 4-spline wavelet scaling functions and quasi-wavelet used to solve the two-dimensional unsteady diffusion equation are compared. The proposed method has potential applications in many fields including marine science.     

湿饱和流中的Richardson数和不稳定的研究

本文推导了干、湿饱和流中的Ri数方程.选取华北的一次降水过程,利用模式输出资料计算并对比分析干空气、湿饱和流中Ri方程中的各项,发现干过程中水平风速的垂直切变对干Ri数的变化和不稳定的影响占主导地位,而在湿过程中,湿饱和Ri的分子(Brunt_Vaisala frequency, 简称BVF)对Ri数的变化及不稳定的影响可达到与分母相当的量级. 既然BVF在湿过程中对不稳定有非常重要的影响,我们寻找了更符合湿过程中实际情况的BVF表达式,以对Ri进行修正,从而更好地判断湿过程的不稳定.研究发现,由于考虑了气块的虚温效应和总的水物质混合比的变化,即考虑降水过程中部分液态水脱离气块,有效地减小了雨区的稳定性,使得修正后的湿Ri公式可能更适合于诊断雨区的不稳定.     

三维波动方程基准面校正方法的应用研究

波动方程基准面校正处理被认为是当地表高程变化剧烈、地表一致性假设又不成立的情况下对于常规高程基准面校正的必要替代. 在二维波动方程基准面校正方面,前人已经作了大量工作,并且在实际应用中获得了良好的效果,应该将其进一步推广到三维. 本文采用三维频率空间域有限差分波场延拓算子以“逐步－累加”的方式实现了三维波动方程基准面校正,并对实际数据进行了处理. 对于西部某三维数据的实际计算结果表明,相对于传统的模型法高程静校正,波动方程基准面校正更合理地实现了基准面校正,有助于提高后续的速度分析精度和成像质量.     

傅里叶有限差分法保幅叠前深度偏移方法

地震数据中饱含有丰富的走时信息和振幅信息. 为解决传统偏移方法中几何扩散和入射角变化引起的振幅误差问题,本文提出了一种实用的波动方程保幅地震偏移方法. 该方法从全声波方程出发进行单程波保幅分解,得到直观、高效率的直接面对地震波传播波场的压力分量进行延拓的保幅偏移单程波方程,进而推导出一个含有6项的傅里叶有限差分法保幅偏移的算子方程;修改边界条件和成像条件,使修改后的边界条件和成像方程中考虑振幅补偿,从而从三方面补偿几何扩散损失和入射角变化对振幅的影响. 脉冲响应测试、单炮记录的数值试验以及SEG/EAGE盐丘模型的叠前偏移结果表明,该方法不但可以使散射能量聚焦、归位,提高成像精度;而且可以输出正确反映地下反射系数的振幅信息,为后续的地震属性分析(如AVO/AVA)提供更真实的地震信息.     

水平层状介质中基于DTA的三维电磁波逆散射快速模拟算法

为提高水平层状介质中三维电磁波散射和逆散射数值模拟的效率,在对角张量近似(DTA)的基础上根据不同回代方式得到了求解积分方程的DTA1和DTA2两种近似. 这两种近似可以作为计算积分方程稳定型双共轭梯度快速Fourier变换(BCGS-FFT)算法的初始猜测值和预条件因子,从而形成效率更高的混合DTA-BCGS算法. 散射实例说明了DTA2的高精度和混合DTA-BCGS算法尤其是混合DTA2-BCGS算法的高效率. 由于DTA2近似程度更高,将DTA2与变型Born迭代反演方法(DBIM)相结合形成了一种对三维异常体进行重构的快速电磁波逆散射技术. 文中的逆散射实例说明所开发的逆散射技术对重构水平层状介质中的任意三维异常体是非常有效的.