隐伏层的形成与识别

利用折射波场的特点，研究隐伏层，建立不同性质中间层的临界方程和隐伏层形成条件的公式．研究表明，隐伏层形成与剖面中每个地层的厚度、波速、以及折射面的倾角有关．厚度薄而波速低的中间层容易形成隐伏层．还用τ－ｐ变换方法来识别和提取隐伏层的视速度与截距时间，并给出了模型的结果．     

Identification of the best architecture of a multilayer perceptron in modeling daily total ozone concentration over Kolkata,India

Autoregressive neural network (AR-NN) models of various orders have been generated in this work for the daily total ozone (TO) time series over Kolkata (22.56°N, 88.5°E). Artificial neural network in the form of multilayer perceptron (MLP) is implemented in order to generate the AR-NN models of orders varying from 1 to 13. An extensive variable selection method through multiple linear regression (MLR) is implemented while developing the AR-NNs. The MLPs are characterized by sigmoid non-linearity. The optimum size of the hidden layer is identified in each model and prediction are produced by validating it over the test cases using the coefficient of determination (R ²) and Willmott’s index (WI). It is observed that AR-NN model of order 7 having 6 nodes in the hidden layer has maximum prediction capacity. It is further observed that any increase in the orders of AR-NN leads to less accurate prediction.     

SOME ASPECTS OF HANDLING VELOCITY INVERSION AND HIDDEN LAYER PROBLEMS IN SEISMIC REFRACTION WORK*

There are two types of masked layers in seismic refraction work: the velocity reversal (low-velocity layer) and the hidden layer (insufficient velocity contrast or layer thickness). On the basis of an analytical formulation of the general case of a masked layer under an overburden of plane and parallel multiple refractors the two limiting cases are discussed: the solution resulting from an uncritical interpretation of the measured time-distance curve and the blind zone solution. Between these two limiting cases there is a variety of possible masked layer solutions. These no-blind zone solutions—as well as the blind zone solution itself—are formulated separately for the velocity inversion and the hidden layer case. For the evaluation of some no-blind zone solution a diagram is presented which can be used for any case of multiple refractors in the overburden of the masked layer. However, it is only for the three- and the four-layer case that a blind zone interpretation by use of diagrams is advisable. Such diagrams are presented together with the basic sets of formulae which contain as parameters only ratios of velocities and layer thicknesses. As the velocity of the masked layer is usually unknown the diagrams are principally constructed to show the dependence on the masked layer velocity. This is useful for estimation of the largest possible error.     

Magnitude scale structure from filtered broadbandP-wave records

A series of Hokkaido events, recorded by the FBV Broadband Seismograph System at the KHC Seismic Station, is used to study the structure of the earthquake magnitude scale on the basis of maximum velocity amplitudesA _vmax of teleseismicPwaves in different period bands. Amplitude-periodband (APB) diagrams are constructed for each event. According to the shape of the APB diagrams the events investigated can be divided into three types: (a) events with largestA _vmax values in the intermediate period range (periods ofA _vmax from 2.2 to 23 sec), (b) events with largestA _vmax values in the short-period range (periods ofA _vmax from 1 to 2 sec), (c) events exhibiting anomalous APB diagrams. Type (a) events seem to represent the process of wave generation that prevails for shallow earthquakes. Type (b) events approach to explosive-like generation of seismic waves. The nature of the exceptionally occurring type (c) events must be clarified in further investigations. The influence of the type of earthquake on the magnitude values estimated on the basis of standard class A and B (short-period and intermediate-band) seismograms is demonstrated. It appears that for estimating correct values of earthquake magnitudes complementary information on the process of seismic wave generation in the focus is necessary. At teleseismic distances this information can be obtained from either APB diagrams or amplitude spectra ofP waves recorded, e.g., by broadband velocity sensing instruments.     

REFRACTION LE LONG DES BANCS MINCES RAPIDES ET EFFET D'ECRAN POUR LES MARQUEURS PROFONDS *

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 k₁ 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.     

A note on the approximate determination of the dispersion of Love waves in a laterally nonhomogeneous layer lying over a homogeneous half-space

Rayleigh's principle and the concept of the local wave number have been utilised for the approximate determination of the dispersion of Love waves propagating in a laterally heterogeneous layer lying over a homogeneous half-space. The shear wave velocity and the rigidity in the surface layer have been assumed to decrease with the increase of the lateral distance from the origin. The range of validity of the dispersion equation obtained by this method has been examined critically. It was found that: (a) for existence of Love waves the minimum value of shear wave velocity in the layer must be less than that in the matter below, and (b) the phase velocity of Love waves decreases with the increase of the lateral distance from the origin.     

