地震尺度下碳酸盐岩储层的岩石物理建模方法(英文)

碳酸盐岩油藏的强非均质性以及孔隙结构的复杂性,使得作为连接油藏参数与地震参数重要桥梁的岩石物理模型,以及作为油藏预测和定量表征最有效工具的流体替换成为岩石物理建模的难点与重点。在碳酸盐岩储层复杂孔隙结构与地震尺度下碳酸盐岩储层非均质性分析基础上,研究采用岩石网格化方法,将地震尺度下非均质碳酸盐岩储层岩石划分为具有独立岩石参数的均质岩石子体,根据岩石孔隙成因与结构特征采用不同岩石物理模型分步计算岩石子块干岩石弹性模量,并根据不同孔隙连通性进行流体替换,计算饱和不同流体岩石弹性模量。基于计算的岩石子块弹性模量,采用Hashin-Shtrikman-Walpole弹性边界计算理论方法实现地震尺度下碳酸盐岩储层弹性参数计算。通过对含有不同类型孔隙组合碳酸盐岩储层模型的弹性模量进行计算与分析,明确不同孔隙对岩石弹性参数的影响特征,模拟分析结果与实际资料认识一致。     

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second-order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite-element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.     

三维非均匀地质模型中的逐段迭代射线追踪

地震射线追踪是地震定位、层析成像、偏移等领域的重要正演环节.随着这些领域研究的深入,针对传统的网格结构和层状结构在描述复杂地质模型遇到的很大困难,我们采用大小不等、形状各异的地质块组成的集合体来描述三维复杂地质模型,并用三角形面片来描述地质块之间的物性间断面,理论上可以描述任意复杂的地质模型.为适应任意非均匀速度分布的地质模型,基于费马原理,本文发展了与之相适应的逐段迭代射线追踪方法.该方法属于弯曲法范畴,对路径点采用一阶显式增量修正,相对于传统的迭代法,高效省时.数值试验表明,联合逐段迭代法和伪弯曲法的射线追踪扰动修正方案在三维复杂非均匀块状模型中有适用性和高效性.     

Predictive modeling of the permafrost thermal regime in Russian railroad subgrade support systems

The goal of a predictive thermotechnical calculation is to model the behavior of the top permafrost boundary under current operational conditions as well as increasing average annual air temperatures that results in degradation of the permafrost layer. Numerical modeling was used to assess the efficient application of construction measures to create sustainable operation of the railroad. The numerical modeling was carried out in the programming complex FEM-models developed by geotechnical engineers of St. Petersburg, Russia under Prof. V. M. Ulitsky’s guidance. The Termoground Program as a part of the FEM-models enables the research of freezing, heaving and thawing in different design solutions. Research was carried out in space resolution for a year cycle. The performed model has shown that the designing measures accepted for permafrost protection from retreat in the subgrade support were generally effective.     

Integrated assessment for establishing an oil environmental vulnerability map: Case study for the Santos Basin region,Brazil

The growth of maritime transport and oil exploitation activities may increase the risk of oil spills. Thus, plans and actions to prevent or mitigate impacts are needed to minimize the effects caused by oil. However, tools used worldwide to support contingency plans have not been integrated, thus leading to failure in establishing priority areas. This investigation aimed to develop indices of environmental vulnerability to oil (IEVO), by combining information about environmental sensibility to oil and results of numerical modeling of spilled oil. To achieve that, a case study concerning to oil spills scenarios in a subtropical coastal area was designed, and IEVOs were calculated and presented in maps, in order to make the information about the areas’ vulnerability more easily visualized. For summer, the extension of coastline potentially affected by oil was approximately 150 km, and most of the coastline presented medium to high vulnerability. For winter, 230 km coastline would be affected, from which 75% were classified as medium to high vulnerability. Thus, IEVO maps allowed a rapid and clearer interpretation of the vulnerability of the mapped region, facilitating the planning process and the actions in response to an oil spill.     

A numerical study of the behavior of fully grouted rockbolts under dynamic loading

