Numerical modeling of nonlinear water waves over heterogeneous porous beds   总被引：1，自引：0，他引：1  

Eric C. Cruz  Qin Chen   《Ocean Engineering》2007,34(8-9):1303-1321

The transformation of nonlinear water waves over porous beds is studied by applying a numerical model based on Chen's [2006. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. Journal of Engineering Mechanics, 132:2, 220–230] Boussinesq-type equations for highly nonlinear waves on permeable beds. The numerical model uses a high-order time-marching solution and fourth-order finite-difference schemes for discretization of first-order spatial derivatives to obtain a computational accuracy consistent with the model equations. By forcing the wave celerity and spatial porous-damping rate of the linearized model to match the exact linear theory for horizontal porous bed over a prescribed range of relative depths, the values of the model parameters are optimally determined. Numerical simulations of the damped wave propagation over finite-thickness porous layer demonstrate the accuracy of both the numerical model and governing equations, which have been shown by prior theoretical analyses to be accurate for both nominal and thick porous layers. These simulations also elucidate on the significance of the higher-order porous-damping terms and the influence of the hydraulic parameters. Application of the model to the simulation of the wave field around a laboratory-scale submerged porous mound provides a measure of its capability, as well as useful insight into the scaling of the porous-resistance coefficients. For application to heterogeneous porous beds, the assumption of weak spatial variation of the porous resistance is examined using truncated forms of the governing equations. The results indicate that the complete set of Boussinesq-type equations is applicable to porous beds of nonhomogeneous makeup.  相似文献   

Integration of local–global upscaling and grid adaptivity for simulation of subsurface flow in heterogeneous formations     

M. Gerritsen  J. V. Lambers 《Computational Geosciences》2008,12(2):193-208

We propose a methodology, called multilevel local–global (MLLG) upscaling, for generating accurate upscaled models of permeabilities or transmissibilities for flow simulation on adapted grids in heterogeneous subsurface formations. The method generates an initial adapted grid based on the given fine-scale reservoir heterogeneity and potential flow paths. It then applies local–global (LG) upscaling for permeability or transmissibility [7], along with adaptivity, in an iterative manner. In each iteration of MLLG, the grid can be adapted where needed to reduce flow solver and upscaling errors. The adaptivity is controlled with a flow-based indicator. The iterative process is continued until consistency between the global solve on the adapted grid and the local solves is obtained. While each application of LG upscaling is also an iterative process, this inner iteration generally takes only one or two iterations to converge. Furthermore, the number of outer iterations is bounded above, and hence, the computational costs of this approach are low. We design a new flow-based weighting of transmissibility values in LG upscaling that significantly improves the accuracy of LG and MLLG over traditional local transmissibility calculations. For highly heterogeneous (e.g., channelized) systems, the integration of grid adaptivity and LG upscaling is shown to consistently provide more accurate coarse-scale models for global flow, relative to reference fine-scale results, than do existing upscaling techniques applied to uniform grids of similar densities. Another attractive property of the integration of upscaling and adaptivity is that process dependency is strongly reduced, that is, the approach computes accurate global flow results also for flows driven by boundary conditions different from the generic boundary conditions used to compute the upscaled parameters. The method is demonstrated on Cartesian cell-based anisotropic refinement (CCAR) grids, but it can be applied to other adaptation strategies for structured grids and extended to unstructured grids.  相似文献   

Accurate local upscaling with variable compact multipoint transmissibility calculations     

James V. Lambers  Margot G. Gerritsen  Bradley T. Mallison 《Computational Geosciences》2008,12(3):399-416

