于桥水库流域农业非点源磷污染控制区划研究

以于桥水库流域28个次级小流域作为区划对象,以影响该区域农业非点源磷污染的6个主要因子,包括降雨侵蚀力因子(R),土壤可蚀性因子(K),地形因子(LS),植被覆盖因子(C)以及土壤有效磷含量(SAP)和化肥磷施用量(FP),作为区划因子,采用系统聚类法进行区划,将于桥水库流域分为4类磷污染控制类型区,分析了各个类型区的生态环境特征、农业非点源磷污染特征和主要影响因子,提出相应的控制措施。为进一步进行更深入的研究提供科学依据。     

基于地震叠后反演与属性分析的页岩气有利储层预测

中国南方页岩气勘查区块多为勘探新区,仅部署少量二维地震和地质钻井,且资料品质普遍较差,页岩气储层预测很难像成熟勘探区一样开展系统性研究工作,如何对现有数据进行充分挖掘是解决勘探新区页岩气储层预测的关键。为此,综合利用不依赖初始模型的叠后约束稀疏脉冲反演方法和多种频率类属性分析技术,开展页岩气有利储层预测。地震反演能够获得多种参数参与储层预测,属性分析能够从多种角度识别储层,通过多种信息综合应用和叠合分析以提高储层预测精确性和有效性,为新区页岩气勘探提供了一种新的思路。该项技术应用于南方某区块,圈定了页岩气的有利目标区,在该区部署的两口水平井均获得高产工业气流。     

Seismic Ground Motion Expected for the Eastern District of Naples

The seismic ground motion of a test area in the eastern district of Naples is computed with a hybrid technique based on the mode summation and the finite difference methods. This technique allows us the realistic modelling of source and propagation effects, including local soil conditions. In the modelling, we consider the 1980 Irpinia earthquake, a good example of strong shaking for the area of Naples, which is located about 90 km from the epicenter.The detailed geological setting is reconstructed from a large number of drillings. The sub-soil is mainly formed by alluvial (ash, stratified sand and peat) and pyroclastic materials overlying a pyroclastic rock (yellow neapolitan tuff), representing the neapolitan bedrock. The detailed information available on mechanical properties of the sub-soil and its geometry warrants the application of the sophisticated hybrid technique.As expected, the sedimentary cover causes an increase of the signal's amplitudes and duration. If thin peat layers are present, the amplification effects are reduced, and the peak ground accelerations are similar to those observed for the bedrock model. This can be explained by the backscattering of wave energy at such layers, that tend to seismically decouple the upper from the lower part of the structure.For SH-waves, the influence of the variations of the S-wave velocities on the spectral amplification is studied, by considering locally measured velocities and values determined from near-by down-hole measurements. The comparison between the computed spectral amplifications confirms the key role of an accurate determination of the seismic velocities of the different layers.The comparison performed between a realistic 2-D seismic response and a standard 1-D response, based on the vertical propagation of waves in a plane layered structure, shows considerable difference, from which it is evident that serious caution must be taken in the modelling of expected ground motion at a specific site.     

Crustal structure across the Grand Banks–Newfoundland Basin Continental Margin – II. Results from a seismic reflection profile

New multichannel seismic reflection data were collected over a 565 km transect covering the non-volcanic rifted margin of the central eastern Grand Banks and the Newfoundland Basin in the northwestern Atlantic. Three major crustal zones are interpreted from west to east over the seaward 350 km of the profile: (1) continental crust; (2) transitional basement and (3) oceanic crust. Continental crust thins over a wide zone (∼160 km) by forming a large rift basin (Carson Basin) and seaward fault block, together with a series of smaller fault blocks eastwards beneath the Salar and Newfoundland basins. Analysis of selected previous reflection profiles (Lithoprobe 85-4, 85-2 and Conrad NB-1) indicates that prominent landward-dipping reflections observed under the continental slope are a regional phenomenon. They define the landward edge of a deep serpentinized mantle layer, which underlies both extended continental crust and transitional basement. The 80-km-wide transitional basement is defined landwards by a basement high that may consist of serpentinized peridotite and seawards by a pair of basement highs of unknown crustal origin. Flat and unreflective transitional basement most likely is exhumed, serpentinized mantle, although our results do not exclude the possibility of anomalously thinned oceanic crust. A Moho reflection below interpreted oceanic crust is first observed landwards of magnetic anomaly M4, 230 km from the shelf break. Extrapolation of ages from chron M0 to the edge of interpreted oceanic crust suggests that the onset of seafloor spreading was ∼138 Ma (Valanginian) in the south (southern Newfoundland Basin) to ∼125 Ma (Barremian–Aptian boundary) in the north (Flemish Cap), comparable to those proposed for the conjugate margins.     

