A method for spatial heterogeneity evaluation on landscape pattern of farmland shelterbelt networks: A case study in midwest of Jilin Province,China

On the basis of landscape ecology, combining the Spot 5 high resolution satellite imagery with GIS, a method evaluating the spatial heterogeneity of shelterbelts distribution at landscape scale is put forward in this paper. The distance coefficients of reasonable and existing landscape indexes of farmland shelterbelt networks were computed, and then through the classification of the distance coefficients, and the establishment of evaluation rules, the spatial heterogeneity of farmland shelterbelts was evaluated. The method can improve the evaluating system of previous studies on shelterbelts distribution, resolve the disadvantages of lacking spatiality of overall evaluation, and make the evaluation results have more directive significance for shelterbelt management. Based on this method, spatial heterogeneity of shelterbelt networks was evaluated in the midwest of Jilin Province, China. The results show that the regions with fewer shelterbelts and no closed network account for 34.7% of the total area, but only 4.9% of the area has relative reasonable pattern of shelterbelt networks. Many problems exist in the distribution pattern of shelterbelts, therefore, much attention should be paid to construct farmland shelterbelts in the study area.     

卫城油田卫22块储层非均质性评价

通过对卫城油田卫22块铸体薄片、扫描电镜以及孔隙度渗透率资料的分析,对研究区储层微观非均质性,储层层内、层间、平面非均质性做了详细的研究。认为研究区储层砂岩以细砂岩粉砂岩为主,孔喉结构为中孔细喉;储层非均质性主要受控于沉积作用和成岩作用,层内渗透率大部分为反韵律序列,各小层内部渗透率的差异性很大,表现为严重非均质型;纵向上各小层砂岩密度相差很大,但各井垂向上渗透率差异不大,以非均质型居多;平面上储层非均质性主要受沉积相带平面分布的影响,水下分流河道和河口砂坝沉积非均质性较弱,前缘沙席表现为较强的平面非均质性。     

Stochastic computational modelling of highly heterogeneous poroelastic media with long‐range correlations

The compaction of highly heterogeneous poroelastic reservoirs with the geology characterized by long‐range correlations displaying fractal character is investigated within the framework of the stochastic computational modelling. The influence of reservoir heterogeneity upon the magnitude of the stresses induced in the porous matrix during fluid withdrawal and rock consolidation is analysed by performing ensemble averages over realizations of a log‐normally distributed stationary random hydraulic conductivity field. Considering the statistical distribution of this parameter characterized by a coefficient of variation governing the magnitude of heterogeneity and a correlation function which decays with a power‐law scaling behaviour we show that the combination of these two effects result in an increase in the magnitude of effective stresses of the rock during reservoir depletion. Further, within the framework of a perturbation analysis we show that the randomness in the hydraulic conductivity gives rise to non‐linear corrections in the upscaled poroelastic equations. These corrections are illustrated by a self‐consistent recursive hierarchy of solutions of the stochastic poroelastic equations parametrized by a scale parameter representing the fluctuating log‐conductivity standard deviation. A classical example of land subsidence caused by fluid extraction of a weak reservoir is numerically simulated by performing Monte Carlo simulations in conjunction with finite elements discretizations of the poroelastic equations associated with an ensemble of geologies. Numerical results illustrate the effects of the spatial variability and fractal character of the permeability distribution upon the evolution of the Mohr–Coulomb function of the rock. Copyright © 2004 John Wiley & Sons, Ltd.     

An adaptive grid technique for minimizing heterogeneity of cells or elements

Geostatistical techniques allow simulation of properties such as porosity or conductivity on a fine scale. Typically, porous media flow modeling is performed at a coarser scale. Upscaling properties from the fine scale to the coarser scale introduces potential errors which are constrained by the degree of homogeneity of the cell or element. Adaptive grid techniques can be used to minimize the heterogeneity in the individual cells or elements, thus minimizing potential upscaling errors. A geostatistical adaptive grid (GAG) algorithm based on local minimization of heterogeneity is introduced. Local minimization allows greater control over the type of distortion permitted. Comparisons are made with a general elastic grid adjustment (GEGA) algorithm based on global minimization of heterogeneity. Several sample problems are used to test and demonstrate the two approaches.     

Anisotropy and depth‐related heterogeneity of hydraulic conductivity in a bog peat. I: laboratory measurements

