首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
This paper part one is set out to lay primary observations of experimental compaction measurements to form the basis for rock physics modelling in paper part two. P- and S-wave velocities and corresponding petrophysical (porosity and density) properties of seven unconsolidated natural sands with different mineralogical compositions and textures are reported. The samples were compacted in a uniaxial strain configuration from 0.5 up to 30 MPa effective stresses. Each sand sample was subjected to three loading cycles to study the influence of stress reduction on acoustic velocities and rock physical properties with the key focus on simulating a complex burial history with periods of uplift. Results show significant differences in rock physical properties between normal compaction and overconsolidation (unloaded and reloaded). The differences observed for total porosity, density, and P- and S-wave velocities are attributed to irrecoverable permanent deformation. Microtextural differences affect petrophysical, acoustic, elastic and mechanical properties, mostly during normal consolidation but are less significant during unloading and reloading. Different pre-consolidation stress magnitudes, stress conditions (isotropic or uniaxial) and mineral compositions do not significantly affect the change in porosity and velocities during unloading as a similar steep velocity–porosity gradient is observed. The magnitude of change in the total porosity is low compared to the associated change in P- and S-wave velocities during stress release. This can be explained by the different sensitivity of the porosity and acoustic properties (velocities) to the change in stress. Stress reduction during unloading yields maximum changes in the total porosity, P- and S-wave velocities of 5%, 25%, and 50%, respectively. These proportions constitute the basis for the following empirical (approximation) correlations: Δϕ ∼ ±5 ΔVP and ΔVP ∼ ±2ΔVS. The patterns observed in the experiments are similar to well log data from the Barents Sea. Applications to rock physics modelling and reservoir monitoring are reported in a companion paper.  相似文献   

2.
Seismic monitoring of reservoir and overburden performance during subsurface CO2 storage plays a key role in ensuring efficiency and safety. Proper interpretation of monitoring data requires knowledge about the rock physical phenomena occurring in the subsurface formations. This work focuses on rock stiffness and elastic velocity changes of a shale overburden formation caused by both reservoir inflation induced stress changes and leakage of CO2 into the overburden. In laboratory experiments, Pierre shale I core plugs were loaded along the stress path representative for the in situ stress changes experienced by caprock during reservoir inflation. Tests were carried out in a triaxial compaction cell combining three measurement techniques and permitting for determination of (i) ultrasonic velocities, (ii) quasistatic rock deformations, and (iii) dynamic elastic stiffness at seismic frequencies within a single test, which allowed to quantify effects of seismic dispersion. In addition, fluid substitution effects connected with possible CO2 leakage into the caprock formation were modelled by the modified anisotropic Gassmann model. Results of this work indicate that (i) stress sensitivity of Pierre shale I is frequency dependent; (ii) reservoir inflation leads to the increase of the overburden Young's modulus and Poisson's ratio; (iii) in situ stress changes mostly affect the P‐wave velocities; (iv) small leakage of the CO2 into the overburden may lead to the velocity changes, which are comparable with one associated with geomechanical influence; (v) non‐elastic effects increase stress sensitivity of an acoustic waves; (iv) and both geomechanical and fluid substitution effects would create significant time shifts, which should be detectable by time‐lapse seismic.  相似文献   

3.
A 3D relocation technique permits precise locations of induced earthquakes. Geostatistical processing using the data of 87 boreholes provides the basis of a precise 3D structure, with a dome geometry. Conventional laboratory mechanical tests performed on deep rock samples (1000 m to 5000 m) define the rock properties at depths similar to those of the seismic events (1<M L<4.2) that range from 1 to 7 km.In the studied period, most (85%), of the events were located above the gas reservoir, with very few located in the reservoir itself. Because the production parameters (50 MPa depletion of the gas pressure reservoir) are homogeneous throughout the gas field, the lateral inhomogeneity of the seismic rupture locations are a consequence of variations in the rheological response of the dome to the deformation induced by gas production.Here a ratio of two is found between the elastic modulus of the seismic rock matrix and the elastic modulus of the aseismic rock matrix. The contrast in strength is at least as great, if not greater. Repeated measured surface deformations involve the whole structure. Spatial and temporal deformations indicate that aseismic deformation is quantitatively the main process of this structural deformation. The heterogeneous stress pattern inferred fromP-axes of induced earthquakes disagrees with the tectonic regional stress field. The radial distribution ofP-axes towards the gas reservoir probably reflects the production induced deformation. The inferred deformation of the dome occurs in response to weak induced stresses.  相似文献   