THE THIN-LAYER PROBLEM IN GEOELECTRICAL PROSPECTING *

For a thin highly-conducting layer with given longitudinal conductance the recurrence formulae for an n-fold horizontally stratified subsoil are established for d.c. resistivity and magnetotelluric soundings. Similarly, a thin low conductivity layer with given transverse resistance is treated in the d.c. case and a non-conducting intermediate bed in magnetotellurics. Model curves for a thin high- or low-conductivity intermediate layer in the three-layer case have been carried out, which may serve as an extension of the well-known three-layer diagrams for a Schlumberger configuration. The corresponding model curves in magnetotellurics are given. By numerical comparison of these curves with real three-layer curves some diagrams have been developed to show the allowed thicknesses of the intermediate layer in the Schlumberger case and in the case of magnetotelluric sounding.     

On results of the surface wave analyses in poroelastic media by means of the Simple Mixture Model and the Biot model

This note is devoted to an overview of the results on the number, velocity and attenuation of surface waves in two-component saturated poroelastic media. The author's numerical results which have been obtained using the Simple Mixture Model are compared to some results of other authors who carried out either an analysis by means of the classical Biot model for poroelastic media or experiments. It is shown that the Simple Mixture Model is a simplification of the Biot model. Due to its simpler form the Simple Mixture Model is more suitable for a complex analysis which depends on several parameters. Thus, the Simple Mixture Model is the only model for which, simultaneously, either the dependence on the frequency, ω, and the dependence on the bulk permeability, π (boundary porous medium/vacuum) or the dependence on the frequency, ω, and the dependence on the surface permeability, α (boundary porous medium/fluid) has been considered.     

TIME-DOMAIN ELECTROMAGNETIC SOUNDING—COMPUTATION OF MULTI-LAYER RESPONSE AND THE PROBLEM OF EQUIVALENCE IN INTERPRETATION*

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 h²ρ=constant for K-type earth sections, where h and ρ are respectively the thickness and resistivity of the intermediate layer.     

Rayleigh波勘探中“之”字形频散曲线“起跳点”频率研究

当层状介质中存在低速层的时候,实际提取到的Rayleigh波频散曲线往往会发生“之”字形回折.已有研究表明,“之”字形回折与各模式的激发能量有关.特别的,“之”字形回折的“起跳点”(发生“之”字形回折的点)与介质参数有一定关系,因此在应用中它反映了介质的一些特征.但是,这一关系是怎样的还没有人进行详细研究.本文计算了三...     

随县-西安剖面地壳结构的初步研究

本文对随县-西安测线的地震测深资料进行了初步分析,讨论了该剖面地壳内主要界面的震相特征。通过对观测走时的分层反演及综合反演,得到三层的地壳模型。地壳的平均速度为6.33km/s,厚度约34km,上地幔顶部的速度为8.10km/s。下地壳是由高速及低速层交替构成,高速层的速度为7.39km/s。构制了主要界面的剖面起伏图,并进行了讨论。     

Effect of non‐uniformity of flow on velocity and turbulence intensities over a cobble‐bed

Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?u_fw_f have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m³/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd.     

The amplitude ratioPcP/P and the core-mantle boundary

Summary Using the Haskell matrix formulation, theoretical reflection coefficient curves have been calculated for a multi-layered core-mantle boundary for comparison with observational data. Two cases are considered, first when the shear velocity in the core is equal to zero and second when the core has a finite rigidity. If the velocity contrast is large between the imbedded layer and the mantle, the reflection coefficient curves for the multi-layered medium are irregular in shape as compared to those for two half-spaces, representing the core and the mantle, respectively. The reflection coefficient curves show an oscillatory character if the imbedded layer is thick and has a high velocity contrast.The observational data consist of short-period vertical-component seismograph records ofP andPcP from nuclear explosions in the Aleutian chain, Nevada, Novaya Zemlya, Kazakh and Sahara. Attenuation and geometrical spreading are taken into consideration. Four different models for the quality factorQ are applied to the observational data. The data are found to be much affected by theQ-model used for the corrections.Based on proposedQ-values, a model for the core-mantle boundary is found, characterized by two low-velocity layers at the bottom of the mantle. The thicknesses are 16.10 km (outer layer) and 19.96 km (inner layer), the compressional wave velocities 12.17 km/sec and 10.94 km/sec and the shear wave velocities are 6.29 km/sec and 5.33 km/sec, respectively. A better fit to this model is found when in addition the shear velocity in the outer core is 2.20 km/sec and the density ratio at the core-mantle boundary is 1.07. In other words, the observations favour a layer of finite rigidity in the outer core rather than a fluid one.     

