首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 23 毫秒
1.
Burial stress on a sediment or sedimentary rock is relevant for predicting compaction or failure caused by changes in, e.g., pore pressure in the subsurface. For this purpose, the stress is conventionally expressed in terms of its effect: “the effective stress” defined as the consequent elastic strain multiplied by the rock frame modulus. We cannot measure the strain directly in the subsurface, but from the data on bulk density and P‐wave velocity, we can estimate the rock frame modulus and Biot's coefficient and then calculate the “effective vertical stress” as the total vertical stress minus the product of pore pressure and Biot's coefficient. We can now calculate the elastic strain by dividing “effective stress” with the rock frame modulus. By this procedure, the degree of elastic deformation at a given time and depth can be directly expressed. This facilitates the discussion of the deformation mechanisms. The principle is illustrated by comparing carbonate sediments and sedimentary rocks from the North Sea Basin and three oceanic settings: a relatively shallow water setting dominated by coarse carbonate packstones and grainstones and two deep water settings dominated by fine‐grained carbonate mudstones and wackestones.  相似文献   

2.
Seismic attenuation mechanisms receive increasing attention for the characterization of fractured formations because of their inherent sensitivity to the hydraulic and elastic properties of the probed media. Attenuation has been successfully inferred from seismic data in the past, but linking these estimates to intrinsic rock physical properties remains challenging. A reason for these difficulties in fluid-saturated fractured porous media is that several mechanisms can cause attenuation and may interfere with each other. These mechanisms notably comprise pressure diffusion phenomena and dynamic effects, such as scattering, as well as Biot's so-called intrinsic attenuation mechanism. Understanding the interplay between these mechanisms is therefore an essential step for estimating fracture properties from seismic measurements. In order to do this, we perform a comparative study involving wave propagation modelling in a transmission set-up based on Biot's low-frequency dynamic equations and numerical upscaling based on Biot's consolidation equations. The former captures all aforementioned attenuation mechanisms and their interference, whereas the latter only accounts for pressure diffusion phenomena. A comparison of the results from both methods therefore allows to distinguish between dynamic and pressure diffusion phenomena and to shed light on their interference. To this end, we consider a range of canonical models with randomly distributed vertical and/or horizontal fractures. We observe that scattering attenuation strongly interferes with pressure diffusion phenomena, since the latter affect the elastic contrasts between fractures and their embedding background. Our results also demonstrate that it is essential to account for amplitude reductions due to transmission losses to allow for an adequate estimation of the intrinsic attenuation of fractured media. The effects of Biot's intrinsic mechanism are rather small for the models considered in this study.  相似文献   

3.
Shales comprise more than 60% of sedimentary rocks and form natural seals above hydrocarbon reservoirs. Their sealing capacity is also used for storage of nuclear wastes. The world's most important conventional oil and gas reservoirs have their corresponding source rocks in shale. Furthermore, shale oil and shale gas are the most rapidly expanding trends in unconventional oil and gas. Shales are notorious for their strong elastic anisotropy, i.e., so‐called vertical transverse isotropy. This vertical transverse isotropy, characterised by a vertical axis of invariance, is of practical importance as it is required for correct surface seismic data interpretation, seismic to well tie, and amplitude versus offset analysis. A rather classical paradigm makes a clear link between compaction in shales and the alignment of the clay platelets (main constituent of shales). This would imply increasing anisotropy strength with increasing compaction. Our main purpose is to check this prediction on two large databases in shaly formations (more than 800 samples from depths of 0–6 km) by extracting the major trends in the relation between seismic anisotropy and compaction. The statistical analysis of the database shows that the simultaneous increase in density and velocity, a classical compaction signature, is quite weakly correlated with the anisotropy strength. As a consequence, compaction can be excluded as a major cause of seismic anisotropy, at least in shaly formations. Also, the alignment of the clay platelets can explain most of the anisotropy measurements of both databases. Finally, a method for estimating the orientation distribution function of the clay platelets from the measurement of the anisotropy parameters is suggested.  相似文献   