The widespread preference for rockbolts in providing adequate rock stability in underground mines and man-made cavities make it necessary to obtain a better understanding of the response of these support systems. It is widely accepted that the fully grouted rockbolts provide better roof stability in areas of very poor roof condition which may be caused by high ground stress conditions. There is little information about the in situ behavior of these systems especially under dynamic loading. Hence, to study the behavior of fully grouted rockbolts under dynamic loading, a numerical modeling study was conducted using the FLAC^3D code. In this study the behavior of three types of fully grouted rockbolts were compared with each other including rockbolts with and without head plate and a yielding type one. The results of analyses indicated that under dynamic loading the fully grouted rockbolts without the head plate are incapable of controlling the rock mass movement. Although fully grouted rockbolts with head plate damp a considerable amount of the dynamic energy through friction as these bolts slide within the grout, but the elongation of the rockbolt is not possible due to the rapid breakage of the rockbolt–grout interface. Yielding rockbolts are the best choice for the absorption of the dynamic stress wave and controlling of the rock mass movement. The obtained results show that the optimal design of yielding rockbolt should be in such a way that after dynamic loading, anchoring part of bolt has a limited movement to prevent stress concentration in the shaft of rockbolt and its breakage.     

AVO analysis with partial stacking to detect gas anomalies in the GÜEPAJÉ-3D project

We applied analysis AVO with partial stacking to an onshore gas field in order to detect gas at the top of the Ciénaga de Oro formation (FCO), by using seismic gathers of the Güepajé-3D project. The technique was supported by petrophysical analysis of the Ayombe-1 well.The sensitivity of seismic response to changes in the saturating fluid was evaluated by fluid substitution technique in an interval of the Ayombe-1 well. As a result, a class I AVO anomaly at the top of the FCO was observed, even though the gas and water responses were similar. After applying AVO with partial stack technique at the top of FCO, AVO class I anomalies were observed in the Ayombe-1 and Güepajé-1 gas wells and none in the dry Güepajé-3 well. In spite of being reported with gas, the results in the Güepajé-2 well indicated any to an incipient class I anomaly related to a weak seismic response associated to gas.When the technique was applied to the seismic volume of the Güepajé-3D Project, a map of the direct gas indicator was obtained. The map shows a high correlation with the structural surface at top of FCO. Lithic factors such as compaction, cementation and carbonate overlaying the FCO and variable thickness of partial and full saturated layers, may cause anomalies in zones unlikely to contain gas.     

The impact of subsurface conceptualization on land energy fluxes

There is a significant body of work demonstrating the importance of hydrologic control on land energy feedbacks. Yet, quantitative data on aquifer conductivity can be difficult to assemble. Furthermore, how subsurface uncertainty propagates into land-surface processes is not well understood. This study analyzes the impact of aquifer characterization on land energy fluxes, using a coupled hydrology–land-surface model. Four gridded subsurface conductivity fields are developed for the Upper Klamath basin using two data sources and different levels of imposed heterogeneity. Each model is forced with the same transient, observed meteorology for 3 years prior to the final year presented here. Results are analyzed to quantify the impact of subsurface heterogeneity on groundwater surface water interactions and spatial patterns in hydrologic variables. Analysis shows that heterogeneity does not fundamentally alter the connection between groundwater and land surface processes. However, differences between scenarios impact the extent and location of the critical zone.     

Predicting Diel,Diurnal and Nocturnal Dynamics of Dissolved Oxygen and Chlorophyll‐a Using Regression Models and Neural Networks

Human‐induced and natural interruptions with continuous streams of observational data necessitate the development of gap‐filling and prediction strategies towards better understanding, monitoring and management of aquatic systems. This study quantified the efficacy of multiple non‐linear regression (MNLR) versus artificial neural network (ANN) models as well as the temporal partitioning of diurnal versus nocturnal data for the predictions of chlorophyll‐a (chl‐a) and dissolved oxygen (DO) dynamics. The temporal partitioning increased the predictive performances of the best MNLR models of diurnal DO by 45% and nocturnal DO by 4%, relative to the best diel MNLR model of diel DO ($r_{{\rm adj}}^{2} = 68.8\%$ ). The ANN‐based predictions had a higher predictive power than the MNLR‐based predictions for both chl‐a and DO except for diurnal DO dynamics. The best ANNs based on independent validations were multilayer perceptron (MLP) for diel chl‐a, generalized feedforward (GFF) for diurnal and nocturnal chl‐a, MLP for diel DO, GFF for diurnal DO, and MLP for nocturnal DO.