Mathematical simulation on the oil slick spreading and dispersion in nonuniform flow fields

The viscous fluid boundary layer equations were adopted to characterize the velocity distribution across the vertical section of the oil slick on moving water. The velocity profile was found to be the combination of linear and parabolic distribution. A numerical model including spreading and dispersion was developed to describe the oil slick’s early movement in the open and ice-covered water. The flume test was conducted to determine the dispersion coefficients (K_x) and the effects of velocity and wave height on the slick’s dispersion were also investigated. In the open water, K_x increased from 4.34 cm²/s to 20.08 cm²/s as the velocity changed from 3 cm/s to 12 cm/s. Acoefficient β that characterized K_x fluctuated at 1.5 when wave heights were between 20 mm and 70 mm. Under ice, the slick didn’t move until the velocity exceeded 6 cm/s and K_x increased from 2.69 cm²/s to 5.64 cm²/s as the velocity changed from 6 cm/s to 12 cm/s. β remained 0.4 when wave heights were between 20 mm and 60 mm. The model performed very well in predicting the slick’s position and length during the gravitation-inertia phase for the two situations when relevant dispersion coefficients were input.     

MT3D-通用的三维地下水污染物运移数值模型

定量研究污染物在地下水中的运移过程通常采用数值模拟方法。MT3D是一套基于有限差分方法的污染物运移模拟软件，近年来在国外水文地质和水环境模拟等领域的研究中已经得到较为广泛的认可。MT3D比较全面地考虑了污染物在地下水中的对流、弥散和化学反应等过程，可以灵活处理各种复杂的源汇项和边界条件，能够准确模拟承压、无压和越流含水层中的污染物运移过程。MT3D具有模块化的程序结构、灵活的求解方法以及全面的模拟功能，非常适合实际问题的研究，值得在国内推广使用。     

A Fully Coupled Chemo-Poroelastic Analysis of Pore Pressure and Stress Distribution around a Wellbore in Water Active Rocks

Water active rocks consist of minerals that hold water in their crystalline structure and in pore spaces. Free water from drilling fluid can be attracted by the formation depending on the potential differences between pore space and drilling fluid. The fluid movement into the formation or out of the formation can lead to a change in effective stress, thus causing wellbore failures. In all previous studies it is found that the solute transport from or to the formation is primarily controlled by diffusion process and the effect of advection on solute transfer is negligible for a range of very low permeable shale formations (>10⁻⁵ mD). In this study a range of permeable shale formations (10⁻⁵ to 10⁻³ mD) commonly encountered in drilling oil and gas wells are considered to investigate the solute transfer between drilling fluid and formation due to advection. For this purpose a finite element model of fully coupled chemo-hydro-mechanical processes was developed. Results of this study revealed that the solute transfer between the drilling fluid and the shale formation is controlled primarily by permeability of the shale formations. For the range of shale formations studied here, there exists a threshold permeability below which the solute transfer is dominated by diffusion process and above which by fluid in motion (fluid flow). Results from the numerical experiments have shown that when the permeability of shales is greater than this threshold permeability, the chemical potential gradient between the pore fluid and drilling fluid reaches equilibrium faster than that when the permeability of shales is below this threshold value. Also it has been found that when advection is taken into account, effective radial and tangential stresses decrease around the wellbore, particularly near the wellbore wall where the solute concentration has reached near equilibrium.     

Contaminant transport with non‐equilibrium processes in unsaturated soils and implicit characteristic Galerkin scheme

A non‐equilibrium sorption—advection—diffusion model to simulate miscible pollutant transport in saturated–unsaturated soils is presented. The governing phenomena modelled in the present simulation are: convection, molecular diffusion, mechanical dispersion, sorption, immobile water effect and degradation, including both physical and chemical non‐equilibrium processes. A finite element procedure, based on the characteristic Galerkin method with an implicit algorithm is developed to numerically solve the model equations. The implicit algorithm is formulated by means of a combination of both the precise and the traditional numerical integration procedures. The stability analysis of the algorithm shows that the unconditional stability of the present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results illustrate good performance of the present algorithm in stability and accuracy, and in simulating the effects of all the mentioned phenomena governing the contaminant transport and the concentration distribution. Copyright © 2000 John Wiley & Sons, Ltd.     