Anisotropy and heterogeneity of hydraulic conductivity (K) are suspected of greatly affecting rates and patterns of ground‐water seepage in peats. A new laboratory method, termed here the modified cube method, was used to measure horizontal and vertical hydraulic conductivity (K_h and K_v) of 400 samples of bog peat. The new method avoids many of the problems associated with existing field and laboratory methods, and is shown to give relatively precise measurements of K. In the majority of samples tested, K_h was much greater than K_v, indicating that the bog peat was strongly anisotropic. Log₁₀K_h, log₁₀K_v, and log₁₀ (K_h/K_v) were found to vary significantly with depth, although none of the relationships was simple. We comment on the scale dependency of our measurements. Copyright © 2002 John Wiley & Sons, Ltd.     

Constituents of Vertical-component Coda Waves at Long Periods

Aki (1969) first modeled coda waves of a local earthquake as a superposition of scattered surface waves. This paper attempts to clarify the constituents of surface-wave coda at long periods at very long lapse times. For a large earthquake of magnitude 7 or larger, vertical component oscillation in periods from 90 to 180 s persists for more than 20 hours from the earthquake origin time. Although the early portion of the coda envelope is successfully modeled by assuming incoherent scattered Rayleigh waves by heterogeneities distributed all over the Earth, the later potion of the observed coda envelope (roughly later than 35,000 s) has systematically larger amplitude than theoretical prediction. To clarify the cause of this discrepancy, we studied the constituents of vertical-component seismograms of three large earthquakes recorded by the F-net in Japan using the f-k power spectral analysis. We found that the direct and scattered fundamental-mode Rayleigh waves of velocity about 3.7 km/s are dominant in the earlier part of each envelope. It justifies the use of a scattering model of the fundamental Rayleigh waves for synthesizing the envelope. At lapse times later than 20,000 s–35,000 s, higher modes with phase velocities around 20 km/s become dominant. The transition time to the dominance of higher modes is found to become earlier for a deeper focus earthquake. The small coda attenuation factor from (1.90±0.23) × 10⁻³ to (2.38±0.32) × 10⁻³ estimated from later coda envelopes recorded at IRIS stations distributed worldwide also agrees with the attenuation factor of spheroidal modes according to PREM. We may interpret that higher mode waves are uniformly distributed at large lapse time due to large velocity dispersion and/or scattering and they dominate over the fundamental mode waves because of smaller attenuation in the lower mantle. The coda attenuation measurement proposed by Aki is found to be useful even for long periods and at very large lapse times.     

生态空间理论及其在生物入侵研究中的应用

生态空间理论（theoryofspat ialecology）是对生态系统空间关系进行研究的一种理论，主要包括尺度、空间格局和镶嵌动态等。模型是研究空间关系的重要手段，其中莱文斯（Levins）模型、细胞自动机（cel lularautomata）模型和反应-扩散（React ion Di f fusion）模型是 3个基本的模型。尺度是生态学研究中的核心问题，对于不同的研究对象，应在相应的尺度上进行思考。空间异质性是不同尺度上的生态学过程的结果，同时它也影响各种生态学过程。随着人们对生物入侵的日益关注，空间生态学理论尤其是一些空间模型越来越多的应用于入侵和竞争的研究中。     

Feasibility of waveform inversion of Rayleigh waves for shallow shear-wave velocity using a genetic algorithm

Conventional surface wave inversion for shallow shear (S)-wave velocity relies on the generation of dispersion curves of Rayleigh waves. This constrains the method to only laterally homogeneous (or very smooth laterally heterogeneous) earth models. Waveform inversion directly fits waveforms on seismograms, hence, does not have such a limitation. Waveforms of Rayleigh waves are highly related to S-wave velocities. By inverting the waveforms of Rayleigh waves on a near-surface seismogram, shallow S-wave velocities can be estimated for earth models with strong lateral heterogeneity. We employ genetic algorithm (GA) to perform waveform inversion of Rayleigh waves for S-wave velocities. The forward problem is solved by finite-difference modeling in the time domain. The model space is updated by generating offspring models using GA. Final solutions can be found through an iterative waveform-fitting scheme. Inversions based on synthetic records show that the S-wave velocities can be recovered successfully with errors no more than 10% for several typical near-surface earth models. For layered earth models, the proposed method can generate one-dimensional S-wave velocity profiles without the knowledge of initial models. For earth models containing lateral heterogeneity in which case conventional dispersion-curve-based inversion methods are challenging, it is feasible to produce high-resolution S-wave velocity sections by GA waveform inversion with appropriate priori information. The synthetic tests indicate that the GA waveform inversion of Rayleigh waves has the great potential for shallow S-wave velocity imaging with the existence of strong lateral heterogeneity.