感潮河网二维水流-输运耦合数学模型

感潮河网水流不仅受径流、潮汐动力的双重影响,还常受水闸调度等人类活动的影响,水动力条件极为复杂。针对水闸调度影响下感潮河网复杂水动力及其伴随物质输移扩散模拟,提出了考虑水闸调度并集成输运对流项的水流-输运通量耦合求解器,建立了基于非结构网格Godunov格式的二维水流-输运耦合数学模型。采用具有时空二阶精度的MUSCL-Hancock预测-校正格式,结合变量重构限制器技术,在保证计算精度的同时避免了数值振荡。运用斜底三角单元网格表达水闸线状地形,并通过地形调整模拟了水闸启闭过程。算例研究表明,该模型计算精度高,可有效模拟水闸调度影响下感潮河网水流运动及污染物输运过程,具有较好的推广应用价值。     

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

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

Transition from wavelets to ripples in a laboratory flume with a diverging channel

An experimental investigation on the initiation and development of bed forms on a bed of fine silica sand was conducted under alluvial flow conditions in a laboratory flume with a diverging channel. The main aims of the study were to assess： i） the steepness of bed forms in the transition stage of development; and ii） the threshold height of wavelets （ηt） that triggered the start of ripple development. Detailed bed profile measurements were carried out using an acoustic Doppler probe, traversed longitudinally over the sediment bed at various experimentation times. The bed form dimensions were extracted from such bed profile records and analysed for the wavelet, transition and equilibrium stages. It was found that the steepness of ripples in the transition and equilibrium stages were similar, confirming predictions of previous mathematical model simulations. A lognormal distribution fitted the wavelet length data. The wavelet threshold height was estimated as ηt ≈ 7 mm, or ηt≈ 80 in wall units. Such a height magnitude suggested that ripple development could be triggered by the wavelets reaching the outer flow zone of a turbulent boundary layer. The ηt value obtained corresponded generally to the intersection point between two predictive equations for bed form dimensions. A formulation was developed to predict ηt as a function of the sediment grain size, which was confirmed for the fine sand used in this study.     

三维波动方程时空域混合网格有限差分数值模拟方法

常规高阶和时空域高阶有限差分方法广泛应用于三维标量波动方程的数值模拟,这两种差分方法仅利用笛卡尔坐标系中的坐标轴网格点构建三维Laplace差分算子,相应的差分离散波动方程本质上仅具有2阶差分精度,模拟精度低.本文将三维笛卡尔坐标系中非坐标轴网格点分为两类:坐标平面内的非坐标轴网格点和坐标平面外的非坐标轴网格点,系统推导出了两类非坐标轴网格点构建三维Laplace差分算子的方法,进而提出了一种利用坐标轴网格点和非坐标轴网格点共同构建三维Laplace差分算子的混合网格有限差分方法,并利用时空域频散关系和泰勒展开建立差分系数方程,推导出了差分系数的通解.相比常规高阶和时空域高阶差分格式的2阶差分精度,时空域混合网格差分离散波动方程理论上能够达到任意偶数阶差分精度,模拟精度显著提高,同时稳定性更强.频散分析表明:相比常规高阶和时空域高阶差分格式,在计算效率基本相同时,时空域混合网格差分格式能更有效地减小数值频散,减弱数值各向异性,模拟精度更高;在模拟精度基本相当时,混合网格差分格式能采用更大的时间采样间隔,计算效率更高.数值模拟实例进一步验证了混合网格差分格式在提高模拟精度和计算效率方面的...     

《华南地质与矿产》征稿简则

《华南地质与矿产》为国土资源部武汉地质矿产研究所主办的综合性地学学术刊物(季刊,国内外公开发行),以报导地调科研成果为主要任务,促进学术和科技信息交流为目的。本刊栏目主要包括基础地质、矿床地质、环境地质、同位素地质、水文地质与工程地质、灾害地质、农业地质、测试方法与技术、矿产品开发和利用,以及综合评述和动态报导等。     

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

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

南海海底地震仪异常数据的分析和处理

海底地震仪(Ocean Bottom Seismometer, OBS)数据处理至关重要,是获取深部地壳结构的基础与前提.2006年实施OBS2006-2测线时,有2台OBS(OBS03,OBS06)数据出现异常,无法使用.由于海上航次花费巨大,采集到的数据弥足珍贵.本文采用数据格式检查、邻近台站对比分析、重采样等方法,成功地对这2台OBS数据进行了解编处理,得到了这两个台站的综合地震记录剖面;利用上述方法对2011年实施的OBS973-3测线中的异常台站OBS03进行了分析处理,同样得到了OBS03台站的综合地震剖面;通过查看两次海上实验班报发现,OBS2006-2测线之OBS06与OBS973-3测线之OBS03内部Sedis编号相同,为同一台记录仪器,再一次验证上述处理方法正确可行;然后对OBS2006-2测线2个台站进行震相识别与走时拾取后,利用前人纵波速度模型开展了射线追踪与走时模拟.此次对异常OBS数据的重新处理工作,不仅为OBS探测提供了宝贵的数据处理经验,而且将提高OBS2006-2测线地壳结构的可靠性和约束性,具有重要的研究意义.     

适用于声波方程数值模拟的时间-空间域隐式有限差分算子优化方法

声波方程数值模拟已广泛应用于理论地震计算,同时构成了地震逆时偏移成像技术的基础.对于有限差分法而言,在满足一定的稳定性条件时,普遍存在着因网格化而形成的数值频散效应.如何有效地缓解或压制数值频散是有限差分方法研究的关键所在.有限差分格式分为显式有限差分和隐式有限差分.隐式有限差分能够进一步压制数值频散效应.因此本文提出了给定频率范围满足时间-空间域隐式有限差分频散关系的方法,并根据震源频率、波速和网格间距确定波数范围,在此基础上建立方程确定了相应的隐式有限差分系数,使得差分系数能在更大频率范围符合波场传播规律.通过频散分析和正演模拟,验证了本文方法的有效性.     

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法（DG-FEM）作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔（SSP Runge-Kutta）算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层（NPML）的思想,基于复频移（CFS）拉伸坐标变换推导了一种新的PML吸收边界条件（简称为CFS-NPML）,该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.     

Excel在测量计算及报表格式中的应用

作为微软办公软件之一的Excel在各行各业都用得比较普遍而且频繁,在测量作业中的计算更不用说,很多简易平差都可以在上面完成.本文主要就其报表格式及常用的几个测量计算函数予以简述.