The effect of viscosity ratio on the dispersal of fracturing fluids into groundwater system

In the development of oil and gas reservoirs, the transport of miscible fluids in porous rocks is a key issue for oil and gas recovery. The simplified unidirectional flow model is employed to investigate the effects of the viscosity ratio on dispersion in semi-infinite homogenous media. The viscosity is supposed to be unsteady due to changing component concentration over time. In cases of both a viscosity ratio larger and less than 1, the pollutant concentration and flow velocity are computed at different initial conditions and viscosity ratios. The analytical solutions are then developed by introducing new variables and transforming the equation governing advection–diffusion equation in semi-infinite homogeneous media with a continuous source. A comparison of the numerical solution with the analytical solution revealed a similarity over 98%, highlighting the usability of the analytical solution. If the viscosity ratio is larger than 1, flow velocity declines exponentially and concentration attenuates with transporting time. In addition, the timescale plays a significant role and the effects become more prominent in long-term transport. In the case of a viscosity ratio less than 1, both the timescale and viscosity ratio variables have little influence on the changing speed of the concentration profile. This work helps to predict the position and the time required to reach the harmless pollutant concentration when monitoring fracturing fluids transportation into groundwater system and would be especially useful in designing and interpreting laboratory experiments studying the miscible flow.     

A discussion on analytical and numerical modelling of the land subsidence induced by coal seam gas extraction

Coal seam gas (CSG) is an increasingly important source of natural gas all over the world. Although the influence of conventional oil and gas extraction on surface subsidence has been widely recognized and studied, few studies are carried out on the surface subsidence in coal seam gas fields and its impact on surface infrastructure and the environment. This paper discusses modelling of the surface subsidence associated with coal seam gas production by applying both analytical and numerical methods. By comparison of results from the numerical model and two analytical models, i.e. the disc-shaped reservoir model and the uniaxial compaction model, the analytical solutions cannot describe the complex process of water and gas extraction and have the limitations to predict the surface subsidence, while the numerical model can be better used in prediction of subsidence. After applying the numerical model in numerical analysis, the deformation characteristics of coupled fluid flow, and the effects of permeability change of coal seam, associated overlying and underlying layers, and depressurization rates on surface subsidence are investigated. The results demonstrate that the proposed model can simulate the production of water and gas from coal seams and the associated surface subsidence.     

含水层热量运移中自然热对流作用的数值模拟

较严密地推导了能够描述含水层热量运移中自然热对流作用的水流方程和热量运移方程,并在此基础上建立了一个三维非线性模型.根据模型的特点,提出综合采用Douglas Brian ADI差分法、特征线法返回技术及改进Picard迭代法求解,并给出了计算步骤.模型及数值方法用于模拟上海第二承压含水层中的一个群井储能试验.试验中的自然热对流作用得到了较好的体现,反映出模型合理可靠,方法有效,结果令人满意.     

波浪作用下沙坝剖面形成过程的数值模拟

建立了波浪、底部离岸流、泥沙运动和沙滩剖面演变耦合数学模型,模拟了实验室中波浪作用下沙坝剖面的形成过程,讨论了各个物理参数对剖面上的水动力和地形变化的影响。通过参数率定,较好地复演了不规则波作用下沙滩地形由均匀斜坡向沙坝剖面的演变,波高、底部离岸流、输沙率和剖面演变的计算结果与实测数据吻合良好。采用数值实验和误差分析方法,讨论了破碎波水滚倾角、泥沙扩散系数和床面休止角等物理参数对计算结果的影响。研究表明,所建立的模型能较好地描述波浪作用下沙坝剖面形成过程中的重要物理机制。     

Newspapers as a validation proxy for GIS modeling in Fujairah,United Arab Emirates: identifying flood-prone areas

