Probabilistic solution of floodplain inundation equation

Uncertainty in bed roughness is a dominant factor in providing a sufficiently accurate simulation of floodplain flows. This study describes a method to compute the transition probability density distribution of time-varying water elevations where the evolutionary process is based on a conventional one-dimensional storage cell model with governing stochastic differential equation. By including the random inputs (or noise terms) of bed roughness and initial water depth, time-dependent and spatially varying probability density function of the water surface leads to a Fokker–Planck equation. The model’s performance is evaluated by applying it to shallow water flow with a horizontal bed. Sensitivity of model predictions to variations in the bed friction parameters is shown. By comparing the result of the proposed method with that of conventional Monte Carlo simulation, the advantage of the former as a method for density function prediction is confirmed.     

用共方位角偏移实现叠前二维资料三维化成像

波动方程偏移保持了波场动力学特征,依此本文应用共方位角叠前深度偏移技术来实现二维资料三维化处理.本方法能使共反射点偏移成像与叠前数据插值一步完成,较好地解决二维资料处理中的成像点不能准确归位、侧面反射波难处理、主测线与联络线不闭合等问题.     

Comparing the performance of a 2‐D finite element and a 2‐D finite volume model of floodplain inundation using airborne SAR imagery

The performances of a finite volume model (SFV) and finite element model (TELEMAC‐2D) in reproducing inundation on a 16 km reach of the river Severn, United Kingdom, are compared. Predicted inundation extents are compared with 4 airborne synthetic aperture radar images of a major flood event in November 2000, and these are used to calibrate 2 values of Manning's n for the channel and floodplain. The four images are shown to have different capacities to constrain roughness parameters, with the image acquired at low flow rate doing better in determining these parameters than the image acquired at approximately peak flow. This is assigned to the valley filling nature of the flood and the associated insensitivity of flood extent to changes in water level. The level of skill demonstrated by the models, when compared with inundation derived using a horizontal water free surface, also increases as flow rate drops. The two models show markedly different behaviours to the calibration process, with TELEMAC showing less sensitivity and lower optimum values for Manning's n than SFV. When the models are used in predictive mode, calibrated against one image and predicting another, SFV performs better than TELEMAC. Copyright © 2007 John Wiley & Sons, Ltd.     

高分辨有限体积法模拟二维声波传播

论文使用高分辨有限体积法模拟二维声波在多层介质中的传播现象.该方法需要求解双曲型偏微分方程的黎曼解,数值解在强速度变化的分界面上满足数值通量守恒,可以有效地计算数值解的不连续性.根据数值实验给出的波场快照和接收点上的波形变化,可以清晰地观察到声波在介质分界面上的多次反射波,说明该方法不但可以得到高分辨率的波场快照,而且...     

二维SH波方程的半解析解及其数值模拟

本文以波动理论为基础, 半解析化求解地震勘探中常用的SH波方程. 获得的主要结果包括: 给出了二维均匀介质中SH波方程的解析解; 利用Cagniard-de Hoop方法详细推导了二维双层介质中SH波方程的解析解, 获得了透射波的解析解表达式. 同时, 基于SH波方程的解析表达式, 给出了包含各种波(如直达波、反射波、首波以及透射波)的解析解和波形图. 对于比较复杂的积分型解析解, 利用数值积分方法给出了数值结果, 并与优化的近似解析离散化方法(ONADM)和4阶Lax-Wendroff修正方法(LWC)的数值结果进行了比较, 以验证解析解的正确性. 本文的研究成果有望在检验波动方程数值新方法的有效性、波传播理论分析等方面得到应用.     

基于第二代小波自适应网格的二维声波方程波传模拟

本文利用第二代小波多尺度分解和快速变换的特点,构造自适应计算网格.对初始计算网格上的数值解进行第二代小波变换,得到数值解对应的小波系数空间.小波系数的大小表示相邻网格上数值变化率,小波系数大的区域网格点上的数值解变化梯度大.当小波系数大于等于预设的阈值时,在小波系数对应的网格点周围插入新的计算网格点,通过阈值可以实现网格的细化,得到多尺度下层层嵌套的细化自适应网格;由有限差分法得到相应网格点的空间导数.比较数值算例得到的波场快照和计算时间,验证了该方法的有效性.     

2.5D finite-difference solution of the acoustic wave equation

