首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
At Delft Geotechnics the technique of ground-penetrating radar is in use for the detection of buried objects such as pipes. In order to give our ‘measurements in the field’ a more quantitative interpretation, a series of experiments has been started under well-defined conditions. A cylindrical vessel contains water, simulating wet soil. A pulsed transmitting antenna (TA) is mounted above the water-surface irradiating horizontal underwater pipes. The reflected pulses are detected by a horizontal receiving dipole (AP, i.e. ‘air-probe’) in the air. The reflecting objects used in the experiments are steel pipes, PVC rods and hollow PVC pipes filled with air or water. The depth of the pipes varies from 0.25 to 2 m. The strength of the reflected pulse depends on the type of pipe, its diameter, its depth, the electromagnetic properties of the water and also on the strength and polarization of the incident E-field. The latter is (mainly) parallel to the axis of the pipe in the present experiments. The experimental results have been compared with calculated results using the theory of a continuous plane wave, incident on an infinitely long cylinder in a homogeneous dielectric medium with the same dielectric constant as water. In a previous paper an experiment was described in which a movable receiving dipole in the water measured the transverse, mutually orthogonal Eø- and Eθ-components of the pulses emitted by the TA. The amplitudes of Eθ versus depth, measured in that paper, are used in the calculations as amplitudes of the incident field. The attenuation of the scattered field is accounted for by assuming exponential damping. The relative amplitudes of the scattered field for different pipes lying at the same depth are reproduced very well by this simple theory. Also the relative values of the scattered amplitudes for a given pipe lying at different depths are described neatly, provided the refraction of the scattered field at the water-air interface is accounted for.  相似文献   

2.
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.  相似文献   

3.
The effective relative dielectric constant ?e, r and the effective conductivity σe have each been determined as a function of frequency in the range 1–3000 MHz at volumetric water contents of up to approximately 0.74 for clays, 0.83 for a peat and 0.56 for a silt. At frequencies above about 25 MHz (depending on soil type), ?e, rincreases with water content for all samples. However, at lower frequencies, ?e, ronly increases with water content as long as the wet density also increases, which is the case for water contents up to a critical value lying between 0.35 and 0.48. At higher water contents, ?e, rand the wet density decrease with increasing water content. Consequently, curves of ?e, rversus frequency for two wet samples with different water contents, at least one of them higher than the critical value, are seen to cross at about 25 MHz. Below the critical value the curve of the sample with the lower water content is below the other curve at all freqencies applied. At a given frequency, σe has a maximum as a function of water content. This is tentatively explained by assuming that σe is the sum of pore water conductivity (increasing with water content until all salts in the soil are dissolved into the water and then decreasing) and surface water conductivity (increasing with wet density and therefore increasing with water content up to the critical value and then decreasing). At frequencies higher than 1000 MHz, ?e, rdepends only weakly on salinity (which is represented by the measured conductivity). It shows an increasing dependence if the frequency is decreased towards 1 MHz. The highest values of ?e, rand σe, measured in this work, occur for a sample of wet, nearly saturated silt originating from a location below sea-level near to the Dead Sea, Israel: ?e, rdecreases continuously from a value of about 104 at 3 MHz to about 102 at 200 MHz, while σe rises from about 4 S/m to 5 S/m at these respective frequencies. The dependence of the wavelength on the loss-tangent is strong and the wavelength is considerably smaller than it would be in a dielectric. This is the only sample for which σe increases with water content, even if the latter is above its critical value. Therefore it is assumed that the pore water conductivity is greater than the surface water conductivity if the volumetric water content is lower than 0.564, the maximum value applied. The measurements give evidence for the presence of a relaxation at about 3 MHz for all samples examined.  相似文献   

4.
Soil water repellency may be characterized in terms of the delayed infiltration time of a water droplet resting on the soil surface, which is, water drop penetration time (WDPT), or repellency persistence. Such repellency persistence varies nonlinearly with soil water content (θg), although no models have been proposed to reproduce the variation of WDPT with θg in soils. Dynamic factor analysis (DFA) is used to identify two common patterns of unexplained variability in a scattered dataset of WDPT versus θg measurements. A four‐parameter lognormal distribution was fitted to both common patterns obtained by DFA, and these were combined additively in a weighted multiple linear bimodal model. We show how such an empirical model is capable of reproducing a large variety of WDPT versus θg curve shapes (N = 80) both within a wide range of measured WDPTs (0–17 000 s) and for samples with organic matter content ranging from 21·7 to 80·6 g (100 g)?1. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

