粘弹性浅圆弧形山谷地形对地震动谱特性的影响

基于凹陷地形对SV波的散射进行了分析,通过与解析解的对比,对所采用的数值算法进行了检验,然后利用钟形脉冲作为SH波输入,获得了浅圆山谷地形的地表响应,进而获得了不同观测点对自由场地点的谱比,并利用谱比研究了浅圆弧形山谷地形特征尺寸(深宽比)、波的入射角度以及介质阻尼效应等几个因素对地震动特性的影响.结果显示,浅圆弧形山谷地形的特征尺寸和地震波入射角度对地震动谱特性影响较大,影响程度与入射地震波频率有关.其中,入射角度的影响更显著.     

东亚赤道异常区电离层CT诊断──实验及初步结果

国际上第一个低纬电离层层析（ＣＴ）探测台链于１９９３年１１月建成，并开始数据采集．该ＣＴ台链沿１２０°Ｅ子午圈由６个卫星信标差分多普勒接收站组成，北起我国上海，南至菲律宾的马尼拉，纬度覆盖范围约２０°．该台链专为电高层赤道异常区电子密度大尺度结构ＣＴ成像而设计，其目的是监测研究赤道异常形态与动力学、电离层等离子体泡、行优以及暴时磁层－电离层－热层耦合过程等重要的赤道与低纬区电离层现象．本文介绍该低纬电离层ＣＴ实验的有关技术问题，包括接收台站布设、实验设备、数据采集软件和数据预处理方法．文中还给出利用ＣＴ技术，由实测斜电子含量数据反演电离层电子密度二维分布图像的几个实例．结果表明，利用低纬电离层ＣＴ台链可以有效地监测电离层赤道异常及其它低纬电高层电子密度大尺度结构的时空变化．     

基于遥感反射率的太湖优势藻识别方法

我国淡水湖库频发水华,不同类群形成的水华特征、危害及其治理方法差异显著,因此,如何区分不同藻种的遥感反射率特征,获取湖泊优势种信息是一个亟待解决的科学问题。研究基于室内藻种培养实验,培养了富营养化湖泊中的典型蓝藻和绿藻藻种,其中,蓝藻包括铜绿微囊藻（Microcystis aeruginosa）、假鱼腥藻（Pseudanabaenasp．）和束丝藻（Aphanizomenonsp．）,绿藻包括小球藻（Chlorellasp．）以及四尾栅藻（Scenedesmus quadricauda）;基于实测的遥感反射率,经归一化处理后,分析了不同藻种的遥感反射率特征,构建了DI(difference index)指数以及ADI(algae distinguish index)指数,建立了藻种分类模型,利用验证集数据进行检验,整体识别精度达77.55%,Kappa系数为0.7178。将分类方法应用于太湖野外实测遥感反射率数据集中,结果与实测的生物量数据有较好的匹配;将模型应用于太湖OLCI(ocean and land colour instrument)影像数据,获得了2019年12月和2020...     

时间-空间高阶精度矩形交错网格隐式有限差分声波正演模拟

有限差分方法因其操作简单、计算消耗低而成为地震勘探领域中最为常用的数值模拟方法之一, 然而用离散的显式差分算子数值逼近地震波动方程中的连续导数容易导致数值频散, 并且基于正方形网格离散形式的有限差分方法对不同地质模型的适应性较低.针对一阶变密度声波方程的数值模拟, 本文发展了一种适用于矩形网格离散形式的时间高阶空间隐式有限差分格式, 可以有效压制时间和空间频散, 同时灵活的网格剖分增强了其应用的广泛性.基于本文矩形交错网格时间高阶空间隐式有限差分格式的时空域频散关系和变量替换的思想, 首先采用泰勒级数展开方法求解不同方向的非轴上时间差分系数及轴上空间差分系数, 使本文差分格式可以获得任意偶数阶时间和空间精度.为了进一步提高本文差分格式在更大波数区域的空间模拟精度, 我们采用线性优化方法来求取新的轴上空间差分系数用于一阶变密度声波方程的波场迭代求解中.频散、稳定性分析及数值模拟算例表明: 相比于传统十字形空间域隐式有限差分格式, 本文矩形交错网格时间高阶空间隐式有限差分格式在精度、稳定性和效率方面均具有优势.

     

波场模拟有限差分参数优化选取