4.
Transverse isotropy with a vertical axis of symmetry is a common form of anisotropy in sedimentary basins, and it has a significant influence on the seismic amplitude variation with offset. Although exact solutions and approximations of the PP-wave reflection coefficient for the transversely isotropic media with vertical axis of symmetry have been explicitly studied, it is difficult to apply these equations to amplitude inversion, because more than three parameters need to be estimated, and such an inverse problem is highly ill-posed. In this paper, we propose a seismic amplitude inversion method for the transversely isotropic media with a vertical axis of symmetry based on a modified approximation of the reflection coefficient. This new approximation consists of only three model parameters: attribute A, the impedance (vertical phase velocity multiplied by bulk density); attribute B, shear modulus proportional to an anellipticity parameter (Thomsen's parameter ε−δ); and attribute C, the approximate horizontal P-wave phase velocity, which can be well estimated by using a Bayesian-framework-based inversion method. Using numerical tests we show that the derived approximation has similar accuracy to the existing linear approximation and much higher accuracy than isotropic approximations, especially at large angles of incidence and for strong anisotropy. The new inversion method is validated by using both synthetic data and field seismic data. We show that the inverted attributes are robust for shale-gas reservoir characterization: the shale formation can be discriminated from surrounding formations by using the crossplot of the attributes A and C, and then the gas-bearing shale can be identified through the combination of the attributes A and B. We then propose a rock-physics-based method and a stepwise-inversion-based method to estimate the P-wave anisotropy parameter (Thomsen's parameter ε). The latter is more suitable when subsurface media are strongly heterogeneous. The stepwise inversion produces a stable and accurate Thomsen's parameter ε, which is proved by using both synthetic and field data.  相似文献   

5.
Yong Il  Lee  Dong Hyun  Lim 《Island Arc》2008,17(1):152-171
Abstract The Gyeongsang Basin is a non‐marine sedimentary basin formed by extensional tectonism during the Early Cretaceous in the southeastern Korean Peninsula. The sediment fill starts with the Sindong Group distributed along the western margin of the basin. It consists of three lithostratigraphic units: the Nakdong (alluvial fan), Hasandong (fluvial) and Jinju (lacustrine) formations with decreasing age. Sindong Group sandstones are classified into four petrofacies (PF) based on their detrital composition: PF‐A consists of the lower Nakdong Formation with average Q73F12R15; PF‐B the upper Nakdong and lower Hasandong formations with Q66F15R18; PF‐C the middle Hasandong to middle Jinju formations with Q49F29R22; and PF‐D the upper Jinju Formation with Q26F34R41. The variations of detrital composition influenced the diagenetic mineral assemblage in the Sindong Group sandstones. Illite and dolomite/ankerite are important diagenetic minerals in PF‐A and PF‐B, whereas calcite and chlorite are dominant diagenetic minerals in PF‐C and PF‐D. Most of the diagenetic minerals can be divided into early and late diagenetic stages of formation. Early diagenetic calcites occur mostly in PF‐C, probably controlled by arid to semiarid climatic conditions during the sandstone deposition, no early calcite being found in PF‐A and PF‐B. Late‐stage calcites are present in all Sindong Group sandstones. The calcium ions may have been derived from shale diagenesis and dissolution of early stage calcites in the Hasandong and Jinju sandstones. Illite, the only diagenetic clay mineral in PF‐A and lower PF‐B, is inferred to be a product of kaolinite transformation during deep burial, and the former presence of kaolinite is inferred from the humid paleoclimatic conditions during the deposition of the Nakdong Formation. Chlorites in PF‐C and PF‐D are interpreted to be the products of transformation of smectitic clay or of precipitation from alkaline pore water under arid to semiarid climatic conditions. The occurrence of late‐stage diagenetic minerals largely depended on the distribution of early diagenetic minerals, which was controlled initially by the sediment composition and paleoclimate.  相似文献   

