首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 406 毫秒
1.
The generalised three layer boundary value problem with a transition layer sand-witched between an isotropic overburden and dipping anisotropic substratum is discussed assuming that plane electro-magnetic waves are incident normally over the air-earth inter-face. The tangential electric (Ey) and magnetic (Hx) fields and the expression for surface-impedance (Ey/Hx) have been evaluated at the earth's surface. Through numerical analysis it is shown that changes in the values of the parameters m (coefficient of anisotropy), h (thickness of the transition layer), α (angle of inclination of the dipping beds), and b (conductivity ratio between substratum and upper layer) modify the amplitude and phase-variation curves (with skindepth) significantly.  相似文献   

2.
The detectability of an intermediate layer in a three-layer earth model in the time domain has been investigated. The calculations were made for the perpendicular loop (designated system II) and vertical-coplanar (designated system III) electromagnetic (EM) sounding systems. The primary excitation employed is a train of half-sinusoidal and square waveforms of alternating polarity. The time-domain response has been determined by Fourier transformation of the matched complex mutual coupling ratios into the time domain and by linear digital filtering. Top and bottom layers have equal resistivity. EM responses have been computed for conductive and resistive intermediate layer with a wide range of thickness and for two values (500 m and 1000 m) of loop-separation. For the detectability analyses, the root mean square (rms) difference between three-layer and homogeneous-earth responses is adapted. The threshold value for detectability is defined as an rms difference of 10% and the measurement error is arbitrarily assumed to be of the order of 3%. It is observed that the perpendicular-loop system is better than the vertical-coplanar system in detecting thin intermediate layers (either conductive or resistive). For a loop separation of 1000 m and half-sinusoidal pulse excitation, the detectable thickness ratio (h2/h1) is 0.10 by system II for the conducting middle layers; for square pulse excitation the corresponding thickness ratios are 0.06 for system II and 0.12 for system III. For a loop separation of 1000 m and half-sinusoidal pulse excitation in detecting the resistive intermediate layers, the corresponding thickness ratios are 0.9 for system II and 2.25 for system III; while for square pulse excitation the thickness ratios are 0.55 for system II and 1.55 for system III. Results in the frequency domain and time domain (for half-sinusoidal and square pulsed field) have also been presented for systems II and III for detecting conducting layers by considering an earth model where p1≠ p3 and p3 > p1 (p is the resistivity). The loop separa- tion used is 1000 m. Direct comparisons between the frequency domain and time-domain results clearly demonstrate the superiority of frequency-domain systems for detecting con- ducting intermediate layers.  相似文献   

3.
Abstract

A cylindrical annulus containing a conducting fluid and rapidly rotating about its axis is a useful model for the Earth's core. With a shear flow U 0(s)∮, magnetic field B 0(s)∮, and temperature distribution T o(s) (where (s, ∮, z) are cylindrical polar coordinates), many important properties of the core can be modelled while a certain degree of mathematical simplicity is maintained. In the limit of rapid rotation and at geophysically interesting field strengths, the effects of viscous diffusion and fluid inertia are neglected. In this paper, the linear stability of the above basic state to instabilities driven by gradients of B 0 and U 0 is investigated. The global numerical results show both instabilities predicted by a local analysis due to Acheson (1972, 1973, 1984) as well as a new resistive magnetic instability. For the non-diffusive field gradient instability we looked at both monotonic fields [for which the local stability parameter Δ, defined in (1.4), is a constant] and non-monotonic fields (for which Δ is a function of s). For both cases we found excellent qualitative agreement between the numerical and local results but found the local criterion (1.6) for instability to be slightly too stringent. For the non-monotonic fields, instability is confined approximately to the region which is locally unstable. We also investigated the diffusive buoyancy catalysed instability for monotonic fields and found good quantitative agreement between the numerical results and the local condition (1.9). The new resistive instability was found for fields vanishing (or small) at the outer boundary and it is concentrated in the region of that boundary. The resistive boundary layer plays an important part in this instability so it is not of a type which could be predicted using a local stability analysis (which takes no account of the presence of boundaries).  相似文献   

