拟空间域弹性波方程交错网格有限差分逆时偏移

弹性波逆时偏移是一种高精度的复杂构造地震成像方法.然而,在传统的基于矩形网格离散化的逆时偏移中,介质界面通常会产生畸变.另外,因使用双程波动方程进行波场延拓,其产生的反射波会在成像过程中产生偏移假象.为解决这些问题,本文提出了一种拟空间域弹性波方程高阶交错网格有限差分格式,并给出了差分格式的稳定性条件,进而实现了高精度的拟空间域弹性波方程有限差分逆时偏移.模型实验表明,若在计算拟空间域采样间隔时引入速度界面信息,则拟空间域弹性波方程高阶交错网格有限差分逆时偏移能够避免常规弹性波方程逆时偏移中弯曲界面形态畸变问题;此外基于该方法进行波场延拓时可有效压制弯曲界面的假散射现象,并能有效压制层间反射波,因此可以减少剖面上的偏移假象,从而显著提高成像的质量.     

电性各向异性地层频率域电磁响应模拟

在现有计算机内存和运算速度条件下,电性各向异性地层中频率域电磁响应的数值模拟是一个难点. 为探索各向异性介质模型电磁响应的有效算法,同时为进一步研究各向异性地层测量响应的校正方法奠定基础,从Maxwell电磁响应方程出发,使用由三个彼此垂直的发射线圈和三个彼此垂直的接收线圈组成的线圈结构模型,应用交错网格有限差分法推导了三维频率域电磁响应的差分计算格式,采用预条件双共轭梯度法对所形成的复系数矩阵进行了迭代求解,从而得到电性各向异性地层模型下的线圈测量响应. 应用这些方法进行了多种条件下电磁响应的计算,计算结果表明,所形成的计算方法是准确可靠的,预条件双共轭梯度法的运算效率也较高;电性各向异性地层的电磁响应与地层倾角（或井斜角）、仪器测量方位、地层的水平电阻率和垂直电阻率有关;根据水平方向磁场强度分量可以确定地层界面;在电磁响应计算过程中,背景电导率的取值要尽量接近整个计算区域的平均电导率.     

基于最小二乘的三维正交介质弹性波高精度正演模拟

本文借助线性滑动理论,将裂隙密度、纵横比等不同物性参数等效为正交介质刚度矩阵,建立两组正交直立裂隙介质正演模型.利用基于最小二乘法优化的高精度交错网格高阶有限差分法模拟弹性波在正交介质中的传播过程.模拟结果表明介质各向异性强度与裂隙物性直接相关,正演波场和共炮点显示相比普通交错网格,最小二乘优化方法可以压制数值频散现象,明显提高模拟精度,子波主频较高和模型波速较低时改善效果更为突出.     

用于声波正演模拟的时空域高精度交错网格有限差分方法

地震正演模拟是逆时偏移和全波形反演中的核心问题之一,因为它们都需要高效、高精度地模拟波场正向和反向传播。为了提高数值模拟的精度,人们广泛采用高阶有限差分方法,但是大多数方法仅在空间上具有更高的精度,在时间上只有二阶精度。首先系统介绍时空域高精度交错网格有限差分方法的基本原理,然后利用模型验证方法的有效性,结果表明:时空域高精度交错网格有限差分方法拥有比常规交错网格有限差分方法更低的数值频散。      

黏滞声波地震波场模拟的SG-ONADM混合方法

交错网格方法（SG）和最优近似解析离散化方法（ONADM）是两类典型的地震波场数值模拟方法.两类方法各有其优势,相对于ONADM方法,SG方法在单个时间层内的计算更为简单;相对于SG方法,ONADM方法可以在较大空间步长条件下有效压制数值频散.结合两种方法的优势,本文提出了一种新的地震波场模拟方法（SG-ONADM）.该方法对控制方程中的一阶偏导数采用SG方法给出的一阶偏导数近似公式,高阶偏导数采用ONADM方法给出的高阶偏导数逼近公式.理论分析及数值算例表明,SG-ONADM方法保留了两种方法的优势,不仅能在较大空间步长条件下有效压制数值频散,同时具有较低的内存需求量;在对同一计算区域进行波场模拟时,SG-ONADM方法的计算效率要高于SG方法和ONADM方法.最后,我们使用SG-ONADM方法进行黏滞声波波场模拟,研究了黏滞声波在复杂介质中的传播.     

分层分类和多指标结合的西北农牧交错带植被信息提取