We propose a new single-phase local upscaling method that uses spatially varying multipoint transmissibility calculations. The method is demonstrated on two-dimensional Cartesian and adaptive Cartesian grids. For each cell face in the coarse upscaled grid, we create a local fine grid region surrounding the face on which we solve two generic local flow problems. The multipoint stencils used to calculate the fluxes across coarse grid cell faces involve the six neighboring pressure values. They are required to honor the two generic flow problems. The remaining degrees of freedom are used to maximize compactness and to ensure that the flux approximation is as close as possible to being two-point. The resulting multipoint flux approximations are spatially varying (a subset of the six neighbors is adaptively chosen) and reduce to two-point expressions in cases without full-tensor anisotropy. Numerical tests show that the method significantly improves upscaling accuracy as compared to commonly used local methods and also compares favorably with a local–global upscaling method.  相似文献   

Using Sequential Self-Calibration Method to Identify Conductivity Distribution: Conditioning on Tracer Test data     

Bill X. Hu  Changming He 《Mathematical Geosciences》2008,40(8):845-859

An iterative inverse method, the sequential self-calibration method, is developed for mapping spatial distribution of a hydraulic conductivity field by conditioning on nonreactive tracer breakthrough curves. A streamline-based, semi-analytical simulator is adopted to simulate solute transport in a heterogeneous aquifer. The simulation is used as the forward modeling step. In this study, the hydraulic conductivity is assumed to be a deterministic or random variable. Within the framework of the streamline-based simulator, the efficient semi-analytical method is used to calculate sensitivity coefficients of the solute concentration with respect to the hydraulic conductivity variation. The calculated sensitivities account for spatial correlations between the solute concentration and parameters. The performance of the inverse method is assessed by two synthetic tracer tests conducted in an aquifer with a distinct spatial pattern of heterogeneity. The study results indicate that the developed iterative inverse method is able to identify and reproduce the large-scale heterogeneity pattern of the aquifer given appropriate observation wells in these synthetic cases.  相似文献   

Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area     

P. Tullen  P. Turberg  A. Parriaux 《Engineering Geology》2006,87(3-4):195-204

Three methods were combined to determine the groundwater recharge and transfer processes of a landslide prone area. First, the radiomagnetotelluric method was used to investigate the distribution of electrical resistivity (ρ) of the subsurface and build a three-dimensional model of permeability (k), through an experimental relation between ρ and k. Second, this structural model of permeability and additional climatologic data were used to fix boundary and recharge conditions to perform a three-dimensional and transient numerical simulation of the groundwater flow. Finally 18-Oxygen time series observed at the main springs were used to validate the model. This association of methods led to an improved characterization of the groundwater flow system at local scale and a better understanding of the role of this system on the landslide phenomenon. This structured approach is thought to be useful to design specific remediation strategies to drain the unstable mass.  相似文献   

Determination of Area-Averaged Water Vapour Fluxes with Large Aperture and Radio Wave Scintillometers over a Heterogeneous Surface – Flevoland Field Experiment     

W. M. L. Meijninger  A. E. Green  O. K. Hartogensis  W. Kohsiek  J. C. B. Hoedjes  R. M. Zuurbier  H. A. R. De Bruin 《Boundary-Layer Meteorology》2002,105(1):63-83

A large aperture scintillometer (LAS) andradio wave scintillometer (RWS)were installed over a heterogeneous areato test the applicability of the scintillation method.The heterogeneity in the area, whichconsisted of many plots, was mainly caused bydifferences in thermal properties ofthe crops; the variations in theaerodynamic roughness lengthwere small. The water vapour fluxesderived from the combined LAS-RWSsystem, also known as the two-wavelengthmethod, agreed fairly well with the aggregatedwater vapour fluxes derived from in-situeddy covariance measurements. The water vapourfluxes derived from a stand-alone LASare also presented. It was found that a single LASand an estimate of the area averagedavailable energy (using a simple parameterisationscheme) can provide also reasonablearea-averaged water vapour fluxes.  相似文献   

Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data   总被引：1，自引：0，他引：1  

Harrie-Jan Hendricks Franssen  Jaime Gómez-Hernández  Andrés Sahuquillo 《Journal of Hydrology》2003,281(4):281-295