COMPUTING FIELD STATICS WITH THE HELP OF SEISMIC TOMOGRAPHY*

Field static corrections in general need be applied to all onshore seismic reflection data to eliminate the disturbing effects a weathering layer or near-surface low velocity zone has on the continuity of deep seismic reflections. The traveltimes of waves refracted at the bottom of the low velocity zone (or intermediate refracting interfaces) can often be observed as first breaks on shot records and used to develop a laterally inhomogeneous velocity model for this layer, from which the field static corrections can then be obtained. A simple method is described for computing accurate field statics from first breaks. It is based on a linearization principal for traveltimes and leads to the algorithms that are widely and successfully applied within the framework of seismic tomography. We refine an initial model for the low velocity layer (estimated by a standard traveltime inversion technique) by minimizing the errors between the observed first arrivals on field records and those computed by ray theory through an initial model of the low velocity layer. Thus, one can include more lateral velocity variations within the low velocity layers, which are important to obtain good field static corrections. Traditional first break traveltime inversion methods cannot, in general, provide such refined velocity values. The technique is successfully applied to seismic data from the Amazon Basin. It is based on a simple model for the low velocity layer that consists of an undulating earth surface and one planar horizontal refractor overlain by a laterally changing velocity field.     

Wave propagation in two layered medium under initial stresses

Summary In this paper, the frequency equation for phase velocity of waves propagated in a laminated medium consisting of two eleastic layers of finite thickness under initial stresses, has been obtained. It has been shown that when wave length becomes very small compared to the thickness of each layer, the wave approaches two Rayleigh waves at the two outer surfaces with the possibility of Stoneley waves at the interface. The propagation ofSH-waves in the composite medium under initial stresses has also been discussed. A particular case has been taken to find the velocity of Love wave in the homogeneous half space under initial compressive stresses.Biot's incremental deformation theory has been used.     

Velocity Structure of the West-Bengal Sedimentary Basin, India along the Palashi-Kandi Profile Using a Travel-time Inversion of Wide-angle Seismic data and Gravity Modeling – An Update

The 2-D shallow velocity structure along the north-south Palashi-Kandi profile in the West Bengal sedimentary basin has been updated by travel-time inversion of seismic refraction, wide-angle reflection and gravity data. A six-layer shallow model up to a depth of about 7 km has been derived. The first layer, which has an average velocity of 2.0 kms^?1, represents the alluvium deposit, which rests over the shale formation with average velocity of 3.0 kms^?1. The thin (200 m) Sylhet limestone, observed at a nearby Palashi well, remains hidden in the present data set. Hence a 200-m thin layer with a velocity of 3.7 kms^?1, corresponding to the Sylhet limestone, has been assumed to be present throughout the profile. The fourth layer with a velocity of 4.5–4.7 kms^?1 at a depth of 1.7–2.4 km represents the Rajmahal traps. The ‘skip’ phenomenon and rapid amplitude decay of first arrivals indicate a low-velocity zone (LVZ) in the study area. Using the ‘skip’ phenomena and wide-angle reflection data, identified on seismograms, the LVZ with a velocity of 4.0 kms^?1, indicating the Gondwana sediments, has been delineated below the Rajmahal traps. The next layer with a velocity 5.4–5.6 kms^?1 overlying the crystalline basement (5.8–6.25 kms^?1) may be associated with the Singhbhum group of meta volcanic rock that has been exposed in the western part of the basin. The basement lies at a variable depth of 4.9 to 6.8 km. The overall uncertainties of various velocity and boundary nodes are ± 0.15 kms^?1 and ± 0.5 km, respectively. The elevated basement feature in the north might have acted as a structural barrier for the deposition of Sylhet limestone during the Eocene epoch. The seismically derived shallow structure correctly explains the observed Bouguer gravity anomaly along the profile. The addition of reflections in the present analysis provides a stronger control on the depths and velocities of basement and overlying sedimentary formations, compared to the earlier model derived mainly by the first arrival seismic data.     

Hydrological characteristics of vegetated river flows: a laboratory flume study

Abstract

Laboratory flume experiments were undertaken to measure the vertical profiles of mean flow velocity for three different flow discharges and four different stem densities of Hydrilla verticillata. The data were used to calculate three parameters, namely Manning's roughness coefficient, the Reynolds number and the Froude number. In addition, empirical equations were obtained for the vertical distribution of measured flow velocity within the transitional zone and above the plant canopy. The results show that: (a) the vertical distribution of measured flow velocity exhibits three zone profiles; (b) Manning's roughness coefficient decreases with increasing depth-averaged flow velocity; (c) the relationship between Manning's roughness coefficient and the depth-averaged flow velocity is within the smooth left inverse curve; (d) Manning's roughness coefficient significantly changes with increasing density of Hydrilla; (e) the Froude number is independent of the density of Hydrilla; and (f) both the Reynolds number and the Froude number increase with increasing depth-averaged flow velocity.