4.
使用基于机器学习构建的贵州西部多流域水库区地震目录和震相报告,采用波速比一致性约束的双差层析成像方法联合反演得到了该地区的地震位置和三维vp、vs、vp/vs结构。结果显示,研究区内的地壳速度结构具有明显的不均匀性,不同大地构造变形分区展现出不同的速度结构特征。在10km以上深度,横穿黔西中部NW向的威宁构造变形区显示出显著的低速异常条带,揭示了威宁—水城深大断裂带的影响深度和范围。受库区岩性和流体渗透的影响,0km深度的速度结构显示普遍的低波速和高波速比特征,包括夹岩、平寨、光照和马马崖水库区域。重定位后,地震的空间分布勾勒出大量隐伏断层的几何展布特征,结合三维速度结构,推测该区域的地震活动与水库周围的断层活化有关。  相似文献   

5.
Concepts from the Theory of Interacting Continua are employed to develop constitutive relations for thermoelastic fluid-saturated porous media; these constitutive relations are shown to be equivalent to those of Biot in isothermal limit. For uniaxial compaction, the constitutive theory is simplified to give formation compaction as a function of changes in fluid pressure and temperature. The formation compaction due to thermal cooling is equal to 3hη[(1 + v)/3(1 − v)]ΔT.  相似文献   

6.
We review earthquake distributions associated with hydrocarbon fields in the context of pore pressure diffusion models, poroelastic stress transfer and isostasy theory. These three mechanisms trigger or induce seismic instabilities at both local scale (D5 km) and at regional scale (D20 km). The modeled changes in stress are small (1 MPa), whatever the tectonic setting. Each mechanism corresponds to different production processes. (1) Local hydraulic fracturing due to fluid injection induces seismic-slip on cracks (M L3) within the injected reservoir through decreasing the effective stress. (2) Pure fluid withdrawal causes pore pressure to decrease within the reservoir. It triggers adjustments of the geological structure to perturbations related to the reservoir response to depletion. Poroelastic mechanisms transfer this stress change from the reservoir to the surrounding levels whereM L5 seismic instabilities occur either above or below the reservoir. (3) Massive hydrocarbon recovery induces crustal readjustments due to the removal of load from the upper crust. It can induce larger earthquakes (M L6) at greater distance from the hydrocarbon fields than the two other mechanisms.Due to the mechanical properties of the shallow rock matrices involved, seismic slip triggered either by mechanism (1) or (2), is a second-order process of the main elastoplastic deformation. for a minimum of 80% of commercially productive basins, most of the local deformation is reported as aseismic, i.e., there is no evidence forM L3 earthquakes. Nevertheless, the induced stresses vary as a function of time in a manner that depends on the hydraulic diffusivity (i.e., permeability) of the reservoir and surrounding rocks. Because small earthquakes (M L3) indicate changes in stress and pore pressure, monitoring of seismicity is a means of assessingin situ reservoir behavior.The less constrained seismic response to hydrocarbon recovery is the possible connection between local fluid manipulations, triggered earthquakes and major regional earthquakes. Positive feedback mechanisms suggest that the region of seismic hazard changes is much larger than the area where hydrocarbons are extracted. These observations and models testify that fluid movement and pore pressure changes (increase or decrease) play important roles in the mechanics of earthquakes and in the triggering of natural earthquakes.  相似文献   

7.
The Tian Shan is a vast range that spans several countries in Asia. Understanding its evolutionary history may provide valuable insights into intracontinental orogenic dynamics. In this study, we explored the crustal characteristics of the Tian Shan and their relationships to the tectonic evolution of the region. A new H-stacking method that combines the P receiver function and gravity anomalies was used to estimate the thickness and ratio of P- to S-wave velocities (vP/vS) for 91 broadband seismic stations in the central and western Tian Shan. Our results revealed significant lateral variations in crustal thickness and vP/vS. A ~45-km-thick crust and an intermediate-high vP/vS (~1.74–1.84) were found in the Kazakh Shield and Tarim Basin, which we interpreted to indicate a mafic crystalline basement and lower crust. The central Tian Shan varied greatly in crustal thickness (40–64 km) and vP/vS ratio (1.65–2.00), which may be due to crustal shortening, mafic underplating, and crustal melting. In contrast, we observed a relatively thin crust (42–50 km) with an intermediate vP/vS ratio (~1.78) in the western Tian Shan. The differences in the crustal structures between the western and central Tian Shan imply that the Talas-Fergana Fault may be trans-lithospheric.  相似文献   