This paper presents the extension of the self-calibrating method to the coupled inverse modelling of groundwater flow and mass transport. The method generates equally likely solutions to the inverse problem that display the variability as observed in the field and are not affected by a linearisation of the state equations. Conditioning to the state variables is measured by an objective function including, among others, the mismatch between the simulated and measured concentrations. Conditioning is achieved by minimising the objective function by gradient-based methods. The gradient contains the partial derivatives of the objective function with respect to: log conductivities, log storativities, prescribed heads at boundaries, retardation coefficients and mass sources. The derivatives of the objective function with respect to log conductivity are the most cumbersome and need the most CPU-time to be evaluated. For this reason, to compute this derivative only advective transport is considered. The gradient is calculated by the adjoint-state method. The method is demonstrated in a controlled, synthetic study, in which the worth of concentration data is analysed. It is shown that concentration data are essential to improve transport predictions and also help to improve aquifer characterisation and flow predictions, especially in the upstream part of the aquifer, even in the case that a considerable amount of other experimental data like conductivities and heads are available. Besides, conditioning to concentration data reduces the ensemble variances of estimated transmissivity, hydraulic head and concentration.  相似文献   

地表非均匀性对区域平均水分通量参数化的影响   总被引：3，自引：3，他引：3  

刘晶淼  丁裕国  周秀骥  汪方 《气象学报》2003,61(6):712-717

次网格尺度地表非均匀性对于网格区平均通量具有重要影响。若将网格区视为均一地表 ,并不能真实描述地 气通量交换过程 ,且可造成很大误差。文中从理论上证明 ,区域平均水分通量的变化率可分解为两部分 :第一部分为区域水分通量的算术平均变化率 ;第二部分为非均匀性所引起的水分通量变化率扰动 ,它与区域内土壤水分空间分布的变差系数有关。数值试验表明 ,地表土壤水分的水平空间变差系数集中反映了区域内土壤水分分布的非均匀程度 ,不同土壤对同样的非均匀程度其敏感性是不同的。变差系数愈大 ,非均匀性愈强 ,在相同的土壤水分平均值下 ,不同土壤类型对地表非均匀程度的敏感性并不相同。例如沙土和粘土受非均匀性的影响就可相差数十倍。  相似文献   

Borehole and surface-based hydrogeophysics     

Roger Guérin 《Hydrogeology Journal》2005,13(1):251-254

  相似文献   

Pore structure and heterogeneity of shale gas reservoirs and its effect on gas storage capacity in the Qiongzhusi Formation     

Shangbin Chen  Zhuo Gong  Xueyuan Li  Huijun Wang  Yang Wang  Yingkun Zhang 《地学前缘(英文版)》2021,12(6):192-208

Fine characterization of pore systems and heterogeneity of shale reservoirs are significant contents of shale gas reservoir physical property research.The research on micro-control factors of low productivity in the Qiongzhusi Formation(Fm.)is still controversial.The lower Cambrian Qiongzhusi Fm.in the Qujing,Yunnan was taken as the object to investigate the influence of mineral compositions on the phys-ical properties of the reservoir and the heterogeneity of shale,using the algorithm to improve the char-acterization ability of Atomic Force Microscopy(AFM).The results showed that:(1)The pores are mainly wedge-shaped pores and V-shaped pores.The pore diameter of the main pore segment ranges from 5 to 10 nm.Mesopores are mainly developed in the Qiongzhusi Fm.shale in Well QD1,with the average pore diameter of 6.08 nm.(2)Microscopic pore structure and shale surface properties show strong hetero-geneity,which complicates the micro-migration of shale gas and increases the difficulty of identifying high-quality reservoirs.(3)The increase of clay mineral content intensifies the compaction and then destroys the pores.Conversely,brittle minerals can protect pores.The support and protection of brittle minerals to pores space depend on their content,mechanical properties and diagenesis.(4)Compression damage to pores,large microscopic roughness and surface fluctuations and strong pore structure heterogeneity are the reasons for the poor gas storage capacity of the Qiongzhusi Fm.,which will lead to poor productivity in the Qiongzhusi Fm.  相似文献