Citation Shi, J.Z., Li, Y.-H., Hughes, J.M.R., and Zhao, M., 2013. Hydrological characteristics of vegetated river flows: a laboratory flume study. Hydrological Sciences Journal, 58 (5), 1047–1058.

Editor Z.W. Kundzewicz     

Near-bed mass flux profiles in aeolian sand transport: high-resolution measurements in a wind tunnel

Vertical profiles of the streamwise mass flux of blown sand in the near-bed (< 17 mm) region are analysed from high-resolution measurements made using an optical sensor in a wind tunnel. This analysis is complemented by detailed measurements of mass flux and mean velocity profiles throughout the boundary layer depth (0·17 m) using passive, chambered sand traps of small dimensions and armoured thermal anemometers, respectively. The data permit a preliminary analysis of the relations between the observed forms of the profiles of near-bed fluid stress and horizontal mass flux within a carefully conditioned boundary layer. Profiles of mass flux density are found to be characterized by three regions of differing gradient with transitions at about 2 mm and 19 mm above the bed. The exponential decay of mass flux with height is confirmed for elevations above 19 mm, and when plotted as a function of u_*²/g (a parameter of mean vertical trajectory height in saltation), the gradient of mass flux in this region scales with the wake-corrected friction velocity (u), where u > 0·30 m s⁻¹. A separate near-bed region of more intense transport below 19 mm is identified which carries 80 per cent of the total mass flux. This region is evident in some previous field and wind tunnel data but not in profiles simulated by numerical models. Ventilated passive sand traps underestimate mass flux in this region by 37 per cent. At slow or moderate wind speeds a third significant region below 2 mm is observed. These regions are likely to be related to grain populations in successive saltation, low-energy ejections and intermittent bed contact, respectively. Optical measurements reveal locally high grain concentrations at some elevations below 5 mm; these heights scale with transport rate, mass flux gradient and wind speed. Copyright © 1999 John Wiley & Sons, Ltd.     

A note on the heat transfer by the axisymmetric thermal convection in a rotating fluid annulus

Abstract

Some new measurements are presented of the axisymmetric heat transport in a differentially heated rotating fluid annulus. Both rigid and free upper surface cases are studied, for Prandtl numbers of 7 and 45, from low to high rotation rates. The rigid lid case is extended to high rotation rates by suppressing the baroclinic waves, that would normally develop at some intermediate rotation rate, with the use of sloping endwalls.

A parameter P is defined as the square of the ratio of the (non-rotating) thermal sidewall layer thickness to the Ekman layer thickness. For small P the heat transport remains unaffected by the rotation, but as P increases to order unity the Ekman layer becomes thin enough to inhibit the radial mass transport, and hence the heat flux. No explicit Prandtl number dependence is observed. Also this scaling allows the identification of the region in which the azimuthal velocity reaches its maximum. Direct comparisons are drawn with previous experimental and numerical results, which show what can be interpreted as an inhibiting effect of increasing curvature on the heat transport.     

Artificial Neural Network Prediction of Tropospheric Ozone Concentrations in Istanbul,Turkey

Tropospheric (ground‐level) ozone has adverse effects on human health and environment. In this study, next day's maximum 1‐h average ozone concentrations in Istanbul were predicted using multi‐layer perceptron (MLP) type artificial neural networks (ANNs). Nine meteorological parameters and nine air pollutant concentrations were utilized as inputs. The total 578 datasets were divided into three groups: training, cross‐validation, and testing. When all the 18 inputs were used, the best performance was obtained with a network containing one hidden layer with 24 neurons. The transfer function was hyperbolic tangent. The correlation coefficient (R), mean absolute error (MAE), root mean squared error (RMSE), and index of agreement or Willmott's Index (d₂) for the testing data were 0.90, 8.78 µg/m³, 11.15 µg/m³, and 0.95, respectively. Sensitivity analysis has indicated that the persistence information (current day's maximum and average ozone concentrations), NO concentration, average temperature, PM₁₀, maximum temperature, sunshine time, wind direction, and solar radiation were the most important input parameters. The values of R, MAE, RMSE, and d₂ did not change considerably for the MLP model using only these nine inputs. The performances of the MLP models were compared with those of regression models (i.e., multiple linear regression and multiple non‐linear regression). It has been found that there was no significant difference between the ANN and regression modeling techniques for the forecasting of ozone concentrations in Istanbul.