裂缝诱导双相HTI介质模型及其弹性波传播方程

将Biot双相介质理论与Gurevich裂缝各向异性理论相结合,建立了能够同时考虑实际裂缝性储层孔隙性和各向异性的裂缝诱导双相HTI介质模型。从本构方程、动力学方程和动力学达西定律出发,推导出了裂缝诱导双相HTI介质中弹性波传播的一阶速度-应力方程,并针对方程的刚性问题,给出了利用显式二阶时间积分法数值求解该方程时所需要满足的稳定性条件。该方程能够定量地给出双相HTI介质的波场特征与裂缝参数、背景孔隙介质参数之间的关系,描述弹性波在这种介质中的传播机理。     

Time‐scale analysis in unsaturated porous media under external wave loads

Extreme waves caused by tsunamis and storm surges can lead to soil failures in the near‐shore region, which may have severe impact on coastal environments and communities. Multiphase flows in deformable porous media involve several coupled processes and multiple time scales, which are challenging for numerical simulations. The objective of this study is to investigate the roles of the various processes and their interactions in multiphase flows in unsaturated soils under external wave loading, via theoretical time‐scale analysis and numerical simulations. A coupled geomechanics–multiphase flow model based on conservation laws is used. Theoretical analysis based on coupled and decoupled models demonstrates that transient and steady‐state responses are governed by pore pressure diffusion and saturation front propagation, respectively, and that the two processes are essentially decoupled. Numerical simulations suggest that the compressibility of the pore fluids and the deformation of the soil skeleton are important when the transient responses of the media are of concern, while the steady‐state responses are not sensitive to these factors. The responses obtained from the fully coupled numerical simulations are explained by a simplified time‐scale analysis based on coupled and decoupled models. Copyright © 2010 John Wiley & Sons, Ltd.     

原油粘度变化对水驱油开发动态影响的数学模拟方法

为了研究注水开发油田原油粘度升高对开发效果的影响, 通过对实际油藏原油粘度统计, 回归出了原油粘度增长模型.在三维三相黑油渗流模型的基础上, 建立了一个原油粘度随含水和压力变化的油藏渗流数学模型, 并采用有限差分方法建立了相应的数值模型, 采用超松弛法对该模型进行了求解, 用Fortran90语言开发了一个新的数值模拟器.应用该模拟器模拟了不同的原油粘度变化规律对水驱效果的影响, 并与常规模拟器的结果进行了对比.结果表明: 初始水油粘度比为1∶10、含水达到98%时, 粘度增长指数由0增加到0.02, 对应的原油采出程度由44.80%降低到34.29%.目前商业软件中忽略了原油粘度随含水升高而增加的因素, 使得预测的采收率明显偏高.      

Adaptive hierarchical upscaling of flow in heterogeneous reservoirs based on an a posteriori error estimate

This paper treats the upscaling of the absolute permeability in a heterogeneous reservoir. By replacing the fine scale permeability tensor with an upscaled, or effective permeability tensor, a modelling error is introduced. An a posteriori error estimate on this modelling error is formulated and tested. An implementation of the theory, based on domain decomposition coupled with a hierarchical representation of the absolute permeability field, is given. As hierarchical basis functions we have chosen the Haar system, which leads to a wavelet representation of the permeability. The wavelet representation offers a natural upscaling technique which resembles the highcut filters commonly used in signal analysis. This procedure represents an adaptive upscaling method. The numerical results show that this method conserves both the dissipation and the mean velocity in the problem fairly well. The a posteriori error estimate on the modelling error coupled with domain decomposition methods constitutes a powerful modelling tool.     

双复杂介质条件下的反Q滤波偏移延拓算子研究

地下介质的非弹性吸收衰减效应造成地震资料中深层信号频带窄、能量弱、成像困难,因此需要对地层的非弹性吸收衰减进行补偿。在已有的补偿方法基础上,基于单程波方程,采用与品质因子和频率相关的复速度代替介质真实速度,将粘滞传播算子在频率-空间域采用连分式展开,得到相移项和补偿项,针对连分式的误差,在频率-波数域进行补偿校正,并用逐步累加法解决起伏地表问题。该方法提高了成像精度,有效补偿了由于地层非弹性吸收造成的振幅和频率损失。模型试算和实际资料处理表明,该方法具有一定的适应性和实用性。     

A novel approach to leveraging social media for rapid flood mapping: a case study of the 2015 South Carolina floods