4.
Magnetotelluric response is studied for an inhomogeneous medium having conductivity varying linearly with depth as σ(z) =σ1+αz. For a medium having conductivity increasing linearly with depth, the phase of the impedance approaches 60° at long periods and the apparent resistivity becomes log (ρa) = 2 log (1.36/α1/3) — 1/3 log (T'). The asymptote of log (ρa, T'→∞) when plotted against log (T') has a constant gradient —1/3 and has an intercept on the log (T') axis, which equals 6 log (1.36/α1/3). When a homogeneous layer with a moderate thickness overlies an inhomogeneous half-space, this layer does not affect the asymptote, but it affects the cut-off period and pushes this toward the long period direction. For a medium having conductivity decreasing linearly with depth, the impedance is equivalent to that of a Cagniard two-layer model; the intercept period related to the thickness is T'0=σ1(h2/2)2. Homogeneous multilayer approximations to an inhomogeneous layer are also investigated, and it is shown that the fit to the model variation depends on the number of layers and the layer parameters chosen.  相似文献   

5.
Hydrocarbon source systems and formation of gas fields in Sichuan Basin   总被引：1，自引：0，他引：1  
The formation of large and middle gas fields in Sichuan Basin is investigated based on source wntrolling theory and hydrocarbon source systems. It is indicated that Є1, Sl, P1, P2 and T3 are the main source beds and Є1/Z2 d, C2 h/S1, P1/P2, P2 ch/P2, T1,2/P,T3 x /T3 x are important hydrocarbon source systems in the basin. All these source systems are the prospective formations and exploration spaces of large and middle gas fields. It is also emphasized that hydrocarbon generation intensity is the most important geochemical factor to estimate large and middle gas fields. Project supported by the “85–102” Chinese National Key Science and Technology Project.  相似文献   

6.
This paper discusses the asymptotic behaviour of the electromagnetic fields received on the sea‐bed (target response), as well as the fields distributed inside a thin resistive target, generated by a horizontal electric dipole above the sea‐bed in marine controlled‐source electromagnetics for hydrocarbon exploration. It is found that the guided wave supported by a thin resistive target can be expressed as a single‐mode exponential function. A simple closed‐form expression is derived to relate the single‐mode wavenumber of the guided wave to the model parameters: the resistivity and thickness of the target layer, the sea‐bed resistivity and the frequency. When the air‐wave is removed, the guided wave is dominant among the fields received on the sea‐bed at far offset. Hence the wavenumber of the guided wave can be calculated from the fields measured on the sea‐bed. The closed‐form expression can then be used to invert the target property from the calculated wavenumber and hence, can be considered as a hydrocarbon indicator.  相似文献   

7.
Refraction along thin high velocity layers and along basement is investigated in two cases. a) high velocity layer just on the basement. b) high velocity layer higher above. Period and attenuation of refracted waves are givers as a function of the layer thickness H. Refracted arrivals along thin high velocity layers are visible at significant distances if the layer thickness is not smaller than A/6, where A is the longitudinal wavelength in high velocity medium. The pseudoperiod is proportional to the layer thickness H. The attenuation at large distance follows an x-ne-k1x law, where n is close to I and k1 is inversely proportional to H. Refracted arrivals along the basement are observable even in the case of thin high velocity layers situated in the overburden; their intensity is smaller and their pseudo-period larger than when no layer exists in the overburden. The intensity of the basement arrival decreases and the pseudoperiod increases with increasing layer thickness. The pseudoperiod and the attenuation of refracted arrivals along high velocity layers and along the basement are also highly dependent on acoustic contrasts. Both arrivals from a high velocity layer and from the basement can be recorded simultaneously, provided the frequency spectrum of the seismic chain is sufficiently broad. In all cases layer arrivals show a character very different from basement arrivals.  相似文献   

