Integrated Two-Dimensional Modeling of Fluid Flow and Compaction in a Sedimentary Basin     

Zhangxin Chen  Richard E. Ewing  Hao Lu  Stephen L. Lyons  Serguei Maliassov  Michael B. Ray  Tong Sun 《Computational Geosciences》2002,6(3-4):545-564

In this paper we employ mixed finite elements and numerically study an integrated two-dimensional model of fluid flow and compaction in a sedimentary basin. This model describes a single phase incompressible flow in a two-dimensional section of a sedimentary basin with vertical compaction. At each time step, an iterative algorithm is used to solve this model. The determination of the grid movement is based on the mass conservation and movement of sediments in the basin, while the mixed method is utilized to solve the fluid flow over the moving grid. Numerical experiments are presented to verify this iterative algorithm and show representative solutions for the model under consideration.  相似文献   

湖南印支期过铝质花岗岩的形成: 岩浆底侵与地壳加厚热效应的数值模拟   总被引：36，自引：0，他引：36  

王岳军  Y.H.ZHANG  范蔚茗  席先武  郭锋  林舸 《中国科学D辑》2002,32(6):491-499

在合理构建华南印支期地质-物理模型的基础上, 利用FLAC软件, 模拟了该区印支期过铝质富钾花岗岩形成与基性岩浆底侵, 及陆壳变形叠置加厚两种动力学背景的可能联系. 模拟结果表明, 220 Ma±的基性岩浆底侵能导致地壳含水矿物相岩石的深熔, 但除非印支期存在大规模基性岩浆的底侵作用, 否则其热效应持续时间和热效应波及范围难以形成具大岩基规模的湖南印支期花岗岩. 在陆壳叠置加厚模型中, 地壳的叠置加厚可导致中下地壳界面温度升高到700℃以上, 引起片麻质岩石熔融, 当加厚因子达1.3, 白云母矿物脱水熔融产生的熔体达到熔体流动临界比例(≥20%), 从而形成花岗岩基; 结合印支期挤压逆冲推覆构造和同期基性火山岩极少出露的地质事实, 认为陆壳变形加厚可能是湖南印支期构造岩浆作用形成的主导因素.  相似文献   

Improved NMO correction with a specific application to shallow-seismic data   总被引：7，自引：0，他引：7  

Jan H. Brouwer 《Geophysical Prospecting》2002,50(2):225-237

The application of traditional NMO‐correction techniques in the processing of seismic data may result in severe distortion. This distortion is observed as a decrease in frequency content (NMO stretch) or even a total loss of data for steep velocity gradients. It is shown that the specific nature of the observed distortion is introduced by the particular implementation of the traditional NMO‐correction technique. It is also shown that other techniques that have been proposed in the literature suffer from similar artefacts. An alternative approach is suggested, based on the correction of tapered blocks of seismic data, followed by a coherence filter to compensate for the specific artefacts thus introduced. Correction of seismic amplitudes (e.g. geometrical spreading and attenuation) is implemented as an integral part of the NMO‐correction method thus allowing for both dynamic and kinematic reconstruction of interfering events. From the application of this method to synthetic and field data, it is concluded that the proposed technique may be particularly useful in the processing of shallow seismic (shear‐wave) data. Although not illustrated in this paper, it is emphasized that application of the proposed method is by no means restricted to shear waves; application to P‐wave data will prove equally useful, although the stretch effects involved may not be as severe as in the examples given in this paper.  相似文献   

Analysis of the resolution function in seismic prestack depth imaging   总被引：1，自引：0，他引：1  

L.-J. Gelius  I. Lecomte  H. Tabti 《Geophysical Prospecting》2002,50(5):505-515

We consider the problem of estimating subsurface quantities such as velocity or reflectivity from seismic measurements. Because of a limited aperture and band-limited signals, the output from a seismic prestack reconstruction method is a distorted or blurred image. This distortion can be computed using the concept of resolution function, which is a quantity readily accessible in the Fourier space of the model. The key parameter is the scattering wavenumber, which at a particular image point is defined by the incident and scattered ray directions in a given background model. Any location in any background model can be considered. In general, the resolution function will depend on the following four quantities: the background velocity model, the frequency bandwidth, the wavefield type and the acquisition geometry.
We first establish the resolution function for a general scattering model assuming local reaction. We then adapt this result for two well-known scattering models: Born and Kirchhoff. For each of these approximations the corresponding resolution function is derived and discussed. Finally, by employing a simple synthetic data example we demonstrate the ability of the resolution function to predict the image distortions.  相似文献   

Wave velocities and attenuation of shaley sandstones as a function of pore pressure and partial saturation   总被引：3，自引：0，他引：3  

Nam H. Pham  †  José M. Carcione  Hans B. Helle  Bjørn Ursin 《Geophysical Prospecting》2002,50(6):615-627

We obtain the wave velocities and quality factors of clay‐bearing sandstones as a function of pore pressure, frequency and partial saturation. The model is based on a Biot‐type three‐phase theory that considers the coexistence of two solids (sand grains and clay particles) and a fluid mixture. Additional attenuation is described with the constant‐Q model and viscodynamic functions to model the high‐frequency behaviour. We apply a uniform gas/fluid mixing law that satisfies the Wood and Voigt averages at low and high frequencies, respectively. Pressure effects are accounted for by using an effective stress law. By fitting a permeability model of the Kozeny– Carman type to core data, the model is able to predict wave velocity and attenuation from seismic to ultrasonic frequencies, including the effects of partial saturation. Testing of the model with laboratory data shows good agreement between predictions and measurements.  相似文献   