Rapid flood mapping is critical for local authorities and emergency responders to identify areas in need of immediate attention. However, traditional data collection practices such as remote sensing and field surveying often fail to offer timely information during or right after a flooding event. Social media such as Twitter have emerged as a new data source for disaster management and flood mapping. Using the 2015 South Carolina floods as the study case, this paper introduces a novel approach to mapping the flood in near real time by leveraging Twitter data in geospatial processes. Specifically, in this study, we first analyzed the spatiotemporal patterns of flood-related tweets using quantitative methods to better understand how Twitter activity is related to flood phenomena. Then, a kernel-based flood mapping model was developed to map the flooding possibility for the study area based on the water height points derived from tweets and stream gauges. The identified patterns of Twitter activity were used to assign the weights of flood model parameters. The feasibility and accuracy of the model was evaluated by comparing the model output with official inundation maps. Results show that the proposed approach could provide a consistent and comparable estimation of the flood situation in near real time, which is essential for improving the situational awareness during a flooding event to support decision-making.     

Stochastic analysis of bounded unsaturated flow in heterogeneous aquifers: Spectral/perturbation approach

This paper describes a stochastic analysis of steady state flow in a bounded, partially saturated heterogeneous porous medium subject to distributed infiltration. The presence of boundary conditions leads to non-uniformity in the mean unsaturated flow, which in turn causes non-stationarity in the statistics of velocity fields. Motivated by this, our aim is to investigate the impact of boundary conditions on the behavior of field-scale unsaturated flow. Within the framework of spectral theory based on Fourier–Stieltjes representations for the perturbed quantities, the general expressions for the pressure head variance, variance of log unsaturated hydraulic conductivity and variance of the specific discharge are presented in the wave number domain. Closed-form expressions are developed for the simplified case of statistical isotropy of the log hydraulic conductivity field with a constant soil pore-size distribution parameter. These expressions allow us to investigate the impact of the boundary conditions, namely the vertical infiltration from the soil surface and a prescribed pressure head at a certain depth below the soil surface. It is found that the boundary conditions are critical in predicting uncertainty in bounded unsaturated flow. Our analytical expression for the pressure head variance in a one-dimensional, heterogeneous flow domain, developed using a nonstationary spectral representation approach [Li S-G, McLaughlin D. A nonstationary spectral method for solving stochastic groundwater problems: unconditional analysis. Water Resour Res 1991;27(7):1589–605; Li S-G, McLaughlin D. Using the nonstationary spectral method to analyze flow through heterogeneous trending media. Water Resour Res 1995; 31(3):541–51], is precisely equivalent to the published result of Lu et al. [Lu Z, Zhang D. Analytical solutions to steady state unsaturated flow in layered, randomly heterogeneous soils via Kirchhoff transformation. Adv Water Resour 2004;27:775–84].     

Surface wave propagation in a transversely isotropic elastic layer overlying a liquid saturated porous solid half-space and lying under the uniform layer of liquid

Dispersion of Rayleigh-type surface wave is studied in a homogeneous transversely isotropic elastic layer overlying a nondissipative liquid-saturated porous solid half-space and lying under a uniform layer of homogeneous liquid. The frequency equation in the form of ninth-order determinant is obtained.Special cases have been deduced by reducing the depth of the layers to zero and by changing the transverse isotropic layer to an isotropic layer. Dispersion curves for the phase velocity have been plotted for a particular model.     

Lattice-Boltzmann simulations of the capillary pressure–saturation–interfacial area relationship for porous media

Hysteresis in the relationship between capillary pressure (P_c), wetting phase saturation (S_w) and nonwetting–wetting interfacial area per volume (a_nw) is investigated using multiphase lattice-Boltzmann simulations of drainage and imbibition in a glass bead porous system. In order to validate the simulations, the P_c–S_w and a_nw–S_w main hysteresis loops were compared to experimental data reported by Culligan et al. [Culligan KA, Wildenschild D, Christensen BS, Gray WG, Rivers ML, Tompson AB. Interfacial area measurements for unsaturated flow through porous media. Water Resour Res 2004;40:W12413]. In general, the comparison shows that the simulations are reliable and capture the important physical processes in the experimental system. P_c–S_w curves, a_nw–S_w curves and phase distributions (within the pores) show good agreement during drainage, but less satisfactory agreement during imbibition. Drainage and imbibition scanning curves were simulated in order to construct P_c–S_w–a_nw surfaces. The root mean squared error (RMSE) and mean absolute error (MAE) between drainage and imbibition surfaces was 0.10 mm⁻¹ and 0.03 mm⁻¹, respectively. This small difference indicates that hysteresis is virtually nonexistent in the P_c–S_w–a_nw relationship for the multiphase system studied here. Additionally, a surface was fit to the main loop (excluding scanning curves) of the drainage and imbibition P_c–S_w–a_nw data and compared to the surface fit to all of the data. The differences between these two surfaces were small (RMSE = 0.05 mm⁻¹ and MAE = 0.01 mm⁻¹) indicating that the P_c–S_w–a_nw surface is adequately represented without the need for the scanning curve data, which greatly reduces the amount of data required to construct the non-hysteretic P_c–S_w–a_nw surface for this data.