8.
Velocities of compressional and shear waves in limestones   总被引：2，自引：1，他引：2  
Carbonate rocks are important hydrocarbon reservoir rocks with complex textures and petrophysical properties (porosity and permeability) mainly resulting from various diagenetic processes (compaction, dissolution, precipitation, cementation, etc.). These complexities make prediction of reservoir characteristics (e.g. porosity and permeability) from their seismic properties very difficult. To explore the relationship between the seismic, petrophysical and geological properties, ultrasonic compressional‐ and shear‐wave velocity measurements were made under a simulated in situ condition of pressure (50 MPa hydrostatic effective pressure) at frequencies of approximately 0.85 MHz and 0.7 MHz, respectively, using a pulse‐echo method. The measurements were made both in vacuum‐dry and fully saturated conditions in oolitic limestones of the Great Oolite Formation of southern England. Some of the rocks were fully saturated with oil. The acoustic measurements were supplemented by porosity and permeability measurements, petrological and pore geometry studies of resin‐impregnated polished thin sections, X‐ray diffraction analyses and scanning electron microscope studies to investigate submicroscopic textures and micropores. It is shown that the compressional‐ and shear‐wave velocities (Vp and Vs, respectively) decrease with increasing porosity and that Vp decreases approximately twice as fast as Vs. The systematic differences in pore structures (e.g. the aspect ratio) of the limestones produce large residuals in the velocity versus porosity relationship. It is demonstrated that the velocity versus porosity relationship can be improved by removing the pore‐structure‐dependent variations from the residuals. The introduction of water into the pore space decreases the shear moduli of the rocks by about 2 GPa, suggesting that there exists a fluid/matrix interaction at grain contacts, which reduces the rigidity. The predicted Biot–Gassmann velocity values are greater than the measured velocity values due to the rock–fluid interaction. This is not accounted for in the Biot–Gassmann velocity models and velocity dispersion due to a local flow mechanism. The velocities predicted by the Raymer and time‐average relationships overestimated the measured velocities even more than the Biot model.  相似文献   

9.
One of the fundamental questions in modern volcanology is the manner in which a volcanic eruption is triggered; the intrusion of fresh magma into a reservoir is thought to be a key component. The amount by which previously ponded reservoir magma interacts with a newly intruded magma will determine the nature and rate of eruption as well as the chemistry of erupted lavas and shallow dykes. The physics of this interaction can be investigated through a conventional monitoring procedure that incorporates the simple and much used Mogi model relating ground deformation (most simply represented by Δh) to changes in volume of a magma reservoir. Gravity changes (Δg) combined with ground deformation provide information on magma reservoir mass changes. Our models predict how, during inflation, the observed Δg/Δh gradient will evolve as a volcano develops from a state of dormancy through unrest into a state of explosive activity. Calderas in a state of unrest and large composite volcanoes are the targets for the methods proposed here and are exemplified by Campi Flegrei, Rabaul, Krafla, and Long Valley. We show here how the simultaneous measurement of deformation and gravity at only a few key stations can identify important precursory processes within a magma reservoir prior to the onset of more conventional eruption precursors.  相似文献   