参照《中国植被》中的植被分类体系,结合野外考察结果,建立了适合中国西北农牧交错带的植被分类体系。以覆盖研究区的多幅Landsat影像为基础,按“分层分类,逐层验证”的思路,实现了对研究区植被信息的提取。提取时,先利用完全约束的最小二乘模型对遥感影像进行混合像元分解方法,将整个研究区划分为植被区和非植被区;在植被区,基于光谱特征、纹理特征和地形特征,构建CART决策树,获得了乔木林、灌丛和草原等7种主要植被型组;在植被型组内,基于不同植被类型NDVI的季节差异特征,构建NDVI差值比值指数 (NDVI_DR),将乔木林和灌丛区分为常绿和落叶植被型,使用温度植被干旱指数(TVDI),将草原进一步区分为荒漠草原、典型草原和草甸草原3种类型,从而得到各个植被型的空间分布范围。经验证,最终分类的总体精度能达到79.51%,kappa系数为0.773。采用的分类方法充分利用了遥感数据既有的光谱信息和纹理信息,同时辅以地形信息。实践结果表明,分层分类和多种指标相结合的方法可以有效实现对影像跨幅的、以复杂镶嵌结构为主要特征的农牧交错带植被信息提取,精度较高,技术可行。     

长春市黑土城乡交错区重金属地球化学特征及生态风险评价

以位于东北典型黑土区的长春市城乡交错区为研究区,以土壤-作物系统中Cu、As、Pb和Zn为研究对象,采用Pearson相关性分析及土壤和作物综合质量指数对土壤-作物系统地球化学特征及生态风险进行评价。结果显示,表层土壤重金属存在明显的积累现象,作物籽实重金属含量均未超标,作物籽实对土壤中Zn和Cu的吸收富集能力显著高于As和Pb。不同重金属在土壤-作物系统中的迁移能力不同可能与作物种类、pH值、有机质、N、P和K₂O含量有关。作物重金属生物富集系数空间分布特征与作物籽实重金属含量空间分布特征相似,与土壤重金属含量空间分布特征存在显著差异。IICQ介于0.13～4.44,均值为0.62。研究区耕地土壤环境质量整体良好,仅中部存在小范围呈圈层分布的风险区域。     

交错网格下的Boussinesq水波数值模型

针对一组近似到二阶完全非线性,四阶色散的Boussinesq方程,在交错网格下建立了数学模型.计算中时间层不交错,模型的求解利用混合四阶Adams-Bashforth-Moulton格式的有限差分法.数值模拟了波浪在潜堤上的演化过程,再现了波浪的浅化、反射以及非线性波能量传递等现象.对数值计算结果采用Friouer变换...     

Kelvin—Voigt黏弹性介质地震波场数值模拟与衰减特征

利用高阶交错网格有限差分模拟Kelvin-Voigt黏弹性介质中传播的地震波,同时将完全匹配层吸收边界条件引入到其边界处理中。数值模拟结果表明,完全匹配层吸收边界效果好,高阶有限差分能模拟得到的黏弹性介质波场精度较高。对模拟的黏弹性波场进行分析,表明介质的粘滞性使地震反射波的能量变弱,高频衰减明显,并比低频衰减得快,主频向低频方向移动,有效频带变窄,即降低了地震波的分辨率;并且反射转换波比反射纵波要衰减得快;而且还随着传播距离的增加,其峰值频率也逐渐降低。通过数值模拟分析具有不同的粘滞系数介质对地震波的吸收和衰减,结果表明随着粘滞系数的增大,地下介质对地震波的吸收衰减更明显。     

应用二阶完全非线性Boussinesq方程模拟破碎波浪

建立了基于高阶Boussinesq水波方程的一维波浪破碎数值模型。基于一组具有二阶完全非线性特征的Boussinesq水波方程,建立了交错网格下的高精度差分格式,推导了适用于该组方程的永形波解析解,其和松弛造波技术相结合实现了数值波浪水槽中(强)非线性波浪的无反射入射。通过模拟封闭容器内水体晃动问题对数值格式进行了验证,通过模拟孤立波在斜坡海岸上的浅化过程说明了将方程从弱非线性发展到完全非线性的必要性。采用涡粘方法处理波浪破碎,利用物理模型实验数据,分析了模型中各波浪破碎参数对数值结果的影响并对参数进行了率定。应用该模型对规则波在斜坡海岸上的传播、变浅以及破碎过程进行了数值模拟研究,数值结果同实验数据吻合良好,验证了模型的有效性。