5.
Rock physics analysis plays a vital role in time‐lapse seismic interpretation because it provides the link between changes in rock and fluid properties and the resulting seismic data response. In this case study of the Schiehallion Field, we discuss a number of issues that commonly arise in rock physics analyses for time‐lapse studies. We show that:
  • 1 Logarithmic fits of dry bulk (Kdry) and shear (Gdry) moduli vs. effective pressure (Peff) are superior to polynomial fits.
  • 2 2D surface fits of Kdry and Gdry over porosity (φ) and effective pressure using all the core data simultaneously are more useful and accurate than separate 1D fits over φ and Peff for each individual core.
  • 3 One average set (facies) of Kdry(φ, Peff) and Gdry(φ, Peff) can be chosen to represent adequately the entire Schiehallion reservoir.
  • 4 Saturated velocities and densities modelled by fluid substitution of Kdry(φ, Peff), Gdry(φ, Peff) and the dry bulk density ρdry(φ) compare favourably with well‐log velocities and densities.
  • 5 P‐ and S‐wave impedance values resulting from fluid substitution of Kdry(φ, Peff), Gdry(φ, Peff) and ρdry(φ) show that the largest impedance changes occur for high porosities and low effective pressures.
  • 6 Uncertainties in Kdry(φ, Peff) and Gdry(φ, Peff) derived for individual cores can be used to generate error surfaces for these moduli that represent bounds for quantifying uncertainties in seismic modelling or pressure–saturation inversion.
  相似文献   

6.
Published accounts of outburst floods from glacier-dammed lakes show that a significant number of such floods are associated not with drainage through a tunnel incised into the basal ice—the process generally assumed—but rather with ice-marginal drainage, mechanical failure of part of the ice dam, or both. Non-tunnel floods are strongly correlated with formation of an ice dam by a glacier advancing from a tributary drainage into either a main river valley or a pre-existing body of water (lake or fiord). For a given lake volume, non-tunnel floods tend to have significantly higher peak discharges than tunnel-drainage floods. Statistical analysis of data for floods associated with subglacial tunnels yields the following empirical relation between lake volume V and peak discharge Qp : Qp = 46V0.66 (r2 = 0.70), when Qp is expressed in metres per second and V in millions of cubic metres. This updates the so-called Clague–Mathews relation. For non-tunnel floods, the analogous relation is Qp = 1100V0.44 (r2 = 0.58). The latter relation is close to one found by Costa (1988) for failure of constructed earthen dams. This closeness is probably not coincidental but rather reflects similarities in modes of dam failure and lake drainage. We develop a simple physical model of the breach-widening process for non-tunnel floods, assuming that (1) the rate of breach widening is controlled by melting of the ice, (2) outflow from the lake is regulated by the hydraulic condition of critical flow where water enters the breach, and (3) the effect of lake temperature may be dealt with as done by Clarke (1982). Calculations based on the model simulate quite well outbursts from Lake George, Alaska. Dimensional analysis leads to two approximations of the form QpVqf(hi, θ0), where q = 0.5 to 0.6, hi is initial lake depth, θ0 is lake temperature, and the form of f (hi, θ0) depends on the relative importance of viscous dissipation and the lake's thermal energy in determining the rate of breach opening. These expressions, along with the regression relations, should prove useful for assessing the probable magnitude of breach-type outburst floods.  相似文献   

7.
论确定地震震源深度的地震宏观方法及其解答   总被引:1,自引:0,他引:1  
作者基于在近场震源不能视作点源的认识,分析导出了相当于线、面源幅射场的深度公式,进而对源的类型不予作假定,将源指标的几何扩散率n作为待定参数,导出了一个物理意义明确且普遍化的震源深度公式。该公式的诸种解法中,以计算方法较准确,精度高;作图求解有较直观的优点,但准确性差,精度低;图算法仅作为获取粗略值和考察数据的均匀性用。用该公式对我国发生的一些不同震级的地震的震源深度进行了计算,经对比,效果很好。利用计算所得震源类型,初步研究了这些地震的震源几何学,并且首次给出了宏观地震震源深度测定值的标准误差。  相似文献   

8.
Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm3 undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (Ks), saturated moisture content (θs), residual moisture content (θr), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. Ks of both horizons, θs of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and Ks and θs of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

9.
In the theoretical part of the present paper, formulas have been analyzed for a magnetic dipole in a homogeneous and unbounded medium. The magnetic field is elliptically polarized in the region between the quasistatic zone and the far field. Since the position and the shape of the polarization ellipses depend on the complex wave number, k, it is possible to determine k by measuring the polarization ellipses. From k, the conductivity and the dielectric constant of the medium are easily calculated. The functions required for the measuring method have been computed and plotted in graphs. In the experimental part it was examined how far the theory may be applied to measurements of propagation through rock at frequencies ranging from 100-1000 kHz. These measurements showed that reasonably defined mean values of rock parameters can be given only if the deviations of the field from the theoretically expected field are not too high. These deviations have been named field distortions and have been examined by means of statistical methods (variance ratio tests). Gallery cavity and inhomogeneity or anisotropy of the medium account for these distortions.  相似文献   