10.
A sequence of moderate shallow earthquakes (3.5M L5.3) was located within the Vercors massif (France) in the period 1961–1984. This subalpine massif has been a low seismic area for at least 5 centuries. During the period 1962–1963, 12 shallow earthquakes occurred in the neighborhood (10 km) of the Monteynard reservoir, 30 km south of the city of Grenoble. The latest fourM L4.0 earthquakes occurred in 1979–1984 either at larger distance (35 km) or greater depth (10 km) from the reservoir. Two triggering mechanisms are suggested for this sequence: (i) the direct effect of elastic loading through either increased shear stress or strength reducing by increased pore pressure at depth; (ii) the pore pressure diffusion induced by poroelastic stress change due to the reservoir filling.The weekly water levels, local balanced geological cross sections, and focal mechanisms argue for two types of mechanical connection between the earthquake sequence and the filling cycles of the Monteynard reservoir. The seismic sequence started with the 1962–1963 shallow earthquakes that occurred during the first filling of the reservoir and are typical of the direct effect of elastic loading. The 1979 deeper earthquake is located at a 10 km depth below the reservoir. This event occurred 16 years after the initial reservoir impoundment, but one month after the previous 1963 maximum water level was exceeded. Moreover the yearly reservoir level increased gradually in the period 1962–1979 and has decreased since 1980. Accordingly we suggest that the gradual diffusion of water from reservoir to hypocentral depths decreases the strength of the rock matrices through increased pore pressure. The transition between the two types of seismic response is supported by the analysis ofM L3.5 earthquakes which all occurred in the period 1964–1971, ranging between 10 and 30 km distance from the reservoir. The three other delayed earthquakes of the 1961–1984 seismic sequence (M L4 during the 1979–1984 period) are all located 35 km away from the reservoir. Based on the seismic activity, the estimates for the hydraulic diffusivities range between 0.2–10 m2/s, except for the first event that occurred 30 km north of the reservoir, the filling just started. The lack ofin situ measurements of crustal hydrological properties in the area, shared by most of the Reservoir-Induced-Seismicity cases, prevents us from obtaining absolute evidence for the triggering processes. These observations and conceptual models attest that previous recurrence times for moderate natural shocks (4.5M L5.5) estimated within this area using historical data, could be modified by 0.1–1 MPa stress changes. These small changes in deviatoric stress suggest that the upper crust is in this area nearly everywhere at a state of stress near failure. Although the paucity of both number and size of earthquakes in the French subalpine massif shows that aseismic displacements prevail, our study demonstrates that triggered earthquakes are important tools for assessing local seismic risk through mapping fault zones and identifying their possible seismic behavior.  相似文献   

11.
An approximation is developed that allows mapped 4D seismic amplitudes and time‐shifts to be related directly to the weighted linear sum of pore pressure and saturation changes. The weights in this relation are identified as key groups of parameters from a petroelastic model and include the reservoir porosity. This dependence on groups of parameters explains the inherent non‐uniqueness of this problem experienced by previous researchers. The proposed relation is of use in 4D seismic data feasibility studies and inversion and interpretation of the 4D seismic response in terms of pore pressure and water saturation changes. A further result is drawn from analysis of data from the North Sea and West Africa, which reveals that the relative interplay between the effects of pore pressure and saturation changes on the seismic data can be simplified to the control of a single, spatially variant parameter CS/CP. Combining these results with those from published literature, we find that CS/CP = 8 appears to be a generality across a range of clastic reservoirs with a similar mean porosity. Using this CS/CP value, an in situ seismic‐scale constraint for the rock stress sensitivity component of the petroelastic model is constructed considering this component carries the largest uncertainty.  相似文献   

12.
Pore-pressure depletion causes changes in the triaxial stress state. Pore-pressure depletion in a flat reservoir, for example, can be reasonably approximated as uniaxial compaction, in which the horizontal effective stress change is smaller than the vertical effective stress. Furthermore, the stress sensitivity of velocities can be angle-dependent. Therefore, time-lapse changes in reservoir elastic anisotropy are expected as a consequence of production, which can complicate the interpretation of the 4D seismic response. The anisotropic 4D seismic response caused by pore-pressure depletion was investigated using existing core velocity measurements. To make a direct comparison between the anisotropic 4D seismic response and the isotropic response based only on vertical velocities, pseudoisotropic elastic properties were utilized, and the two responses were compared in terms of a dynamic rock physics template. A comparison of the dynamic rock physics templates indicates that time-lapse changes in reservoir elastic anisotropy have a noticeable impact on the interpretation of 4D seismic data. Changes in anisotropy as a result of pore-pressure depletion cause a time-lapse amplitude variation with offset response as if there is a reduction in VP/VS (i.e., pseudoisotropic VP/VS decreases), although the vertical VP/VS increases. The impact of time-lapse changes in anisotropy on the amplitude variation with offset gradient was also investigated, and the time-lapse anisotropy was found to enhance changes in the amplitude variation with offset gradient for a given case.  相似文献   

13.
The theme of the 2003 EAGE/SEG imaging workshop concerned the contrast between different philosophies of ‘model building’: whether an explicit, user‐determined model should be imposed throughout the processing, with user updates at each step; or alternatively, whether user intervention should be kept to a minimum so as to avoid preconceived bias, and instead to allow the data itself to guide some heuristic process to converge to an optimal solution. Here we consider a North Sea study where our initial approach was to build the subsurface model using interpreted horizons as a guide to the velocity update. This is common practice in the North Sea, where the geology ‘lends itself’ to a layer‐based model representation. In other words, we encourage preconceived bias, as we consider it to be a meaningful geological constraint on the solution. However, in this instance we had a thick chalk sequence, wherein the vertical compaction gradient changed subtly, in a way not readily discernible from the seismic reflection data. As a consequence, imposing the explicit top and bottom chalk horizons, with an intervening vertical compaction gradient (of the form v(x, y, z) =v0(x, y) +k(x, y).z), led to a misrepresentation of the subsurface. To address this issue, a gridded model building approach was also tried. This relied on dense continuous automatic picking of residual moveout in common‐reflection point gathers at each iteration of the model update, followed by gridded tomography, resulting in a smoothly varying velocity field which was able to reveal the underlying local changes within the chalk.  相似文献   

