三角形网格下求解二维浅水方程的和谐Godunov格式

为保证计算格式的和谐性,通过特殊的底坡源项处理技术,在三角形网格上建立了求解二维浅水流动方程的具有空间二阶精度的Godunov格式。应用准确Riemann解求解法向数值通量,用改正的干底Riemann解处理动边界问题。经典型算例和钱塘江河口涌潮计算验证,表明模型健全,分辨率高,具有较大的推广应用价值。     

多层弹性平面问题解的精确刚度矩阵法

从平面弹性力学的基本方程出发,利用Fourier积分变换等数学手段,推导出了单层平面问题的刚度矩阵,然后按有限元法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程和Fourier积分逆变换,得到在任意静荷载作用下多层弹性平面问题的精确解。由于刚度矩阵不含有正指数项,计算时不会出现溢出现象,从而克服了传递矩阵法的缺点。推导过程中摒弃了应力函数的选择,使得问题的求解更加合理化,同时也为进一步研究此类问题如温度场、动力学等奠定了理论基础。计算实例证明了推导结果的准确性。     

基于Origin的渗透系数衰减方程在地热水回灌中的应用

孔隙型热储层的回灌过程中,由于物理堵塞、生物化学堵塞、悬浮物堵塞和气体堵塞,热储层的渗透性下降,回灌能力变弱.为了定量查明渗透系数的衰减,以辽南地区某回灌试验为例,使用Origin软件计算出研究区渗透系数衰减方程,在此基础上计算了回灌井灌压-时间关系曲线和回扬周期.研究表明,利用Origin软件求取渗透系数衰减方程具有...     

一种快速求解程函方程的间断伽辽金方法

作为一种在勘探地震学里近期发展起来的有限元类数值求解程函方程的算法,间断伽辽金方法通过拓展基函数空间提高数值精度,但存在计算效率低的问题.为了解决该问题,提出一种并行快速扫描间断伽辽金方法.算法在快速扫描过程中对节点的计算顺序进行Cuthill-McKee变换,将能够并行计算的节点分为不同任务集合,进而实现不同任务集的...     

蒸发波导环境下风浪对电磁波传播影响的数值模拟研究

蒸发波导是一种特殊的大气波导,在其中传播的电磁波信号会被陷获在近海大气层中,实现超视距传播。受海表面温度、湿度、风速、微波频率等因素的影响,海洋蒸发波导环境中的微波传播特性起伏变化很大,规律十分复杂。以往的工作主要通过计算这些气象因素对蒸发波导条件下大气折射率剖面的影响来分析它们对路径损失的作用,其结果与实验数据仍有较大差异。本文在一定的蒸发波导条件下,利用一维分形海面模型产生海面"地形",将其作为抛物方程电磁波传播模型的边界条件进行计算,得到相应的路径损失,并与传统计算方法进行对比,分析了不同蒸发波导高度、不同频率及不同接收天线高度时的数值模拟情况,可为舰艇通信系统或者雷达系统的设计提供相应的依据。     

海州湾底层鱼类物种多样性与生物量的关系

随着生物多样性的丧失,全球许多海洋生态系统的功能发生了显著变化。为了解生物多样性和生态系统功能之间的关系,科学的海洋生态保护与管理显得尤为重要。本研究根据2013−2022年春季在海州湾及其邻近海域进行的底拖网调查数据,基于结构方程模型（SEM）探究了底栖鱼类群落中环境因素、生物多样性（物种丰富度与均匀度）和生态系统功能指标（以总生物量表示）之间的关系。结果表明：物种丰富度和生物量之间存在显著的正相关,而均匀度和生物量存在显著的负相关。环境因素中盐度对物种丰富度和生物量均有显著的影响,而对于温度来说,夏季与冬季的温度比年平均温度对生物量的影响更强烈。本研究表明,生态位互补效应和选择效应这两种机制可能同时对维持海州湾底层鱼类群落中的生物多样性−生物量关系发挥了作用,此外这种关系还依赖于它们所生存的环境和栖息条件。     

关于频率电磁测深几个问题的探讨(八)——频率电磁测深与瞬变电磁测深的关系

从探测原理出发,说明人工源频率电磁测深与瞬变电磁测深的关系密切,这是一个问题的两个方面,可以互相转换。两种探测方法实质是一样的,具有共性;但毕竟从频率域和时间域两个侧面反映问题,又具有各自特性。两种方法的原理所依据的麦克斯韦方程,其正演公式可通过傅立叶变换相互转换;从实际探测所发收的信号特征,场区划分与探测深度等主要方面阐明两种探测方法异同,以加深对探测理论与方法的认识,结合实际择优应用。     

裂缝参数对考虑滑脱效应页岩气水平井产能的影响

低渗透页岩气藏中,气体渗流时会受滑脱效应的影响。建立了考虑滑脱效应的气、水两相页岩气藏渗流数学模型,并建立了理想地质模型,采用数值模拟方法,研究了水力压裂的不同裂缝参数对水平井产能的影响。模拟结果表明:裂缝条数、长度和间距是影响页岩气井产能的重要参数,而裂缝宽度和渗透率对产能的影响相对较弱;页岩气井的产能随着裂缝条数和裂缝长度的增加而增大;水平井的水平段长度及裂缝条数一定时,可通过增大裂缝间距来减少裂缝间的相互干扰。     

求解弹性波方程的辛RKN格式

将弹性波方程变换至Hamilton体系, 构造适用于弹性波模拟的高效显式二阶辛Runge-Kutta-Nyström (RKN)格式, 运用根数理论得到此格式的阶条件方程组. 通过给定系数的限定条件, 得到方程的对称解. 为了使时间离散误差达到极小,提出数值频率与真实频率比较,通过Taylor展开,得到关于辛系数的限定方程,求解方程组得到最小频散辛RKN格式. 对比分析时间演进方程的稳定性,得到使库朗数达到极大值的限定方程,求解方程组得到最稳定辛RKN格式. 发现此两种格式为同一格式. 新得到的辛RKN格式不依赖于空间离散方法,为了对比的需要,选取有限差分法进行空间离散. 在频散、稳定性分析中,与常见辛格式对比,从理论上分析了本文提出的格式在数值频散压制、稳定性提升等方面的优势, 数值实验进一步证实了理论分析的正确性.     

求解声波方程的径向基函数无网格法

地震勘探广泛应用于油气、煤田勘探。地震波场数值模拟是整个地震勘探数据处理技术的基石。将径向基函数（RBF）引入地震声波波场数值模拟中,在空间上用径向基函数无网格法来构造二阶导数,而在时间上采用简单的二阶差分公式,并重点讨论了形状参数c对该方法精度的影响,总结c经验取值范围为2~4倍平均数据点间距。设计不同模型,利用径向基函数无网格法进行声波波场模拟,并与空间四阶时间二阶的有限差分计算结果进行对比,结果表明：同样精度下,径向基函数每个波长所取的数据点数远小于空间四阶矩形网格有限差分每个波长所取的网格点数,即径向基函数的空间采样率更低,这表明径向基函数具有更小的数值频散。