The finite‐difference method applied to the full 3D wave equation is a rather time‐consuming process. However, in the 2.5D case, we can take advantage of the medium symmetry. By taking the Fourier transform with respect to the out‐of‐plane direction (the symmetry axis), the 3D problem can be reduced to a repeated 2D problem. The third dimension is taken into account by a sum over the corresponding wave‐vector component. A criterion for where to end this theoretically infinite sum derives from the stability conditions of the finite‐difference schemes employed. In this way, the computation time of the finite‐difference calculations can be considerably reduced. The quality of the modelling results obtained with this 2.5D finite‐difference scheme is comparable to that obtained using a standard 3D finite‐difference scheme.     

基于静主元消元法的频率域波动方程正演(英文)

频率域波动方程正演是求解一个大型线性稀疏方程组问题,其受到计算效率和内存存储问题的限制。常规的高斯消元法不能满足大型数据的并行计算,本文提出基于静主元消元法(GESP)进行稀疏矩阵LU分解和多炮有限差分正演,该方法不仅提高了稳定性,更有利于单频点内LU分解的分布式并行计算。通过Marmousi模型模拟试验,单频波场和转化到时间域地震剖面的试验表明模拟精度和计算效率得到提高,节约并充分利用内存,为波形反演奠定基础。     

Quantifying peakflow attenuation/amplification in a karst river using the diffusive wave model with lateral flow

The aim of this paper is to quantify peakflow attenuation and/or amplification in a river, investigating lateral flow from the intermediate catchment during floods. This is a challenge for the study of the hydrological response of permeable/intermittent streams, and our contribution refers to a modelling framework based on the inverse problem for the diffusive wave model applied in a karst catchment. Knowing the upstream and downstream hydrographs on a reach between two stations, we can model the lateral one, given information on the hydrological processes involved in the intermediate catchment. The model is applied to 33 flood events in the karst reach of the Iton River in French Normandy where peakflow attenuation is observed. The monitored zone consists of a succession of losing and gaining reaches controlled by strong surface‐water/groundwater (SW/GW) interactions. Our results show that despite a high baseflow increase in the reach, peakflow is attenuated. Model application shows that the intensity of lateral outflow for the flood component is linked to upstream discharge. A combination of river loss and overbank flow for highest floods is proposed for explaining the relationships. Our approach differentiates the role of outflow (river loss and overbank flow) and that of wave diffusion on peakflow attenuation. Based on several sets of model parameterization, diffusion is the main attenuation process for most cases, despite high river losses of up to several m³/s (half of peakflow for some parameterization strategies). Finally, this framework gives new insight into the SW/GW interactions during floods in karst basins, and more globally in basins characterized by disconnected river‐aquifer systems.     

Solution of the nonlinear transport equation using modified Picard iteration

The transport and fate of reactive chemicals in groundwater is governed by equations which are often difficult to solve due to the nonlinear relationship between the solute concentrations for the liquid and solid phases. The nonlinearity may cause mass balance errors during the numerical simulation in addition to numerical errors for linear transport system. We have generalized the modified Picard iteration algorithm of Celia et al.⁵ for unsaturated flow to solve the nonlinear transport equation. Written in a ‘mixed-form’ formulation, the total solute concentration is expanded in a Taylor series with respect to the solution concentration to linearize the transport equation, which is then solved with a conventional finite element method. Numerical results of this mixed-form algorithm are compared with those obtained with the concentration-based scheme using conventional Picard iteration. In general, the new solver resulted in negligible mass balance errors (< ∥10⁻⁸∥%) and required less computational time than the conventional iteration scheme for the test examples, including transport involving highly nonlinear adsorption under steady-state as well as transient flow conditions. In contrast, mass balance errors resulting from the conventional Picard iteration method were higher than 10% for some highly nonlinear problems. Application of the modified Picard iteration scheme to solve the nonlinear transport equation may greatly reduce the mass balance errors and increase computational efficiency.     

MAST-2D diffusive model for flood prediction on domains with triangular Delaunay unstructured meshes