14.
An elastic point source model proposed by Mogi for magma chamber inflation and deflation has been applied to geodetic data collected at many volcanoes. The volume of ground surface uplift or subsidence estimated from this model is closely related to the volume of magma injection into or withdrawal from the reservoir below. The analytical expressions for these volumes are reviewed for a spherical chamber and it is shown that they differ by the factor 2(1-v), where v is Poisson's ratio of the host rock. For the common estimate v=0.25, as used by Mogi and subsequent workers, the uplift volume is 3/2 the injection volume. For highly fractured rocks, v can be even less and the uplift volume can approach twice the injection volume. Unfortunately, there is no single relation between the inflation of magma reservoirs and the dilation or contraction of host rocks. The inflation of sill-like bodies, for instance, generates no overall change in host rock volume. Inflation of dike-like bodies generates contraction such that, in contrast with Mogi's result, the uplift volume is generally less than the injection volume; for v=0.25, the former is only 3/4 of the latter. Estimates of volumes of magma injection or withdrawal are there-fore greatly dependent on the magma reservoir configuration. Ground surface tilt data collected during the 1960 collapse of Kilauea crater, one of the first events interpreted with Mogi's model and one of the largest collapses measured at Kilauea, is not favored by any one of a variety of deformation models. These models, however, predict substantially different volumes of both magma withdrawal and ground surface subsidence.  相似文献   

15.
Theoretical constraints on the stress-dilation relation for a deforming Coulomb material requirev ifC=0 andv sin-1( m / m ) always, wherev is the dilation angle, is the friction angle,C is cohesion, m is the maximum shear stress, and m is the mean effective stress. Recent laboratory measurements of friction and dilatancy of simulated fault gouge show that small amplitude shear-load cycling causes compaction and consolidation. Comparison of the data with theory indicates that such load cycling produces: (1) increased coefficient of friction (or friction angle), (2) increased cohesion, and (3) increased dilatancy rate (or dilation angle). Under certain conditions of load cycling without significant plastic shear strain accumulation ( p <0.005) we find thatv exceeds both and, in contrast to theory, sin-1( m / m ). This result is interpreted in terms of enhanced cohesion and overconsolidation, which lead to residual stresses within the gouge. An analogy is drawn between these special loading conditions and those extant on natural faults. In particular, our results imply that jostling and minor stress variations associated with microearthquakes may produce strengthening of fault gouge and changes in the fault zone's stress-dilatancy relation. Hence, compaction associated with microseismicity may lead to subsequent dilation of fault gouge, even for faults with large displacement rates and large net offsets (e.g., San Andreas). In regions where such dilation persists over sufficient displacements (on the order of the critical slip distance for seismic faulting) it may tend to inhibit unstable slip.  相似文献   

16.
For converted waves stacking requires a true common reflection point gather which, in this case, is also a common conversion point (CCP) gather. We consider converted waves of the PS- and SP-type in a stack of horizontal layers. The coordinates of the conversion points for waves of PS- or SP-type, respectively, in a single homogeneous layer are calculated as a function of the offset, the reflector depth and the velocity ratio vp/vs. Knowledge of the conversion points enables us to gather the seismic traces in a common conversion point (CCP) record. Numerical tests show that the CCP coordinates in a multilayered medium can be approximated by the equations given for a single layer. In practical applications, an a priori estimate of vp/vs is required to obtain the CCP for a given reflector depth. A series expansion for the traveltime of converted waves as a function of the offset is presented. Numerical examples have been calculated for several truncations. For small offsets, a hyperbolic approximation can be used. For this, the rms velocity of converted waves is defined. A Dix-type formula, relating the product of the interval velocities of compressional and shear waves to the rms velocity of the converted waves, is presented.  相似文献   