8.
Proposed is a new definition of earthquake response spectra, which takes account of the number of response cycles N. The Nth largest amplitude of absolute acceleration response of a linear oscilator with natural period T and damping ratio h, which is subjected to ground motion at its base, is defined as SA(T, h, N). By defining a reduction factor η(T, h, N) as SA(T, h, N)/SA(T, h, 1), characteristics of η(T, h, N) were investigated based on 394 components of strong motion records obtained in Japan. Two practical empirical formulae to assess the reduction factor η(T, h, N) are proposed.  相似文献   

9.
Computations of the time-domain electromagnetic response of a multi-layered earth have been carried out for different source-receiver coil systems. The primary excitation is a train of half-sinusoidal waveforms of alternating polarity. The conversion into the time-domain involves Fourier series summation of the matched complex mutual coupling ratios of the layered earth models computed by a digital linear filter method. As an example, the response of a perpendicular coil system on the ground surface for two source-receiver separations has been presented for a five-layer earth model. This has been compared with the responses of homogeneous, two-layer, three-layer, and four-layer models. Next, the investigations have been extended to study the problems of equivalence of three-layer models, the intermediate layer of which is either conductive or resistive. For an intermediate conductive layer (H-type), the studies show that in the early portion of the signal the interpretation of a true three-layer earth is possible to some extent, whereas the ambiguity due to equivalence persists in the late samples. On the other hand, for an intermediate resistive layer (K-type), the three-layer earth and its equivalent model cannot be distinguished from each other over the entire sampling period. On the basis of a computational approach, equivalence has been empirically established as √h/ρ=constant for H-type earth-sections, and as h2ρ=constant for K-type earth sections, where h and ρ are respectively the thickness and resistivity of the intermediate layer.  相似文献   

10.
To investigate the characteristics of the anisotropic stratum, a multi‐azimuth seismic refraction technique is proposed in this study since the travel time anomaly of the refraction wave induced by this anisotropic stratum will be large for a far offset receiver. To simplify the problem, a two‐layer (isotropy–horizontal transverse isotropy) model is considered. A new travel time equation of the refracted P‐wave propagation in this two‐layer model is derived, which is the function of the phase and group velocities of the horizontal transverse isotropic stratum. In addition, the measured refraction wave velocity in the physical model experiment is the group velocity. The isotropic intercept time equation of a refraction wave can be directly used to estimate the thickness of the top (isotropic) layer of the two‐layer model because the contrast between the phase and group velocities of the horizontal transverse isotropic medium is seldom greater than 10% in the Earth. If the contrast between the phase and group velocities of an anisotropic medium is small, the approximated travel time equation of a refraction wave is obtained. This equation is only dependent on the group velocity of the horizontal transverse isotropic stratum. The elastic constants A11, A13, and A33 and the Thomsen anisotropic parameter ε of the horizontal transverse isotropic stratum can be estimated using this multi‐azimuth seismic refraction technique. Furthermore, under a condition of weak anisotropy, the Thomsen anisotropic parameter δ of the horizontal transverse isotropic stratum can be estimated by this technique as well.  相似文献   

11.
We analysed in detail three earthquakes recorded in a small-aperture accelerometric array in Mexico City, using the correlation of the records as a function of time along the accelerogram and frequency. Ground response is strongly conditioned by the fundamental period of the soft soils at the site of the array (T0). Energy at periods longer than 2T0 is guided by the crustal structure (with a thickness of 45 km). The wave field at periods between T0 and 2T0 also consists of surface waves but guided by the upper 2–3 km of volcanic sediments in central Mexico. For periods smaller than T0, ground motion is uncorrelated among the stations. Our results indicate that seismic response of Mexico City, including its very long duration, results from deeply guided surface waves (between 2 and 45 km depth) interacting with the very local response of the soft surficial clay layer.  相似文献   

