非均匀介质热蠕变流动的数值求解

针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格－粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov Galerkin方法,求解能量方程时采用了流线迎风Petrov Galerkin方法,网格－粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格－粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性.     

无单元法用于地震波波动方程模拟与成像

波动方程方法是解决地震正反演问题的基本工具之一.无单元法作为一种新兴的偏微分方程数值计算方法,已经在材料力学、热传导等领域取得了显著的成功.由于抛弃了单元的概念及采用滑动最小二乘的拟合方法,使得无单元法具有前处理简单、精度高、独立变量解高次连续等优点.本文首先介绍无单元法求解波动方程的原理,指出影响其精度的主要因素.在算例的基础上详细讨论了无单元法用于实际波动问题的效果,并进一步尝试利用无单元法进行地震波数值模拟和反演成像的研究.模型计算的结果表明,无单元法能够较好的处理地震模拟和成像问题,精度和稳定性是令人满意的.     

方位各向异性粘弹性介质波场数值模拟

当地震信号通过复杂地球介质时,地层除了表现为各向异性,还表现为内在的粘弹性特征.因此,为准确描述地震波在地球介质中的传播特征,理想的地球介质模型应该能够模拟岩石的各向异性特征和衰减特征.本文给出了各向异性粘弹性介质模型的波动方程及其差分格式,并利用有限差分法实现了地震波波场数值模拟.结果表明了该介质模型中地震波场特征与各向异性主轴方位和介质的粘滞性参数之间的关系.     

利用毛管压力曲线分形分维方法研究流动单元

利用取心井铸体薄片获得的图像资料和毛管压力曲线,通过图像分形几何学方法以分维数的形式定量地表征出了复杂的微观孔隙喉道结构特征,发现能够很好地划分和评价孔隙岩石中油、气、水的渗流差异,可以用于储层微观流动单元表征。文中阐述了岩石微观孔隙喉道结构分形的理论基础、计算方法和应用于表征流动单元的依据。建立了中国西部砾岩低渗透油藏微观孔隙喉道分维数与孔隙度、渗透率之间计算图版,据此在油藏中利用常规测井资料获得的孔隙度、渗透率参数计算微观孔隙喉道分维数,开展全油藏流动单元划分与评价,取得了良好的效果。研究结果表明,利用毛管压力曲线分形分维方法研究储层微观流动单元是一种很有效的途径。     

河南木厂河区银、金、铜矿床的构造控矿特征

河南木厂河地区,位于秦岭-大别山中央造山带的东部,桐柏-大别构造亚带东段。区内以长期活动的NWW-近EW向深大断裂及派生的NWW向,SN向断裂以及平行于东部郯庐断裂带的NE向断裂为主要特征,褶皱表现为形态复杂的宽缓线形。通过调查研究,确认了木厂河矿区存在三个构造岩石单元:核杂岩单元(CC)、超高压单元(UHP)、高压单元(HP).各单元之间发育有大型的拆离断层,有利于成矿元素富集。本区矿床的形成与大别山超高压变质带折返后的伸展拆离作用密切相关。矿体形成部位主要在核杂岩单元之上,超高压单元之下的拆离断层带。拆离断层对本区的矿化起了重要的控制作用。      

二维复杂介质中地震波走时和射线的计算方法及其应用

将Podvin和Lecomte的精确局部格点走时计算方法和Schneider等人的动态规划方法结合起来,可得到一种快速、有效的有限差分波前计算方法。使用该方法对各种类型地震波的走时和射线的计算进行了讨论,并给出了这种有限差分走时计算在叠前深度偏移中的应用实例。     

超参数对GRU-CNN混合深度学习弹性阻抗反演影响研究

CNN-GRU混合深度学习反演弹性阻抗取得了较好的反演效果.但是,基于深度学习的叠前反演参数众多,包括内部深度学习网络可学习参数和外部超参数等,目前超参数选取对网络性能及计算速度影响尚缺乏系统性研究,这直接影响到了该方法的进一步推广应用.因此,本文在混合深度学习反演弹性阻抗基础上,探讨学习率、Ep-och、batch_...     

中国工业烟粉尘排放时空演化及其影响因素

采用空间分析和空间杜宾模型等方法,研究1999—2017年中国工业烟粉尘排放时空演化特征及其影响因素。结果表明：① 中国工业烟粉尘排放空间分布差异明显,其排放的基尼系数和污染物分布指数均呈现下降态势,空间集中程度有所缓和。② 中国工业烟粉尘排放空间分布呈东北?西南走向,其排放中心不断由东南向西北方向迁移。③ 中国工业烟粉尘排放存在空间相关性和空间溢出效应,其冷热点区空间分布发生显著变化。④ 能源消耗、第二产业比重、人口密度和经济发展水平的提升将会增加工业烟粉尘排放,而外资水平、治理技术水平和环境规制力度的提升则有利于减少工业烟粉尘排放。     

有限尺度李雅普诺夫指数与海表温度梯度相关性的初步分析

为了探究海表温度和海面高度之间的瞬时相关性,介绍了一种卫星测高数据的拉格朗日分析指数--有限尺度李雅普诺夫指数(Finite Size Lyapunov Exponent,FSLE),以黑潮延伸体区域的涡旋和南大西洋的亚南极锋为例,通过对观测、模式结果和融合产品结果的分析,探讨了该指数与海表温度梯度(Sea Surface Temperature Gradient,SSTG)之间的相关性。比较FSLE图像和SSTG图像发现,FSLE与SSTG均呈丝状结构,对海洋表层水体结构描述具有一致性,尤其在温度梯度大和地转流强的区域更为一致。二者的一致性要远好于其他常用方式,比如全流速、OW参数涡旋识别方法和Winding-Angle涡旋识别方法。不同区域FSLE与SSTG之间的相关性表现不同,黑潮延伸体区域相关系数存在显著的季节变化,而南大西洋亚南极锋区域季节内变化突出。     

砾石单元法与SPH耦合模型在数值波浪水槽中的应用研究

为研究以流体粒子描述波浪运动,以固体单元描述砾石运动的两相介质大变形运动,在港口、海岸工程科学研究中具有重要意义。本文提出砾石单元法（GEM）,介绍了光滑粒子动力学方法（SPH）和GEM的基本原理,阐述了GEM与离散单元法（DEM）的异同之处,说明了采用SPH方法与GEM构建波浪砾石耦合运动数学模型的方法和过程。应用SPH方法建立数值波浪水槽,用GEM模拟波浪作用下堆积砾石的滚落、坍塌变形,构建了SPH方法与GEM耦合数学模型。模拟了水槽造波和波浪生成过程和波浪作用下砾石的滚落、坍塌变形,并与物理模型试验成果进行了比较,结果基本吻合。本文提出的GEM法具有模拟单相堆积砾石运动和堆积砾石与流体粒子耦合多相介质运动的功能,是对DEM法的补充和改善。本文提出的堆积力学球概念和拟序排列求解方法是砾石单元法的重要组成部分。