6.
The simplified macro‐equations of porous elastic media are presented based on Hickey's theory upon ignoring effects of thermomechanical coupling and fluctuations of porosity and density induced by passing waves. The macro‐equations with definite physical parameters predict two types of compressional waves (P wave) and two types of shear waves (S wave). The first types of P and S waves, similar to the fast P wave and S wave in Biot's theory, propagate with fast velocity and have relatively weak dispersion and attenuation, while the second types of waves behave as diffusive modes due to their distinct dispersion and strong attenuation. The second S wave resulting from the bulk and shear viscous loss within pore fluid is slower than the second P wave but with strong attenuation at lower frequencies. Based on the simplified porous elastic equations, the effects of petrophysical parameters (permeability, porosity, coupling density and fluid viscosity) on the velocity dispersion and attenuation of P and S waves are studied in brine‐saturated sandstone compared with the results of Biot's theory. The results show that the dispersion and attenuation of P waves in simplified theory are stronger than those of Biot's theory and appear at slightly lower frequencies because of the existence of bulk and shear viscous loss within pore fluid. The properties of the first S wave are almost consistent with the S wave in Biot's theory, while the second S wave not included in Biot's theory even dies off around its source due to its extremely strong attenuation. The permeability and porosity have an obvious impact on the velocity dispersion and attenuation of both P and S waves. Higher permeabilities make the peaks of attenuation shift towards lower frequencies. Higher porosities correspond to higher dispersion and attenuation. Moreover, the inertial coupling between fluid and solid induces weak velocity dispersion and attenuation of both P and S waves at higher frequencies, whereas the fluid viscosity dominates the dispersion and attenuation in a macroscopic porous medium. Besides, the heavy oil sand is used to investigate the influence of high viscous fluid on the dispersion and attenuation of both P and S waves. The dispersion and attenuation in heavy oil sand are stronger than those in brine‐saturated sandstone due to the considerable shear viscosity of heavy oil. Seismic properties are strongly influenced by the fluid viscosity; thus, viscosity should be included in fluid properties to explain solid–fluid combination behaviour properly.  相似文献   

7.
  相似文献   

8.
饱和度对波在土层交界面的反射、透射系数的影响   总被引：3，自引：1，他引：3       下载免费PDF全文
把均匀流体概念引入Biot两相多孔介质动力理论中，用Biot的两相多孔介质模型模拟不完全饱和土层，给出SV波，P波从不完全饱和土层入射到弹性土层时，在土层交界面上反射，透射系数的表达式，结果表明与完全饱和相比，饱和度发生很小的变化就会对交界面上反射，透射系数产生很大的影响。今后应该重视饱和度变化对地震动力响应的影响。  相似文献   