12.
This proposed technique allows sensible and numerically stable behavior in multilayer analytic element models when layers dewater. When saturated thickness approaches zero in an unconfined or fresh/salt interface domain, the domain transitions to a very thin confined domain with a minimum saturated thickness M. M is an adjustable input parameter, so you can make the horizontal flow in dewatered domains negligibly small by making the minimum saturated thickness very small. Vertical flows can pass through a dewatered domain, whether it is near the surface or at depth. For example, recharge may pass through a shallow dewatered layer to a deeper layer that is not dewatered. This approach is examined in detail in an example multilayer model of mine dewatering.  相似文献   

13.
Summary TheWiechert-Herglotz method to investigate seismically the Earth structure, is adapted for refraction prospecting. The case of a single layer bounded by a sloping interface is considered. The velocity is assumed to be an increasing function of depth within the layer and constant beneath. A general procedure of computation is outlined. Special formulae are given in addition for the case in which the dromochroneD 1 corresponding to disturbances propagated only through the layer are of parabolic shape. Table I was calculated using them and assuming an horizontal interface. It may be seen there that valuesh 1 of the layer thickness furnished by these formulae may become considerably greater than those based on D1 rectilinear. Computation ofh 1 may easily be performed by means of the same table.  相似文献   

14.
For the computation of the vertical component Hz of the magnetic field of a horizontal A.C. dipole lying on the earth's surface, a recurrence formula is presented for a horizontally stratified half space, to obtain the (n+ 1)-layer case from the w-layer case. By means of several computed diagrams for the two-layer case, Hz can be determined for different ratios of conductivity of the subsoil and that of the overburden. Thereby the distance from the dipole as well as the layer thickness h are expressed in terms of the wave length A of a plain wave in the overburden. Assuming a sufficiently large conductivity difference, the results show that evidence about the subsurface conditions can be obtained if the distance between the measuring coil and the dipole is of the order of A/3, and if the thickness h of the layer varies within the range A/100 < h < A/6. As an example for the 3-layer case, a nonconducting intermediate layer is assumed.  相似文献   

15.
We present a numerical study for 3D time‐lapse electromagnetic monitoring of a fictitious CO2 sequestration using the geometry of a real geological site and a suite of suitable electromagnetic methods with different source/receiver configurations and different sensitivity patterns. All available geological information is processed and directly implemented into the computational domain, which is discretized by unstructured tetrahedral grids. We thus demonstrate the performance capability of our numerical simulation techniques. The scenario considers a CO2 injection in approximately 1100 m depth. The expected changes in conductivity were inferred from preceding laboratory measurements. A resistive anomaly is caused within the conductive brines of the undisturbed reservoir horizon. The resistive nature of the anomaly is enhanced by the CO2 dissolution regime, which prevails in the high‐salinity environment. Due to the physicochemical properties of CO2, the affected portion of the subsurface is laterally widespread but very thin. We combine controlled‐source electromagnetics, borehole transient electromagnetics, and the direct‐current resistivity method to perform a virtual experiment with the aim of scrutinizing a set of source/receiver configurations with respect to coverage, resolution, and detectability of the anomalous CO2 plume prior to the field survey. Our simulation studies are carried out using the 3D codes developed in our working group. They are all based on linear and higher order Lagrange and Nédélec finite‐element formulations on unstructured grids, providing the necessary flexibility with respect to the complex real‐world geometry. We provide different strategies for addressing the accuracy of numerical simulations in the case of arbitrary structures. The presented computations demonstrate the expected great advantage of positioning transmitters or receivers close to the target. For direct‐current geoelectrics, 50% change in electric potential may be detected even at the Earth's surface. Monitoring with inductive methods is also promising. For a well‐positioned surface transmitter, more than 10% difference in the vertical electric field is predicted for a receiver located 200 m above the target. Our borehole transient electromagnetics results demonstrate that traditional transient electromagnetics with a vertical magnetic dipole source is not well suited for monitoring a thin horizontal resistive target. This is due to the mainly horizontal current system, which is induced by a vertical magnetic dipole.  相似文献   