Prediction and control of blast-induced ground vibration is a matter of concern in mining industry since long. Several approaches ranging from scaled distance regression, different numerical methods to wave superimposition theories have been tried by many researchers for better prediction and control of blast-induced ground vibration. Signature hole analysis is one of the popular simulation methods to predict the ground vibration generated due to production blast. It superimposes the recorded signature hole waveform using a computer program to predict the production blast-induced vibration. The technique inputs the designated time of detonation of each hole and superimposes the waves generated by each hole to predict the nearest value of peak particle velocity and frequency of blast-induced ground vibration. Although a very useful approach, it requires a computer program to simulate the linear superimposition of waveforms. The simulation is not possible for every blast as it takes time and also is difficult for field engineers to simulate every time, whereas it is always easy for blasting engineers to adapt and use an empirical equation/approach for prediction and control of blast-induced ground vibration than simulation. In this paper, an attempt has been made to develop an innovative and simplified analytical approach of signature hole analysis. The simplified sinusoidal wave equation is obtained from recorded signature hole ground vibration waveform properties and is superimposed mathematically according to the multi-hole blast design to predict the production blast-induced ground vibrations. The validation of the developed approach was done in three different sites, and up to 15% more accuracy in prediction of the blast, vibrations are achieved in comparison with signature hole analysis prediction.

     

基于GRAPES的西北地区沙尘暴数值预报模式及其应用研究

介绍了中国气象局兰州干旱气象研究所和中国气象科学研究院数值预报研究中心合作研制的耦合于GRAPES（Global/Regional Assimilation and Prediction Enhanced System）的沙尘暴数值模式GRAPES_SDM。该模式包括沙尘的起沙、传输、吸湿增长、并合、干沉降与云下清洗等详细的物理过程，可以对沙尘暴的起沙和空气中沙尘浓度进行模拟和预报。在此基础上形成了西北地区的沙尘暴数值预报模式系统，并于2005年4月开始在兰州中心气象台试运行，同时还将模式结果与卫星遥感资料反演的沙尘暴监测结果进行了对比验证。以2005年4月17~19日和5月28~29日发生在西北地区的2次强沙尘暴为例，利用GRAPES_SDM对这两次沙尘天气的起沙、传输、扩散直至消散进行了数值模拟。结果表明，模式输出的沙尘时空分布与实况观测和卫星云图监测的沙尘分布范围基本一致，说明该模式系统对西北地区沙尘暴天气的起沙、传输有较好的模拟和预报能力。     

Detecting surface oil slick related to gas hydrate/petroleum on the ocean bed of South China Sea by ENVI/ASAR radar data

For exploring the distribution of oil slick related to gas hydrate or petroleum in the northern part of South China Sea and evaluating the potential of ENVI/ASAR radar data for detecting oil slick, this paper tries to detect the oil slick on ocean surface of some potential areas in South China Sea and map oil slick from anomaly area in radar image. Seven surface oil slicks in Qiongdongnan Basin and four potential targets of oil slicks in Taixi–Dongsha Basin were detected and extracted by using ENIV/ASAR data. The anomalies of the radar image representing oil slick are validated by in-site sampling and laboratory analyzing. The Gas Chromatograph (GC) results of the extracted oil from the water samples gathered from the sites of detected oil slick demonstrate a typical composition of oil rather than that of ocean microbial origins, and Total Organic Carbon (TOC) concentrations in two designed vertical profiles keep steady in different depths, showing that the organic contents in water column are steady from the deep area to ocean surface. These results suggest the detected oil slick is a thin floating oil rather than a thick spill. Influencing factors to the results of oil slick detection including wind speed, current, small islands and capes as well as deep-sea petroleum reservoirs are discussed. The results in this study exhibit that extracting the anomaly signals related to oil slick related to oil/gas seeps using image processing techniques from ENVI/ASAR data is helpful to predict anomaly areas as potential targets for further exploration in this area.     

感潮河段丁坝局部冲刷三维数值模拟

针对局部冲刷和一般冲刷的不同,建立了考虑垂向水流作用的局部冲刷三维紊流泥沙数学模型,并对往复流和单向流作用下的丁坝局部冲刷进行了验证计算和数值模拟。结果表明,往复流作用下丁坝局部冲刷在冲刷过程和冲淤形态上与单向流有明显不同,若采用单向流作用下的局部冲刷模型试验结果来指导往复流作用下丁坝头防护工程的设计和施工是不利的。     