10.
探讨了在电偶源频率电磁测深中利用磁场比值定义的视电阻率ρzy的物理意义及与其相位间的频散关系,给出了在频散关系的基础上电视电阻率资料估算相位的方法,并证明所估算的相位与理论模型计算值是一致的研究了ρzy与ψzy联合反演问题,并研制出具有较好应用效果的实用化软件。  相似文献   

11.
The ion microprobe at Johnson Space Center has been calibrated for in situ water determinations on a 10-μm scale over the range 0.2 wt.% H2O to 1.8, 6.8, and 3.7 wt.%, for basaltic, albitic, and rhyolitic glasses, respectively. The basalt glass calibration curve differs substantially from those of albite and rhyolite glasses, indicating a need to carefully match composition and/or melt structure between H2O standards and unknowns.A value for the diffusivity of water as a function of concentration and time has been calculated from water diffusion profiles measured in rhyolite glasses prepared at 850°C and 700 barsPt(H2O) [1]. Transient diffusion into a semi-infinite medium is described by the equation:?(φ/2)?¸/?φ=?(Dw?¸/?φ)/?φ #x003B8;=1, φ=0, θ→ 0, θ→∞, wherex =distance from the cylinder edge,t =time,C0 =initial concentration,Cs =concentration at the edge,C =concentration at x,θ = C ? C0/Cs ? C0,φ = x/t1/2, andDw =diffusivity of water. An iterative technique has been used to calculate solutions to the diffusion equation as a function ofDw [2]. Comparison of these solutions with the ion probe data indicate that, for0.2wt.% ≤ C ≤ 3.7wt.%H2O,Dw can be described by an exponential function of θ, of the formDw = D0exp(bθ), withD0 (i.e.,Dw at 0.2%) = (0.8?2.2) × 10?8 cm2/s and2 ≤ b ≤ 4.  相似文献   

12.
Summary

The urgent need for planning information on the effect of changes in land use on water resources in East Africa has necessitated the use of intensive methods of experimental catchment area research in order to produce data on the water balance of different vegetation covers in a matter of years rather than decades.

Quantitative water balance studies require an intensive network of raingauges to estimate the volumetric water input with an accuracy comparable with the measurement of outflow. Observations of the soil moisture status and energy balance, in addition to those of rainfall and streamflow are necessary to provide independent checks for “leaks” from the catchments. The successful application of these methods is illustrated from the results of three catchment area experiments in Kenya and Tanzania. The water use of each vegetational complex is characterized by the ratio of the transpiration, E t, to the evaporative demand from an open water surface, E o. This ratio is shown to vary little from year to year despite considerable variation in E t and E o.

An intensive method of analysis of stormflow response, based on the construction of a prediction equation relating stormflow to rainfall quantity and intensity and to antecedent surface soil moistrue condition, is described. Results from the application of the method in one of the catchments are presented in detail.  相似文献   

13.
Abstract

Magnetic instabilities play an important role in the understanding of the dynamics of the Earth's fluid core. In this paper we continue our study of the linear stability of an electrically conducting fluid in a rapidly rotating, rigid, electrically insulating spherical geometry in the presence of a toroidal basic state, comprising magnetic field BMB O(r, θ)1ø and flow UMU O(r, θ)1ø The magnetostrophic approximation is employed to numerically analyse the two classes of instability which are likely to be relevant to the Earth. These are the field gradient (or ideal) instability, which requires strong field gradients and strong fields, and the resistive instability, dependent on finite resistivity and the presence of a zero in the basic state B O(r,θ). Based on a local analysis and numerical results in a cylindrical geometry we have established the existence of the field gradient instability in a spherical geometry for very simple basic states in the first paper of this series. Here, we extend the calculations to more realistic basic states in order to obtain a comprehensive understanding of the field gradient mode. Having achieved this we turn our attention to the resistive instability. Its presence in a spherical model is confirmed by the numerical calculations for a variety of basic states. The purpose of these investigations is not just to find out which basic states can become unstable but also to provide a quantitative measure as to how strong the field must become before instability occurs. The strength of the magnetic field is measured by the Elsasser number; its critical value c describing the state of marginal stability. For the basic states which we have studied we find c 200–1000 for the field gradient mode, whereas for the resistive modes c 50–160. For the field gradient instability, c increases rapidly with the azimuthal wavenumber m whereas in the resistive case there is no such pronounced difference for modes corresponding to different values of m. The above values of c indicate that both types of instability, ideal and resistive, are of relevance to the parameter regime found inside the Earth. For the resistive mode, as is increased from c, we find a shortening lengthscale in the direction along the contour BO = 0. Such an effect was not observable in simpler (for example, cylindrical) models.  相似文献   

14.
Abstract