Transient foundation solution of saturated soil to impulsive concentrated loading     

X. L. Zhou  J. H. Wang  J. F. Lu 《Soil Dynamics and Earthquake Engineering》2002,22(4)

The transient dynamic response of saturated soil under suddenly applied normal and horizontal concentrated loading is studied in this paper. The behavior of saturated soil is governed by Biot's consolidation theory. The general solutions for Biot equations of equilibrium are derived in terms of displacements and variations of fluid volume, using Laplace–Hankel integral transforms. The solutions in the time domain can be evaluated by numerical inverse Laplace–Hankel transforms. Selected numerical results for displacements, stresses, and pore pressures are presented. Comparisons with existing closed-form solutions for the elastic half-space are made to confirm the accuracy of the present solutions. The solutions can be used to study a variety of transient wave propagation problems and dynamical interactions between saturated soil and structures.  相似文献   

The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk     

J. J. Bommer  M. B. Benito  M. Ciudad-Real  A. Lemoine  M. A. Lpez-Menjívar  R. Madariaga  J. Mankelow  P. Mndez de Hasbun  W. Murphy  M. Nieto-Lovo  C. E. Rodríguez-Pineda  H. Rosa 《Soil Dynamics and Earthquake Engineering》2002,22(5)

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.  相似文献   

Fractal properties of medium and seismoelectric phenomena     

V. V. Surkov  S. Uyeda  H. Tanaka  M. Hayakawa 《Journal of Geodynamics》2002,33(4-5)

Electrokinetic phenomena in a water-porous medium with a fractal structure above percolation threshold are theoretically investigated. Fracture zone with space-variable porosity is considered as a model of an earthquake hypocenter zone in which the electrokinetic current results from fluid filtration in a fractal pore network. A critical exponent of the streaming potential coefficient is found to depend on both the transport critical exponent and correlation length critical exponent. In this model, logarithmic dependence of electric field amplitude E on the earthquake magnitude M is derived which is compatible with the one observed by the VAN group. Without fractal properties, this form of dependence contradicts the empirical data. The electromagnetic field far from the hypocenter is calculated, which leads to the prediction of weak magnetic field variations. To explain the observed amplitude of VAN's Seismic Electric Signals (SES), the electric source must be at a distance of about 10 km from the registration point if the medium is homogeneous. Therefore, some conductive channel(s) are needed to explain the long distance selective SES transmission.  相似文献   

Determining Q of near-surface materials from Rayleigh waves     

Jianghai Xia  Richard D. Miller  Choon B. Park  Gang Tian 《Journal of Applied Geophysics》2002,51(2-4)

High-frequency (≥2 Hz) Rayleigh wave phase velocities can be inverted to shear (S)-wave velocities for a layered earth model up to 30 m below the ground surface in many settings. Given S-wave velocity (V_S), compressional (P)-wave velocity (V_P), and Rayleigh wave phase velocities, it is feasible to solve for P-wave quality factor Q_P and S-wave quality factor Q_S in a layered earth model by inverting Rayleigh wave attenuation coefficients. Model results demonstrate the plausibility of inverting Q_S from Rayleigh wave attenuation coefficients. Contributions to the Rayleigh wave attenuation coefficients from Q_P cannot be ignored when Vs/V_P reaches 0.45, which is not uncommon in near-surface settings. It is possible to invert Q_P from Rayleigh wave attenuation coefficients in some geological setting, a concept that differs from the common perception that Rayleigh wave attenuation coefficients are always far less sensitive to Q_P than to Q_S. Sixty-channel surface wave data were acquired in an Arizona desert. For a 10-layer model with a thickness of over 20 m, the data were first inverted to obtain S-wave velocities by the multichannel analysis of surface waves (MASW) method and then quality factors were determined by inverting attenuation coefficients.  相似文献   

Paleolatitudes of the Kerguelen hotspot: new paleomagnetic results and dynamic modeling   总被引：1，自引：0，他引：1  

M. Antretter  B. Steinberger  H. Soffel 《Earth and Planetary Science Letters》2002,203(2):635-650

The Kerguelen Plateau, a Large Igneous Province in the southern Indian Ocean, was formed as a product of the Kerguelen hotspot in several eruptive phases during the last 120 Myr. We obtained new paleolatitudes for the central and northern Kerguelen Plateau from paleomagnetic investigations on basalts, which were drilled during ODP Leg 183 to the Kerguelen Plateau-Broken Ridge. The paleolatitudes coincide with paleolatitudes from previous investigations at the Kerguelen Plateau and Ninetyeast Ridge (the track of the Kerguelen hotspot) and indicate a difference between paleolatitudes and present position at 49°S of the Kerguelen hotspot. We show that true polar wander, the global motion between the mantle and the rotation axis, cannot explain this difference in latitudes. We present numerical model results of plume conduit motion in a large-scale mantle flow and the resulting surface hotspot motion. A large number of models all predict southward motion between 3° and 10° for the Kerguelen hotspot during the last 100 Myr, which is consistent with our paleomagnetic results.  相似文献