9.
None of the standard porosity-velocity models (e.g. the time-average equation, Raymer's equations) is satisfactory for interpreting well-logging data over a broad depth range. Clays in the section are the usual source of the difficulty through the bias and scatter that they introduce into the relationship between porosity and P-wave transit time. Because clays are composed of fine sheet-like particles, they normally form pores with much smaller aspect ratios than those associated with sand grains. This difference in pore geometry provides the key to obtaining more consistent resistivity and sonic log interpretations. A velocity model for Clay–sand mixtures has been developed in terms of the Kuster and Toksöz, effective medium and Gassmann theories. In this model, the total pore space is assumed to consist of two parts: (1) pores associated with sand grains and (2) pores associated with clays (including bound water). The essential feature of the model is the assumption that the geometry of pores associated with sand grains is significantly different from that associated with clays. Because of this, porosity in shales affects elastic compliance differently from porosity in sand-Stones. The predictive power of the model is demonstrated by the agreement between its predictions and laboratory measurements and by its ability to predict sonic logs from other logs over large depth intervals where formations vary from unconsolidated to consolidated sandstones and shales.  相似文献   

10.
In this paper, we consider wave propagation in a layered medium. Using the Baker‐Campbell‐Hausdorff series, we expand the logarithm of a propagator matrix in series of frequency. The series coefficients allow us to extend the effective Backus medium for low frequencies. The proposed technique is applied to vertical propagation in a periodically layered and binary medium as well as for a gradient medium. The velocity dispersion equations are derived for these media. We also consider the layered medium with monoclinic anisotropy. We illustrate the accuracy of the proposed method on synthetic and well‐log data.  相似文献   

11.
The primary objective of this study was to evaluate use of the hydraulic profiling tool‐groundwater sampler (HPT‐GWS) log data as an indicator of water quality (level of dissolved ionic species) in an alluvial aquifer. The HPT‐GWS probe is designed for direct push advancement into unconsolidated formations. The system provides both injection pressure logs and electrical conductivity (EC) logs, and groundwater may be sampled at multiple depths as the probe is advanced (profiling). The combination of these three capabilities in one probe has not previously been available. During field work it was observed that when HPT corrected pressure (Pc) indicates a consistent aquifer unit then bulk formation EC can be used as an indicator of water quality. A high correlation coefficient (R 2 = 0.93) was observed between groundwater specific conductance and bulk formation EC in the sands and gravels of the alluvial aquifer studied. These results indicate that groundwater specific conductance is exerting a controlling influence on the bulk formation EC of the coarse‐grained unit at this site, and probably many similar sites, consistent with Archie's Law. This simple relationship enables the use of the EC and Pc logs, with targeted water samples and a minimum of core samples, to rapidly assess groundwater quality over extended areas at high vertical resolution. This method was used to identify both a brine impacted zone at the base of the aquifer investigated and a groundwater recharge lens developing below storm water holding ponds in the upper portion of the same aquifer. Sample results for trace level, naturally occurring elements (As, Ba, U) further demonstrate the use of this system to sample for low level groundwater contamination.  相似文献   

12.
In this paper a three‐dimensional dynamic analysis program for saturated porous rocks and soils is presented. The theoretical formulations incorporated in the proposed computer program are the extension of Biot's two‐phase theory to non‐linear region. Numerical study using the program is carried out to demonstrate the three‐dimensional features of the program and to show the effects of pore water on the dynamic response of underground openings in saturated rock masses. It is shown that underground openings in saturated porous media could be significantly more vulnerable to the potential damages associated with high motions and shear failure than those in dry media. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

13.
Abstract The Wakino Subgroup is a lower stratigraphic unit of the Lower Cretaceous Kanmon Group. Previous studies on provenance of Wakino sediments have mainly concentrated on either petrography of major framework grains or bulk rock geochemistry of shales. This study addresses the provenance of the Wakino sandstones by integrating the petrographic, bulk rock geochemistry, and mineral chemistry approaches. The proportions of framework grains of the Wakino sandstones suggest derivation from either a single geologically heterogeneous source terrane or multiple source areas. Major source lithologies are granitic rocks and high‐grade metamorphic rocks but notable amounts of detritus were also derived from felsic, intermediate and mafic volcanic rocks, older sedimentary rocks, and ophiolitic rocks. The heavy mineral assemblage include, in order of decreasing abundance: opaque minerals (ilmenite and magnetite with minor rutile), zircon, garnet, chromian spinel, aluminum silicate mineral (probably andalusite), rutile, epidote, tourmaline and pyroxene. Zircon morphology suggests its derivation from granitic rocks. Chemistry of chromian spinel indicates that the chromian spinel grains were derived from the ultramafic cumulate member of an ophiolite suite. Garnet and ilmenite chemistry suggests their derivation from metamorphic rocks of the epidote‐amphibolite to upper amphibolite facies though other source rocks cannot be discounted entirely. Major and trace element data for the Wakino sediments suggest their derivation from igneous and/or metamorphic rocks of felsic composition. The major element compositions suggest that the type of tectonic environment was of an active continental margin. The trace element data indicate that the sediments were derived from crustal rocks with a minor contribution from mantle‐derived rocks. The trace element data further suggest that recycled sedimentary rocks are not major contributors of detritus. It appears that the granitic and metamorphic rocks of the Precambrian Ryongnam Massif in South Korea were the major contributors of detritus to the Wakino basin. A minor but significant amount of detritus was derived from the basement rocks of the Akiyoshi and Sangun Terrane. The chromian spinel appears to have been derived from a missing terrane though the ultramafic rocks in the Ogcheon Belt cannot be discounted.  相似文献   

14.
Lower Cretaceous lacustrine oil shales are widely distributed in southeastern Mongolia. Due to the high organic carbon content of oil shale, many geochemical studies and petroleum exploration have been conducted. Although most of the oil shales are considered to be Early Cretaceous in age, a recent study reveals that some were deposited in the Middle Jurassic. The present study aims at establishing depositional ages and characteristics of the Jurassic and Cretaceous lacustrine deposits in Mongolia. The Lower Cretaceous Shinekhudag Formation is about 250 m thick and composed of alternating beds of shale and dolomite. The Middle Jurassic Eedemt Formation is about 150 m thick and composed of alternating beds of shale, dolomitic marl, and siltstone. The alternations of shale and dolomite in both formations were formed by lake level changes, reflecting precipitation changes. Shales were deposited in the center of a deep lake during highstand, while dolomites were formed by primary precipitation during lowstand. Based on the radiometric age dating, the Shinekhudag Formation was deposited between 123.8 ±2.0 Ma and 118.5 ±0.9 Ma of the early Aptian. The Eedemt Formation was deposited at around 165–158 Ma of Callovian–Oxfordian. The calculated sedimentation rate of the Shinekhudag Formation is between 4.7 ±2.6 cm/ky and 10.0 ±7.6 cm/ky. Shales in the Shinekhudag Formation show micrometer‐scale lamination, consisting of algal organic matter and detrital clay mineral couplets. Given the average thickness of micro‐laminae and calculated sedimentation rate, the micro‐lamination is most likely of varve origin. Both Middle–Upper Jurassic and Lower Cretaceous lacustrine oil shales were deposited in intracontinental basins in the paleo‐Asian continent. Tectonic processes and basin evolution basically controlled the deposition of these oil shales. In addition, enhanced precipitation under humid climate during the early Aptian and the Callovian–Oxfordian was another key factor inducing the widespread oil shale deposition in Mongolia.  相似文献   

15.
We derive the phase velocity dispersion and the scattering for wave vertically propagating in a periodically weak‐contrast horizontally layered medium with arbitrary number of layers in a period. Phase velocity dispersion is defined as the frequency dependence of vertical travel time, and scattering is defined as a reflection coefficient at the interface between the multilayered system and the corresponding Backus medium. Low‐frequency approximation is used to define a dynamic effective medium with frequency‐dependent phase velocity. The results are compared with those obtained earlier for a gradient medium. We show that the low‐frequency weak‐contrast approximation is valid for models with realistic contrasts in elastic properties.  相似文献   

16.
Multiparameter inversion for pre‐stack seismic data plays a significant role in quantitative estimation of subsurface petrophysical properties. However, it remains a complicated problem due to the non‐unique results and unstable nature of the processing; the pre‐stack seismic inversion problem is ill‐posed and band‐limited. Combining the full Zoeppritz equation and additional assumptions with edge‐preserving regularisation can help to alleviate these problems. To achieve this, we developed an inversion method by constructing a new objective function that includes edge‐preserving regularisation and soft constraints based on anisotropic Markov random fields and is intended especially for layered formations. We applied a fast simulated annealing algorithm to solve the nonlinear optimisation problem. The method directly obtains reflectivity RPP values using the full Zoeppritz equation instead of its approximations and effectively controls the stability of the multiparameter inversion by assuming a sectionally constant S‐ and P‐wave velocity ratio and using the generalised Gardner equation. We substituted the inverted parameters, i.e., the P‐wave velocity, the fitting deviation of S‐wave velocity, and the density were inverted instead of the P‐wave velocity, the S‐wave velocity, and the density, and the generalised Gardner equation was applied as a constraint. Test results on two‐dimensional synthetic data indicated that our substitution obtained improved results for multiparameter inversion. The inverted results could be improved by utilising high‐order anisotropic Markov random field neighbourhoods at early stages and low‐order anisotropic Markov random field neighbourhoods in the later stages. Moreover, for layered formations, using a large horizontal weighting coefficient can preserve the lateral continuity of layers, and using a small vertical weighting coefficient allows for large longitudinal gradients of the interlayers. The inverted results of the field data revealed more detailed information about the layers and matched the logging curves at the wells acceptably over most parts of the curves.  相似文献   

17.
马拉乌铅锌矿床为印度尼西亚西加里曼丹（Kalimantan）省吉打邦县近年来发现的铅锌银矿床。该矿矿体顶板为黄铁矿化含碳质泥砂岩,矿体底板为玉髓状硅质岩。矿体呈侧羽状、脉状及透镜状。矿区通常发育的蚀变有硅化、泥化和绿泥石化等。主要矿物组合有方铅矿、闪锌矿、黄铁矿、磁黄铁矿、黄铜矿、斑铜矿、辉锑矿以及自然金、银金矿、辉银矿等,通过对矿石组构等综合分析,表明该矿床的成因是岩浆期后热液充填交代矿床。  相似文献   

18.
Existing and commonly used in industry nowadays, closed‐form approximations for a P‐wave reflection coefficient in transversely isotropic media are restricted to cases of a vertical and a horizontal transverse isotropy. However, field observations confirm the widespread presence of rock beds and fracture sets tilted with respect to a reflection boundary. These situations can be described by means of the transverse isotropy with an arbitrary orientation of the symmetry axis, known as tilted transversely isotropic media. In order to study the influence of the anisotropy parameters and the orientation of the symmetry axis on P‐wave reflection amplitudes, a linearised 3D P‐wave reflection coefficient at a planar weak‐contrast interface separating two weakly anisotropic tilted tranversely isotropic half‐spaces is derived. The approximation is a function of the incidence phase angle, the anisotropy parameters, and symmetry axes tilt and azimuth angles in both media above and below the interface. The expression takes the form of the well‐known amplitude‐versus‐offset “Shuey‐type” equation and confirms that the influence of the tilt and the azimuth of the symmetry axis on the P‐wave reflection coefficient even for a weakly anisotropic medium is strong and cannot be neglected. There are no assumptions made on the symmetry‐axis orientation angles in both half‐spaces above and below the interface. The proposed approximation can be used for inversion for the model parameters, including the orientation of the symmetry axes. Obtained amplitude‐versus‐offset attributes converge to well‐known approximations for vertical and horizontal transverse isotropic media derived by Rüger in corresponding limits. Comparison with numerical solution demonstrates good accuracy.  相似文献   

19.
The important effects of bottom sediments on the seismic response of arch dams are studied in this paper. To do so, a three‐dimensional boundary element model is used. It includes the water reservoir as a compressible fluid, the dam and unbounded foundation rock as viscoelastic solids, and the bottom sediment as a two‐phase poroelastic domain with dynamic behaviour described by Biot's equations. Dynamic interaction among all those regions, local topography and travelling wave effects are taken into account. The results obtained show the important influence of sediment compressibility and permeability on the seismic response. The former is associated with a general change of the system response whereas the permeability has a significant influence on damping at resonance peaks. The analysis is carried out in the frequency domain considering time harmonic excitation due to P and S plane waves. The time‐domain results obtained by using the Fourier transform for a given earthquake accelerogram are also shown. The possibility of using simplified models to represent the bottom sediment effects is discussed in the paper. Two alternative models for porous sediment are tested. Simplified models are shown to be able to reproduce the effects of porous sediments except for very high permeability values. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

20.
An equivalent medium model for wave simulation in fractured porous rocks   总被引：3，自引：0，他引：3  
Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.  相似文献   

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

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