SIMULATION OF POLLUTANTS IN RIVER SYSTEMS USING FINITE DIFFERENCE METHOD

1 INTRODUCTION Numerical computation provides an easily extended and user-friendly environment with computer aided programming for the simulation of pollutants in river systems. In the most of water quality assessment and monitoring problems during water pollution control and environmental impact assessment studies of river systems, mathematical modeling has been playing a key role for the last two decades. The majority of existing water quality models are the mechanistic and are based o…     

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

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

非均匀弹性介质地震波动模拟的Fourier变换方法

一、前言 地震模型技术是应用地震法研究地下储油地层分布的基础工作,它可以用来研究地,震波的传播以及各种波的特性,可以通过修改地质模型,反复进行正演,实现对实际地震剖面的最优逼近,以实现地震反问题。本文在文献[1]的基础上,应用Fourier变换方法研究了用非均匀弹性介质波动方程正演问题的解来模拟地震波的传播,探索了该方法对非均匀弹性介质地震模型的有效性及实用性。     

用地震反射波定量解释煤层厚度的方法

一、引言 工业可采煤层厚度一般为1.5m,如何从地震资料上准确圈定1m厚的煤层边界、并定量解释煤层厚度、准确计算地质储量,国外一些学者研究了薄层地震反射波的动力学特征与薄层厚度的关系后,得出了不同的薄层垂直分辨率标准,Widess用零相位子波作实验时,发现当薄层厚度为λ/8(λ是地震子波的主波长)时,反射波形正好是入射波的导数,并可直观地鉴别顶底反射;Kallweit利用可控震源研究了薄层响应的频谱后,提出了分辨地震波的极限为1/(1.4f)(f为地震子波的上限频率),经主频换算后,此     

STUDY ON SEISMOGENIC PROBABILITIES OF POTENTIAL HYPOCENTER REGIONS IN THE CENTRAL QINGHAI-TIBET PLATEAU USING A TECTONOPHYSICAL METHOD

A deterministic-cum-probabilistic tectonophysical method is used to estimate the hypocentral parameters and seismogenic probabilities of the potential hypocenter regions before and after the Nov. 5. 1988 earthquake (Ms =7) west of Tuotuoheyan in terms of the tectonic feature, stratigraphic structure, geoth-ermal environment, and driving background in the central Qinghai-Tibet plateau. As a result, the earth-quake risk tendency of this area is discussed     

小波变换在识别储层流体性质中的应用

根据小波变换的特点提出一种识别储层流体的新方法. 本文假设测井所得的信号为地层微观孔隙结构、地层流体性质及岩性和岩相等各部分的综合贡献之和，不同的储层中小波分析的能量谱特征不同，因而通过对已知储层模型进行小波能量谱分析，可得到各种储层的特征能量谱，根据这些储层特征能量谱就可以对待划分的储层段进行储层流体性质判别.     

反演地层电阻率的单调同伦法

本文从电法测井的数学模型出发,给出了一种反演地层参数的单调同伦法,把非线性方程组的求解转化成常微分方程初值问题的数值求解,从而给出一种大范围收敛的反演方法.最后的数值结果表明,本文给出的算法是十分有效的.     

中心回线瞬变电磁法2.5维有限单元算法

瞬变电磁法的野外资料解释基本上停留在一维水平，本文从生产实践中常用的中心回线法入手，基于电磁场本身的叠加原理，从麦克斯韦方程组出发，导出了中心回线瞬变电磁2.5维二次场（纯异常）的有限单元计算公式. 该算法采用三角形有限元网格，在尽可能拟合地下电性断面的情况下减少有限元网格的节点数和单元数；用选主元的LU分解法求解线性方程组，做到了在移动场源时只需改变右端项；在反傅氏变换中，使用新的波数选取方案，让波数m随时间t滑动. 最后给出的算例表明，该算法直接计算异常场，计算速度快、精度高.     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.     

台湾费丛峡谷对强地面运动的影响研究──利用单台资料及2．5维SH混合方法

本文通过格林函数反褶积方法，由台湾峡谷附近的记录资料预测峡谷区的强地面运动的时程曲线．峡谷区的理论格林函数应用２．５维ＳＨ混合方法求解．通过格林函数反褶积方法得到的峡谷区费丛１及费丛２台的位移、速度、加速度和实际资料对比，取得了满意结果．计算结果表明，峡谷底部的峰值加速度相对峡谷边缘为最小；在靠近震中的一侧，峡谷的加速度的最大振幅相对比另一侧大．还给出了其余３个台的预测结果，研究了它们的加速度傅里叶谱和反应谱．     

TWO－DIMENSIONAL MATHEMATICAL MODEL FOR STUDYING EROSION AND SEDIMENTATION OF A WANDERING STRETCH IN THE LOWER YELLOW RIVER

ＴＷＯ－ＤＩＭＥＮＳＩＯＮＡＬＭＡＴＨＥＭＡＴＩＣＡＬＭＯＤＥＬＦＯＲＳＴＵＤＹＩＮＧＥＲＯＳＩＯＮＡＮＤＳＥＤＩＭＥＮＴＡＴＩＯＮＯＦＡＷＡＮＤＥＲＩＮＧＳＴＲＥＴＣＨＩＮＴＨＥＬＯＷＥＲＹＥＬＬＯＷＲＩＶＥＲ￥ＺＨＡＮＧＳｈｉｑｉｒ（Ｓｅｎｉｏｒ...     

NUMERICAL STUDY ON THE MIXING OF UNSORTED SEDIMENT PARTICLES DISCHARGED INTO A CROSS-FLOW BY MULTIPHASE PARTICLE-IN-CELL （MP-PIC） METHOD

The mixing characteristics of dredged sediments of variable size discharged into cross-flow are studied by an Eulerian-Lagrangian method. A three-dimensional (3D) numerical model has been developed by using the modified k-ε parameterization for the turbulence in fluid phase/water and a Lagrangian method for the solid phase/sediments. In the model the wake turbulence induced by sediments has been included as additional source and sink terms in the k-ε model; and the trajectories of the sediments are tracked by the Lagrangian method in which the sediment drift velocities in cross-flow are computed by a multiphase particle-in-cell (MP-PIC) method and the diffusion process is approximated by a random walk model. The hydrodynamic behavior of dumped sediment cloud is governed by the total buoyancy on the cloud, the drag force on each particle and velocity of cross-flow. The cross-flow destroys more or less the double vortices occurred in stagnant ambience and dominates the longitudinal movement of sediment cloud. The computed results suggest satisfactory agreement by comparison with the experimental results of laboratory.