The effective porosity θ e for partially penetrated aquifers was determined. The model basin sandy aquifer available in the Centre was used. The values obtained for θ e were in good agreement with the adopted values.  相似文献   

15.
From a great variety of in situ shear wave experiments, i.e., reflection, refraction and borehole surveys in the shallow sediments of the north German plains, several specific properties have been derived. Shear waves (S) differ from compressional waves (P) in that:
  • 1 they are not affected by the degree of water saturation. Thus, they provide a better correlation between the velocity Vs and (solid) lithology;
  • 2 they generally have lower frequencies, but shorter wavelength and, hence, a better resolution of thin layers;
  • 3 they have lower absorption Qs?1 and hence a better penetration in partially saturated and gas-containing sediments than P-waves.
Correlations have been established between Vs and the confining pressure and between reduced Vs values and several lithological parameters like the grain size of sandy material. More lithological and hydrological information is obtained by using S- and P-wave surveys along the same profile. The best information on a sedimentological structure is obtained by the simultaneous observation of Vs, Vp, Qs and Qp.  相似文献   

16.
A key non-linear mechanism in a strong-field geodynamo is that a finite amplitude magnetic field drives a flow through the Lorentz force in the momentum equation and this flow feeds back on the field-generation process in the magnetic induction equation, equilibrating the field. We make use of a simpler non-linear?α?2-dynamo to investigate this mechanism in a rapidly rotating fluid spherical shell. Neglecting inertia, we use a pseudo-spectral time-stepping procedure to solve the induction equation and the momentum equation with no-slip velocity boundary conditions for a finitely conducting inner core and an insulating mantle. We present calculations for Ekman numbers (E) in the range 2.5× 10?3 to 5.0× 10?5, for?α?=α 0cos?θ?sin?π?(r?ri ) (which vanishes on both inner and outer boundaries). Solutions are steady except at lower E and higher values of?α?0. Then they are periodic with a reversing field and a characteristic rapid increase then equally rapid decrease in magnetic energy. We have investigated the mechanism for this and shown the influence of Taylor's constraint. We comment on the application of our findings to numerical hydrodynamic dynamos.  相似文献   

17.
Abstract

The problem of the removal of the degeneracy of the patterns of convective motion in a spherically symmetric fluid shell by the effects of rotation is considered. It is shown that the axisymmetric solution is preferred in sufficiently thick shells where the minimum Rayleigh number corresponds to degree l = 1 of the spherical harmonics. In all cases with l > 1 the solution described by sectional spherical harmonics Yl l (θ,φ) is preferred.  相似文献   

18.
Wide‐angle seismic data acquired by use of air‐guns and ocean bottom seismometers (OBS) contain strong direct water arrivals and multiples, generally considered as noise and thus not included in the modelling. However, a recent study showed that standard ray‐tracing modelling of the water multiples recorded off the Bear Island, North Atlantic, provided a reliable estimate of the velocity distribution in the water layer. Here, we demonstrate that including the amplitudes in the modelling provide valuable information about the VP contrast at the seafloor, as well as the VP/VS ratio and attenuation (QP) of the uppermost sediments. The VP contrast at the seafloor is estimated at about 250 m/s, within a precision of approximately ±30 m/s. The VP/VS ratio in the uppermost sedimentary layer is modelled in the range 2.25–2.50 and the QP factor is estimated at 1000 for the water, 30–50 for the uppermost layer and 40–50 for the second sedimentary layer. The values obtained for the sediments suggest a lithology dominated by silty clays, with porosity below average.  相似文献   

19.
The objective of this paper is to simulate the progress of the soil water content distribution in the soil profile with a water table at the bottom of the soil profile during ponding irrigation. This simulation can be done by solving the two‐dimensional Richards's equation for the assimilation of the advancing water jet, which uses the conditions of the two exponential functional forms k = ks eαψ and θ = θr + (θs − θr) eαψ to represent the hydraulic conductivity and volumetric water content, with ψ the pressure as the third variable. We assume that the ground surface becomes ponded and saturated as soon as the water flux passes the dry ground surface. By the technique of transformation, the analytical solution of these two‐dimensional Richards' equations has enabled figures of volumetric water content distribution to be obtained in successive time periods after irrigation. For the example of loam soil, it can simulate the variation of volumetric water content during and after irrigation in the soil profile. The analytical solutions of this paper reflect the real situation simulated, and can be applied to verify those complicated solutions from other analytical models. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

20.
Among the components of normal dipole fields useful in high-frequency electric prospecting, the electric components E and E y of the normal field of the horizontal electric dipole are considerably complicated in structure. By checking these values by means of computer mathematics in combination with the development of the classical (dynamical) theory of interference soundings, the above expressions alluded to in geophysical literature were shown to contain considerable errors. The present paper is devoted to this question.  相似文献   

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

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