有限差分方法(Finite-difference Method, FD)广泛用于地震波场数值模拟, 但其存在固有的数值频散问题, 影响模拟的计算效率和数值精度.本文主要研究了有限差分方法的空间数值频散误差和网格划分精度以及差分算子的关系, 基于计算量最小准则, 提出了最优化有限差分参数选取流程, 为有限差分数值模拟参数选取提供理论指导.本文主要工作包括: (1) 提出了空间数值频散正变换过程(Forward Space Dispersion Transform, FSDT)方法, 该方法可以高效模拟出不同网格划分精度(波长采样点数)的带有空间数值频散的波场; (2) 提出了波场空间数值频散误差衡量准则, 可以定量地判断出数值模拟导致的波形频散程度, 选取合适的频散误差阈值; (3) 研究了给定空间数值频散误差阈值下, 差分算子系数、差分算子阶数、网格划分精度与计算量之间的关系.文中基于雷米兹交换方法(Remez Exchange Method, RE)和泰勒级数展开方法(Taylor-series Expansion Method, TE)的差分系数, 在空间数值频散误差阈值0.01时, 数值模拟了不同差分算子阶数、网格划分精度与计算量的关系, 并给出了有限差分参数选取的参考值.

     

基于纵横波分离FCT弹性波正演频散压制

有限差分方法被广泛应用于地震波数值模拟和传播.传统有限差分法采用Taylor级数展开实现空间偏导数的差分, 但该方法会因为网格离散化而产生数值频散, 降低地震波模拟的精度.优化差分系数正演方法能在一定程度上压制部分频散, 然而纵、横波速度取值差异较大, 在弹性波有限差分正演模拟中, 在满足纵波最大速度确定的稳定性条件下, 浅层低速横波波场往往会产生明显的频散现象.为了削弱弹性波场正演数值频散, 提高数值模拟精度, 本文首先采用优化差分网格系数降低数值频散, 然后再采用通量校正传输(Flux-Correction Transport, FCT)法来进一步压制弹性波场有限差分数值频散.常规的FCT法是对弹性波场直接进行频散压制, 但由于弹性波场中纵、横波速度差异明显, 横波波场频散明显强于纵波, 为了压制横波波场的数值频散, 往往需要选取较大的频散压制参数, 但这会使频散较弱的纵波产生假象.因此本文提出基于纵横波分离FCT弹性波正演频散压制方法, 对分离之后的纵横波场分别选择合适的频散压制参数进行通量校正, 可以有效压制数值频散, 削弱纵波FCT产生的假象.通过理论分析和数值算例发现, 本文方法能有效削弱弹性波场有限差分数值频散, 相对于常规FCT方法没有假象产生.

     

用交错网格有限差分法计算三维频率域电磁响应

用交错网格有限差分法(SFD),实现了三维频率域电磁场响应 的数值模拟. 该方法适用于任何方向的磁偶极子源. 经与解析方法、积分方程等 其他方法的计算结果对比表明,交错网格有限差分法结合散度校正和不完全乔累斯基分解预 处理的双共轭梯度迭代方法进行正演计算,速度快、精度高、结果稳定,能适应三维复杂介 质的数值模拟,为三维电磁反演奠定了基础.     

利用高阶交错网格有限差分法模拟地震波在非均匀孔隙介质中的传播

烃类储集层是一种复合多相介质，在固体颗粒的空隙中含有气体或液体. 研究弹性波在该类地层中的传播规律对于油气勘探开发，特别对于全波列声波测井有重要意义. 为了提高孔隙弹性介质数值模拟的计算效率，本文采用改进显式交错网格有限差分算法取代常用的空间域四阶和时间域二阶的速度 - 应力有限差分算法，算法的空间域为八阶、时间域为二阶. 虽然计算的时间步长略小于空间域四阶的情形，但高阶有限差分算法可以选择较粗糙的网格，因此补偿了计算的低效；同时高阶交错网格有限差分算法的空间频散性比低阶算法小. 利用该算法计算了一个两层模型的波场，同时还模拟了等效弹性和孔隙弹性模型中波的传播. 结果表明慢波及其影响明显，尽管慢波衰减很快，但被某一界面反射后，转换形成的P波和S波仍以正常的方式传播，且比慢波衰减小.     

求解波动方程的任意差分精细积分法

在地震偏移成像技术中，常常要反向外推波场，因此涉及到波动方程的数值求解问题。本文提出了一种求解偏微分方程的新的半解析方法——任意差分精细积分(ADPI)法。其大体思路是：空间域上作坐标离散，但不采用古典差分法的等分离散方式，而是一种相对自由的、任意的离散法；时间域上则采用子域精细积分的方法，既保留了精细积分法的高精度，又克服了工作量大、占用存储大等缺点。该方法具有精度高、带宽小、稳定性好等多项突出优点，并且可以灵括处理各类边界条件。本文从简单一维、较普遍一维、以及二维3种情形讨论该方法对波动方程的具体应用，结合实际的算例，详细分析了各种算法的可行性和精度特点。