The motion of sediment-water mixtures during intense bedload transport: computer simulations

A new model, which couples fluid and particle dynamics, has been developed to study the motion of the sediment-water mixture during intense bedload transport, including the velocity profiles of both sediment and water, the roughness length of an upper plane bed and the thickness of moving sediment layers. Standard mixing length theory is used to model the motion of water above the boundary between the overlying water and the sediment-water mixture. The turbulent flow within the moving sediment layers is described by a shear stress model, in which the effective viscosity of the flowing water is proportional to the velocity difference between the fluid and the sediment. The particle dynamics method, in which the equations of motion of each of many particles are solved directly, is applied to model the movement of sediment particles. The particle-fluid interaction is expressed by a velocity-squared fluid drag force exerted on each sediment particle. Both computer simulation results and theoretical analysis have shown that the velocities of both sediment and fluid during intense sediment transport decrease exponentially with depth in the top layers of a fast-moving sediment—water mixture. The thickness of the moving sediment layers, obtained from the computer simulation results, is proportional to the shear stress, which agrees with previous experimental observations.     

A particle tracking transport method for the simulation of resident and flux-averaged concentration of solute plumes in groundwater models

A new numerical technique called the convolution-based particle tracking (CBPT) method is developed to simulate resident or flux-averaged solute concentrations in groundwater models. The method is valid for steady-state flow and linear transport processes such as sorption with a linear sorption isotherm and first-order decay. The CBPT method uses particle tracking to take advantage of the ability of particle-based approaches to maintain sharp fronts for advection-dominated transport problems common in groundwater modeling and because of the flexibility of the random walk method to simulate a wide range of possible forms of the dispersion tensor. Furthermore, the algorithm for carrying out the convolution and superposition calculation from particle tracking results is very efficient. We show that from a single particle tracking run, source term variability, sorption, and decay can all be simulated rapidly without rerunning the underlying transport model unless the flow field or dispersion parameters are changed. A series of verification simulations are presented to demonstrate the accuracy and efficiency of the CBPT method compared to more conventional particle tracking approaches.     

轻非水相液体在地下环境中的运移特征与模拟预测研究

通过自行设计制备的二维砂箱模拟轻非水相液体(LNAPL)在湿润多孔介质中的渗漏,得出LNAPL在地下环境中的入渗、迁移及分布规律,并利用多相流体理论对其污染机理进行分析.在实验基础上建立模型来预测渗漏的基本特征,即渗流带中污染锋面扩展速度和透镜体的最终形状,预测结果与实测结果较为接近.这种预测模型对控制和治理地下水中LNAPL污染是至关重要的.     

Numerical Modeling in the Structural Study of Technogenic Rock Array

The relevance of research underlies in the fact that numerical models of technogenic rock arrays can predict and evaluate structural and mechanical properties, which are very difficult or even impossible to study instrumentally. It is possible to simulate the technogenic rock array’s structure by the method of discrete elements, where each particle in the array structure is described separately. The article contains the method of numerical simulation of structure of bulk disjointed rock array, describes the basic initial and boundary conditions of the model; represents the visualization of numerical model of the array created with peripheral dumping (as the array with more textured and visual segregated granulometric composition). The ways of application of technogenic rock array numerical modeling in the mining process are described. Comparative analysis of real and simulated array structure shows the adequacy of the numerical model and the algorithm developed for its creation. In addition, the adequacy of the model is supported by a high convergence of the experimental and calculated dependence of the filtration rate values upon the particle size distribution of technogenic rock arrays.     

淤泥质河口水沙运动数学模型相关问题