A new methodology for the solution of the 2D diffusive shallow water equations over Delaunay unstructured triangular meshes is presented. Before developing the new algorithm, the following question is addressed: it is worth developing and using a simplified shallow water model, when well established algorithms for the solution of the complete one do exist?The governing Partial Differential Equations are discretized using a procedure similar to the linear conforming Finite Element Galerkin scheme, with a different flux formulation and a special flux treatment that requires Delaunay triangulation but entire solution monotonicity. A simple mesh adjustment is suggested, that attains the Delaunay condition for all the triangle sides without changing the original nodes location and also maintains the internal boundaries. The original governing system is solved applying a fractional time step procedure, that solves consecutively a convective prediction system and a diffusive correction system. The non linear components of the problem are concentrated in the prediction step, while the correction step leads to the solution of a linear system of the order of the number of computational cells. A semi-analytical procedure is applied for the solution of the prediction step. The discretized formulation of the governing equations allows to handle also wetting and drying processes without any additional specific treatment. Local energy dissipations, mainly the effect of vertical walls and hydraulic jumps, can be easily included in the model.Several numerical experiments have been carried out in order to test (1) the stability of the proposed model with regard to the size of the Courant number and to the mesh irregularity, (2) its computational performance, (3) the convergence order by means of mesh refinement. The model results are also compared with the results obtained by a fully dynamic model. Finally, the application to a real field case with a Venturi channel is presented.     

Sub-grid scale parameterisation of 2D hydrodynamic models of inundation in the urban area

This paper presents preliminary results from a study considering the parameterisation of coarse-grid 2D flood models to take into account sub-grid scale flow patterns occurring in the urban area. A simulation of a severe flood in an urbanized coastal floodplain is first run using a fine grid that can resolve the flow around and between buildings. Next, the same model is run again using the same underlying topography, although stripped from any buildings, and a set of 7 values of the roughness parameter (Manning’s n), all larger than (or equal to) the value used in the original run. A further set of simulations is carried out using a five-fold increased grid cell size. It is found that while it may be possible to model the overall effects of the buildings using strongly increased roughness parameter values, using a coarse grid otherwise has implications related to the loss of information about the site topography that results in flood flow routes being incorrectly modelled.     

Numerical simulation of 2D wave propagation on unstructured grids using explicit differential operators

We present a numerical method of simulating seismic wave propagation on unstructured 2D grids. The algorithm is based on the velocity–stress formulation of the elastic wave equation and therefore uses a staggered grid approach. Unlike finite-element or spectral-element methods, which can also handle flexible unstructured grids, we use explicit differential operators for the calculation of spatial derivatives in each time step. As shown in previous work, three types of these operators are used, and their particular performance is analysed and compared with standard explicit finite-difference operators on regular quadratic and hexagonal grids. Our investigations are especially focused on the influence of grid irregularity, sampling rate (i.e. gridpoints per wavelength) and numerical anisotropy on the accuracy of numerical seismograms. The results obtained from the various methods are therefore compared with analytical solutions. The algorithm is then applied to a number of models that are difficult to handle using (quasi-)regular grid methods. Such alternative techniques may be useful in modelling the full wavefield of bodies with complex geometries (e.g. cylindrical bore-hole samples, 2D earth models) and, because of their local character, they are well suited for parallelization.     

一种优化的频率域三维声波有限差分模拟方法

为提高频率域有限差分（FD,finite-difference）正演模拟技术的计算精度和效率,基于旋转坐标系统的优化差分格式被广泛应用,但是只应用于正方形网格的情况.基于平均导数法（ADM）的优化差分格式,应用于正方形和长方形网格模拟.这些频率域有限差分算子,各自具有不同的差分格式和对应的优化系数求解表达式.本文基于三维声波方程发展了一种新的优化方法,只要给定FD模板形式,可直接构造频散方程,求取FD模板上各节点的优化系数.此方法的优点在于频率域FD算子的优化系数对应各个节点,可扩展优化其他格式.运用此优化方法,计算得到了不同空间采样间距比情况下27点和7点格式的优化系数.数值实验表明,优化27点格式与ADM 27点格式具有相同的精度,优化7点格式比经典的7点格式具有更小的数值频散.

     

对地震图像进行保边滤波的非线性各向异性扩散算法

地震资料的噪音消除是地震资料处理的一个重要步骤,以往的地震资料滤波,都容易伤害地震资料的高频信号.原始地震图像属于线形纹理图像,地震同相轴表现为线条结构.基于偏微分方程的非线性各向异性扩散技术通过结构张量的特征向量控制扩散方向,利用扩散张量的特征值控制特征方向上的扩散量,对原始地震图像进行图像平滑处理,增强具有一致性的同向轴及其所反映的重要地质结构.     

Estimation of uncertainty propagation in flood inundation mapping using a 1‐D hydraulic model