17.
南堡凹陷C油田沙河街组发育一套披覆在潜山之上的生屑云岩储层,为该油田的主力含油层系,生屑云岩分布不稳定,横向变化较大,先期钻探的5口探井有2口落空。基于周边油田分析,古地貌是生屑云岩储层分布的主要因素。受上覆火成岩的影响,该油藏速度横向变化剧烈,古地貌恢复困难,使预测储层分布范围困难较大。考虑到速度横向变化的影响,利用传统的时深转换方法难以准确地计算构造深度进而恢复古地貌。本文从多口已钻井的时深关系精细对比分析出发,明确了该区的速度影响因素主要有压实效应和火成岩厚度,同时基于对传统时深转换方法的适用性分析,提出了一种基于二元回归的时深转换新方法。实际资料应用表明,利用该方法计算的构造深度与实钻结果吻合较好,能够有效恢复古地貌和预测生屑云岩的分布,该时深转换方法对于速度横向变化剧烈的地区有较好的推广应用价值。   相似文献   

18.
Relative seismic velocity change (dv/v) is important for monitoring changes in subsurface material properties and evaluating earthquake-induced rock slope damage in a geological disaster-prone region. In this paper, we present a rapid damage assessment on three slow-moving rock slopes by measuring dv/v decrease caused by the 2022 ?MS 6.8 Luding earthquake in Southwest China. By applying the stretching method to the cross-correlated seismic wavefields between sensors installed on each slope, we obtain earthquake-induced dv/v decreases of ～2.1%, ～0.5%, and ～0.2% on three slopes at distances ranging from ～86 to ～370 ?km to the epicenter, respectively. Moreover, based on seismic data recorded by 16 sensors deployed on the rock slope at a distance of ～370 ?km away from the epicenter, a localized dv/v decease region was observed at the crest of the slope by calculating the spatial dv/v images before and after the earthquake. We also derive an empirical in situ stress sensitivity of ?7.29?10?8/Pa by relating the dv/v change to the measured peak dynamic stresses. Our results indicate that a rapid dv/v assessment not only can help facilitate on-site emergency response to earthquake-induced secondary geological disasters but also can provide a better understanding of the subsurface geological risks under diverse seismic loadings.  相似文献   

19.
Differential compaction has long been used by seismic interpreters to infer subsurface geology using knowledge of the relative compaction of different types of sediments. We outline a method to infer the gross fraction of shale in an interval between two seismic horizons using sandstone and shale compaction laws. A key component of the method involves reconstruction of a smooth depositional horizon by interpolating decompacted thicknesses from well control. We derive analytic formulae for decompaction calculations using known porosity–stress relations and do not employ discrete layer iterative methods; these formulae were found to depend not only upon the gross fraction of shale but also on the clay content of the shales and the thickness of the interval. The relative merits of several interpolation options were explored, and found to depend upon the structural setting. The method was successfully applied to an oil sands project in Alberta, Canada.  相似文献   

20.
Legacy streamer data and newer 3D ocean‐bottom‐cable data are cross‐matched and analysed for time‐lapse analysis of geomechanical changes due to production in the Valhall Field. The issues relating to time‐lapse analysis using two such distinctly different data sets are addressed to provide an optimal cross‐matching workflow that includes 3D warping. Additionally an assessment of the differences between the imaging using single‐azimuth streamer and multi‐azimuth ocean‐bottom‐cable data is provided. The 3D warping utilized in the cross‐matching procedure is sensitive to acquisition and processing differences but is also found to provide valuable insight into the geometrical changes that occur in the subsurface due to production. As such, this work also provides a demonstration of the use of high‐resolution 3D interpreted warping to resolve the 3D heterogeneity of the compaction and subsidence. This is an important tool for Valhall, and possibly other fields, where compaction and subsidence (and monitoring thereof) are key factors in the reservoir management since the predominant observed production‐induced changes are compaction of the soft, high‐porosity chalk reservoir, due to pore‐pressure reduction, and the resultant overburden subsidence. Such reservoir compaction could have significant implications for production by changing permeabilities and production rates. Furthermore the subsidence effects could impact upon subsea installations and well‐bore stability. Geomechanical studies that have previously been used to model such subsidence and compaction are only constrained by observed surface displacements and measured reservoir pressure changes, with the geological overburden being largely neglected. The approaches suggested herein provide the potential for monitoring and assessment in three dimensions, including the probable heterogeneity and shearing, that is needed for full understanding of reservoir compaction and the resultant effects on the overburden to, for example, mitigate well‐bore failures.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号