为研究淤泥质河口的水沙运动规律,建立了用于模拟淤泥质河口水沙运动的二维数学模型。该模型采用基于无结构三角网格下的有限体积法对方程组进行离散,结合Roe-MUSCL方法及时间方向的预测-校正格式,使模型在时空方向具有二阶计算精度。模型中分别采用不同方法计算粘性和非粘性泥沙的输移源项,并引入粘性泥沙的起动流速和冲刷率计算公式。采用已有的概化水槽试验数据对模型进行了初步验证。然后模拟了1995年10月小潮及大潮期间海河口的潮流运动与泥沙输移过程,计算得到的潮位、潮流速及含沙量过程与实测过程符合较好,结果表明模型能够用来模拟淤泥质河口粘性和非粘性泥沙的不平衡输移过程。同时还比较了泥沙输移源项的不同处理方式对计算结果的影响,计算表明在淤泥质河口水沙运动数学模型中必须同时考虑粘性和非粘性泥沙的输移。     

NUMERICAL GEODYNAMIC MODELLING OF COUPLED MECHANO-THERMO-HYDROLOGICAL PROCESSES AND ITS APPLICATION IN PREDICTIVE EXPLORATION IN THE FENGHUANGSHAN ORE FIELD, TONGLING

Mineralisation is the result of the coupled multi-geodynamic processes in the crust. The coupled mechano-thermo-hydrological （MTH） processes are the basic physical processes that govern the location of the hydrothermal mineralization, which can be simulated in the computer by using of the numerical codes, such as FLAC. The numerical modeling results can be used not only to explain the features of existing ore deposits, but also to predict the fhvorable mineralization locations. This paper has summarized the basic equations describing coupled MHT processes in the water-saturated porous rocks, the principles of FLAC, and its application to the MHT processes related to copper mineralization in the Fenghuangshan ore field. We used the FLAC to simulate the syn-deformation cooling and fluid flowing evolution after the intrusion was emplaced and solidified. The modeling results suggest a most prospective exploration area where the subsequent exploration supported the prediction and the test bore hole disclosed the high quality copper ore bodies in the target, demonstrating a positive role of the numerical MTH modeling in facilitating predictive ore discovery.     

Using global node-based velocity in random walk particle tracking in variably saturated porous media: application to contaminant leaching from road constructions

Precise and efficient numerical simulation of transport processes in subsurface systems is a prerequisite for many site investigation or remediation studies. Random walk particle tracking (RWPT) methods have been introduced in the past to overcome numerical difficulties when simulating propagation processes in porous media such as advection-dominated mass transport. Crucial for the precision of RWPT methods is the accuracy of the numerically calculated ground water velocity field. In this paper, a global node-based method for velocity calculation is used, which was originally proposed by Yeh (Water Resour Res 7:1216–1225, 1981). This method is improved in three ways: (1) extension to unstructured grids, (2) significant enhancement of computational efficiency, and (3) extension to saturated (groundwater) as well as unsaturated systems (soil water). The novel RWPT method is tested with numerical benchmark examples from the literature and used in two field scale applications of contaminant transport in saturated and unsaturated ground water. To evaluate advective transport of the model, the accuracy of the velocity field is demonstrated by comparing several published results of particle pathlines or streamlines. Given the chosen test problem, the global node-based velocity estimation is found to be as accurate as the CK method (Cordes and Kinzelbach in Water Resour Res 28(11):2903–2911, 1992) but less accurate than the mixed or mixed-hybrid finite element methods for flow in highly heterogeneous media. To evaluate advective–diffusive transport, a transport problem studied by Hassan and Mohamed (J Hydrol 275(3–4):242–260, 2003) is investigated here and evaluated using different numbers of particles. The results indicate that the number of particles required for the given problem is decreased using the proposed method by about two orders of magnitude without losing accuracy of the concentration contours as compared to the published numbers.     

各向异性密度流模型及其验证

针对密度流流动各向异性的特性,采用非线性k-ε模式求取雷诺应力,结合RNG模型及浮力修正建立了三维密度流模型.应用该三维模型计算了水库温差异重流,分析了温差异重流演进过程,并就流速分布与水温下泄过程与实测资料进行了比较.结果表明,所建模型由于引入非线性RNG模式,能有效反映密度流的各向异性特性,较好地模拟了流场和温度场强烈耦合的水库温差异重流,能准确预测温度分层过程与异重流潜行时间,具有较高精度.