A need for more accurate flood inundation maps has recently arisen because of the increasing frequency and extremity of flood events. The accuracy of flood inundation maps is determined by the uncertainty propagated from all of the variables involved in the overall process of flood inundation modelling. Despite our advanced understanding of flood progression, it is impossible to eliminate the uncertainty because of the constraints involving cost, time, knowledge, and technology. Nevertheless, uncertainty analysis in flood inundation mapping can provide useful information for flood risk management. The twin objectives of this study were firstly to estimate the propagated uncertainty rates of key variables in flood inundation mapping by using the first‐order approximation method and secondly to evaluate the relative sensitivities of the model variables by using the Hornberger–Spear–Young (HSY) method. Monte Carlo simulations using the Hydrologic Engineering Center's River Analysis System and triangle‐based interpolation were performed to investigate the uncertainty arising from discharge, topography, and Manning's n in the East Fork of the White River near Seymour, Indiana, and in Strouds Creek in Orange County, North Carolina. We found that the uncertainty of a single variable is propagated differently to the flood inundation area depending on the effects of other variables in the overall process. The uncertainty was linearly/nonlinearly propagated corresponding to valley shapes of the reaches. In addition, the HSY sensitivity analysis revealed the topography of Seymour reach and the discharge of Strouds Creek to be major contributors to the change of flood inundation area. Copyright © 2014 John Wiley & Sons, Ltd.     

Detection of multiple fracture sets using observations of shear‐wave splitting in microseismic data

The hydrocarbon industry is moving increasingly towards tight sandstone and shale gas resources – reservoirs that require fractures to be produced economically. Therefore, techniques that can identify sets of aligned fractures are becoming more important. Fracture identification is also important in the areas of coal bed methane production, carbon capture and storage (CCS), geothermal energy, nuclear waste storage and mining. In all these settings, stress and pore pressure changes induced by engineering activity can generate or reactivate faults and fractures. P‐ and S‐waves are emitted by such microseismic events, which can be recorded on downhole geophones. The presence of aligned fracture sets generates seismic anisotropy, which can be identified by measuring the splitting of the S‐waves emitted by microseismic events. The raypaths of the S‐waves will have an arbitrary orientation, controlled by the event and geophone locations, meaning that the anisotropy system may only be partly illuminated by the available arrivals. Therefore to reliably interpret such splitting measurements it is necessary to construct models that compare splitting observations with modelled values, allowing the best fitting rock physics parameters to be determined. Commonly, splitting measurements are inverted for one fracture set and rock fabrics with a vertical axis of symmetry. In this paper we address the challenge of identifying multiple aligned fracture sets using splitting measured on microseismic events. We analyse data from the Weyburn CCS‐EOR reservoir, which is known to have multiple fracture sets, and from a hydraulic fracture stimulation, where it is believed that only one set is present. We make splitting measurements on microseismic data recorded on downhole geophone arrays. Our inversion technique successfully discriminates between the single and multiple fracture cases and in all cases accurately identifies the strikes of fracture sets previously imaged using independent methods (borehole image logs, core samples, microseismic event locations). We also generate a synthetic example to highlight the pitfalls that can be encountered if it is assumed that only one fracture set is present when splitting data are interpreted, when in fact more than one fracture set is contributing to the anisotropy.     

常Q衰减介质分数阶波动方程优化有限差分模拟

本文基于Kjartansson常Q模型理论,推导了常Q衰减介质中黏声波和黏弹性波的速度-应力方程,并采用基于二项式窗函数的优化交错网格有限差分方法进行了数值模拟,同时引入不分裂的复频移卷积完全匹配层(CPML)吸收边界条件,以消除边界反射.使用基于自适应时间步长记忆方法的中心差分近似时间分数阶导数,与常用的短时记忆方法相比,提高了波动方程的离散化精度和计算效率.通过对比均匀模型下声波的数值解与解析解,验证了算法的精确性,并进一步分析了不同品质因子下地震波的频散及衰减特征.对BP盐丘模型的数值模拟结果可以较好地反映本文数值方法对复杂介质的适应性及频散压制效果.      

地震数据的反射波动方程最小二乘偏移

基于反射波动方程,本文提出了一种估计地下反射率分布的地震数据最小二乘偏移方法.高频近似下,非齐次的一次反射波动方程的源项是由反射率与入射波场的时间一阶导数相互作用产生的.根据反射波动方程,利用线性最小二乘反演方法由地震反射数据重建出地下产生反射波的反射源,再结合波场正演计算出的地下入射波场,得到地下反射率分布的估计.在地下反射源的线性最小二乘反演重建中,我们采用迭代求解方法,并以地震波的检波器单向地下照明强度作为最小二乘优化问题中Hessian矩阵的近似.