16.
A non-hydrostatic model in cross-sectional form with an idealized sill is used to examine the influence of sill depth (h s) and aspect ratio upon internal motion. The model is forced with a barotropic tide and internal waves and mixing occurs at the sill. Calculations using a wide sill and quantifying the response using power spectra show that for a given tidal forcing namely Froude number F r as the sill depth (h s) increases the lee wave response and vertical mixing decrease. This is because of a reduction in across sill velocity U s due to increased depth. Calculations show that the sill Froude number F s based on sill depth and across sill velocity is one parameter that controls the response at the sill. At low F s (namely F s ≪ 1) in the wide sill case, there is little lee wave production, and the response is in terms of internal tides. At high F s, calculations with a narrow sill show that for a given F s value, the lee wave response and internal mixing increase with increasing aspect ratio. Calculations using a narrow sill with constant U s show that for small values of h s, a near surface mixed layer can occur on the downstream side of the sill. For large values of h s, a thick well-mixed bottom boundary layer occurs due to turbulence produced by the lee waves at the seabed. For intermediate values of h s, “internal mixing” dominates the solution and controls across thermocline mixing.  相似文献   

17.
We study two asymptotic regimes of unstable miscible displacements in porous media, in the two limits, where a permeability-modified aspect ratio, RL=L/H(kv/kh)1/2, becomes large or small, respectively. The first limit is known as transverse (or vertical) equilibrium, the second leads to the problem of non-communicating layers (the Dykstra–Parsons problem). In either case, the problem reduces to the solution of a single integro-differential equation. Although at opposite limits of the parameter RL, the two regimes coincide in the case of equal viscosities, M=1. By comparison with high-resolution simulation we investigate the validity of these two approximations. The evolution of transverse averages, particularly under viscous fingering conditions, depends on RL. We investigate the development of a model to describe viscous fingering in weakly heterogeneous porous media under transverse equilibrium conditions, and compare with the various existing empirical models (such as the Koval, Todd–Longstaff and Fayers models).  相似文献   

18.
The effects of variable speeds of spreading of submarine slides and slumps on near-field tsunami amplitudes are illustrated. It is shown that kinematic models of submarine slides and slumps must consider time variations in the spreading velocities, when these velocities are less than about 2cT, where is the long period tsunami velocity in ocean of constant depth h. For average spreading velocities greater than 2cT, kinematic models with assumed constant spreading velocities provide good approximation for the tsunami amplitudes above the source.  相似文献   

19.
We have determined the elastic thicknessTe of the oceanic lithosphere along two volcanic chains of the South Central Pacific: Cook-Austral and Society islands. We used a three-dimensional spatial method to model the lithospheric flexure assuming a continuous elastic plate. The model was constrained by geoid height data from the SEASAT satellite.Along the Cook-Austral chain the elastic thickness increases westward, from 2–4 km at McDonald hot spot to 14 km at Rarotonga. At McDonald seamount, however, the data are better explained by a local compensation model. The observed trend shows an increase ofTe with age of plate at loading time. However, the elastic layer under the Cook-Austral appears systematically thinner by several kilometers than expected for “normal” seafloor, suggesting that substantial thermal thinning has taken place in this region. Considering the apparent thermal age of the plate instead of crustal age improves noticeably the results. Along the Society chainTe varies from 20 km under Tahiti to 13 km under Maupiti which is located 500 km westward. When plotting together the Society and Cook-AustralTe results versus age of load, we notice that within the first five million years after loading,Te decreases significantly while tending rapidly to an equilibrium value. This may be interpreted as the effect of initial stress relaxation which occurs just after loading inside the lower lithosphere and suggests that the presently measured elastic thickness under the very young Tahiti load ( 0.8 Ma) is not yet the equilibrium thickness.  相似文献   

20.
  相似文献   

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

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