双相介质地震波场数值模拟的迭积微分算子及其PML边界条件

将基于计算数学中Forsyte 广义正交多项式的迭积微分算子引入到地震波动方程的一阶速度--应力方程的空间微分运算中去,并采用时间错格有限差分算子替代传统的差分算子以匹配高精度的空间迭积微分算子,从而发展一种全新的地震波场正演模拟方法,来解决复杂非均匀介质模型中的波场传播问题.为了大幅衰减人工边界引起的反射,本文将完全匹配层(Perfectly Matched Layer,PML)吸收边界条件引入到所构建的方法中,以解决迭积微分算子法的边界问题.以二维波动方程为例,用迭积微分算子法实现了双相介质的地震波场正演模拟,模拟结果表明,双相介质模型较好地解释了含流体孔隙特性.同时也表明迭积微分算子法是一种非常实用、有效的数值模拟方法.     

覆盖层影响下典型地-井模型瞬变电磁法正演

地-井瞬变电磁法是用于地质找矿的有效方法之一,开展该方法的三维正演研究能对资料处理与解释提供帮助。在时间域有限差分算法的基础上,建立了三维地质模型和覆盖层模型,模拟了该模型下的地-井瞬变电磁响应;引入了参数EA,分析了方形低阻体和覆盖层的响应特征以及覆盖层影响因素。结果显示：低阻体位于发射场源下方时,其响应曲线呈现“双高一低”的极值特征,高峰值分别对应低阻体上下边界,低峰值对应低阻体中心;低阻体位于发射场源旁边时,其响应曲线呈现单极值特征,该特征与覆盖层响应特征相同。当覆盖层电阻率较低时,“趋肤效应”作用使覆盖层的响应强度强于低阻体;但当发射场源离钻孔足够远时,覆盖层影响可以忽略,低阻体响应强度更与覆盖层-低阻体电阻率的比值有关。孔中发射能有效降低覆盖层影响。     

避免异参同效的抽水井表皮层溶质运移方程推导

抽水井表皮层溶质运移控制方程受到异参同效的影响,存在参数估计的多解性问题。提出一个新的抽水井表皮层溶质运移方程(新瞬态Robin边界条件),并以径向收敛示踪试验为例构建其溶质运移模型。当抽水井为完整井时,通过拉普拉斯变换和有限傅立叶余弦变换得到模型的解析解;当抽水井为非完整井时,应用有限元法构建数值解。结果表明,表皮层溶质运移控制方程导致所估计的表皮层宽度(w)和形成层垂向弥散度(α_z)是在0.5 m≤w≤1.0 m和0.08 m≤α_z≤0.10 m范围内的任意组合,偏离实际值。反之,当表皮层佩克莱数(w和表皮层径向弥散度的比值)小于1时,新瞬态Robin边界条件能够准确反映表皮层的影响,消除w与α_z的多解性。参数估计值(w=0.31 m、α_z=0.17 m)唯一,接近实际值。新模型已成功用于野外试验。     

基于S波段双极化雷达的变分法的定量降水估计算法

基于比差分传播相移(K_DP)的降水估计算法R(K_DP)相较于传统基于水平反射率因子(Z_H)的算法R(Z_H)的表现更优。在雷达实际运行中,由于随机误差和后向散射相位(backscattering phase)的影响,可能出现负的K_DP。运用一种基于变分的雷达定量降水估计混合算法(V-RQPE)。该算法用变分拟合方法重构差分相位(Φ_DP),用一种新的稳健的边界条件求解方法,在消除随机误差的同时获得非负的K_DP,进而进行降水估计。随后我们使用2017年5月7日广州S波段雷达的回波数据和地面雨量站观测数据进行验证,同时使用了六种不同的算法进行对比,结果显示,在1小时累计降水估计中,V-RQPE表现最好,在24小时累计降水估计中,V-RQPE和基于变分拟合的K_DP的降水估计算法(R-VKDP)表现最好,实验结果表明变分拟合方法对雷达降水估计能力有显著提升。      

基于PML井间地震纵横波分离的弹性波数值模拟

井间地震方法观测到的波场信息丰富且复杂,其难点是如何识别和分离这些波场信息。这里从纵波波场为无旋场,横波波场为无散场出发,利用一阶速度～应力方程的交错网格高阶有限差分法,对井间弹性波纵、横波场分离问题进行了数值模拟。经数值模拟试验表明,同二阶的位移方程相比,使用等价的波场分离交错网格的一阶速度应力方程,无论是模拟精度还是模拟的稳定性方面都有了较大的提高,同时该横波资料的获取方法为井间地震资料的解释应用,提供了一种新的技术思路。     

淮南煤田新集一矿1煤组水文地质条件分析

在分析矿区及邻矿水文地质资料的基础上,结合所施工的8个水文地质钻孔的单孔抽水和群孔抽水试验,对矿区含(隔)水层的水文地质条件以及井田边界的进行了研究,认为当矿井开采1煤组时,顶板砂岩水和底板岩溶水为其直接充水含水层,砂岩水以静储量为主,它不会对矿井开采1煤组造成大的危害;底板奥灰岩溶水与太灰水之间联系较紧密,前者补给后者,是矿井突水是防治水研究的重点。得出了本矿1煤组的水文地质条件是岩溶裂隙水以底板突水的形式向矿床充水为主,水文地质复杂程度为中等类型的结论,从而为矿井水害的防治与矿井设计提供了依据。     

VTI介质中准P波方程叠前逆时深度偏移

根据具有垂直对称轴的横向各向同性（VTI）介质中的一阶准P波方程,导出了该方程在交错网格中逆时延拓的高阶有限差分格式,给出了其稳定性条件,采用完全匹配层吸收边界条件解决边界反射问题,分别应用下行波最大能量法和归一化互相关成像条件, 实现了VTI介质中准P波方程的叠前逆时深度偏移．各向异性Marmousi模型的试算结果表明,VTI介质准P波方程叠前逆时深度偏移算法不受地下构造倾角和介质横向速度变化的限制,对复杂模型具有良好的成像能力;应用归一化互相关成像条件能得到更好的成像效果．对比该模型的各向异性和各向同性逆时偏移剖面表明,在各向异性地区采集的纵波数据用各向异性偏移算法理论上能得到更好的成像结果．      

有限差分法地震波数值模拟的几个关键问题

地震波数值模拟是研究地震波传播规律的重要手段,对指导地震资料的采集、处理和解释有重要意义.有限差分法由于其简单易于数值计算,在求解偏微分方程和近场波动的问题越来越受到关注.本文详细阐述了与有限差分法地震波数值模拟有关的几个重要问题,对有限差分法地震波数值模拟的发展历史、波动方程的差分离散、震源处理、稳定性分析、边界条件...     

预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

位移控制的子结构地震模拟振动台混合试验方法

为了解决地震模拟振动台承载能力及台面尺寸对大型结构试验的限制,扩展振动台的功能,提出了位移控制子结构地震模拟振动台混合试验方法,包括试验原理、试验过程及数值积分方法,并给出了2种子结构边界条件的模拟形式.通过1个简单框架结构的地震模拟振动台试验和子结构混合加载试验验证了该方法的可行性,并指出了该试验方法的主要技术问题.混合试验方法通过子结构技术和振动台试验相结合,解决了目前的地震模拟振动台试验和拟动力试验在设备规模和加载速度上的局限性.