Small-scale heterogeneities below the Gräfenberg array, Germany from seismic wavefield fluctuations of Hindu Kush events

Small-scale elastic heterogeneities (<5  km) are found in the upper lithosphere underneath the Gräfenberg array, southeast Germany. The results are based on the analysis of broadband recordings of 17 intermediate-depth (201–272  km) events from the Hindu Kush region. The wavefront of the first P arrival and the following 40  s coda are separated into coherent and incoherent (scattered) parts in the frequency range from 0.05 to 5  Hz. The frequency-dependent intensities of the mean and fluctuation wavefields are used to describe the scattering characteristics of the lithosphere underneath the receivers. It is possible to discriminate a weak-fluctuation regime of the wavefield in the frequency range below approximately 1.5–2.5  Hz and a strong-fluctuation regime starting at 2.0–2.5  Hz and continuing to higher frequencies. In order to explain the observed wavefield fluctuations, an approach with seismic scattering at random media-type structures is proposed. The preferred model contains heterogeneities with 3–7 per cent perturbations in seismic velocity and correlation lengths of 0.6–4.8  km in the crust. This is compatible with models from active seismic experiments. Scattering in the lithospheric mantle is not required, but cannot be excluded at weak velocity contrasts (<3 per cent).     

Soil bentonite wall protects foundation from thrust faulting: analyses and experiment

When seismic thrust faults emerge on the ground surface,they are particularly damaging to buildings,bridges and lifelines that lie on the rupture path.To protect a structure founded on a rigid raft,a thick diaphragm-type soil bentonite wall(SBW) is installed in front of and near the foundation,at sufficient depth to intercept the propagating fault rupture.Extensive numerical analyses,verified against reduced–scale(1 g) split box physical model tests,reveal that such a wall,thanks to its high deformability and low shear resistance,"absorbs" the compressive thrust of the fault and forces the rupture to deviate upwards along its length.As a consequence,the foundation is left essentially intact.The effectiveness of SBW is demonstrated to depend on the exact location of the emerging fault and the magnitude of the fault offset.When the latter is large,the unprotected foundation experiences intolerable rigid-body rotation even if the foundation structural distress is not substantial.     

Modelling and inversion of single-ended refraction data from the shot gathers of multifold deep seismic reflection profiling—an approach for deriving the shallow velocity structure

A multifold crustal-scale deep seismic near-vertical reflection profile generates a large number of single-ended shot gathers, which provide redundant data sets because of overlapping coverage of the shallow refractors. We present an approach for deriving the shallow velocity structure by modelling and inversion of single-ended seismic refraction first arrival traveltime data. We apply this method to a data set acquired with a 12-km long spread with 100 m spacing of shots and receivers, of the Neoproterozoic Marwar basin in the NW Indian shield. The approach is shown to be quite successful for delineating the shallow refractor depths, steep dips and velocities, even in the absence of regular reverse refraction profiles. The study reveals two-layered sedimentary formations, Malani volcanics and a complicated basement configuration of the Marwar basin, and provides a measure of resolution and uncertainty of the estimated model parameters. A seismic section of the near-trace gather is found to be qualitatively consistent with the derived structural features of the basin. The relative highs and lows, observed in the Bouguer gravity profile, further corroborate the derived velocity model. The present approach can be especially useful in offshore areas and elsewhere, where the single-ended multifold seismic profiles are the only available data sets.     

Generalized Born inversion of seismic reflection data

Summary . Born inverse methods give accurate and stable results when the source wavelet is impulsive. However, in many practical applications (reflection seismology) an impulsive source cannot be realized and the inversion needs to be generalized to include an arbitrary source function. In this paper, we present a Born solution to the seismic inverse problem which can accommodate an arbitrary source function and give accurate and stable results. It is shown that the form of the generalized inversion algorithm reduces to a Wiener shaping ***filter, which is solved efficiently using a Levinson recursion algorithm. Numerical examples of synthetic and real field data illustrate the validity of our method.