共查询到20条相似文献，搜索用时 109 毫秒
1.
This article describes a proposed work-sequence to generate accurate reservoir-architecture models, describing the geometry of bounding surfaces (i.e., fault locations and extents), of a structurally complex geologic setting in the Jeffara Basin (South East Tunisia) by means of geostatistical modeling. This uses the variogram as the main tool to measure the spatial variability of the studied geologic medium before making any estimation or simulation. However, it is not always easy to fit complex experimental variograms to theoretical models. Thus, our primary purpose was to establish a relationship between the geology and the components of the variograms to fit a mathematically consistent and geologically interpretable variogram model for improved predictions of surface geometries. We used a three-step approach based on available well data and seismic information. First, we determined the structural framework: a seismo-tectonic data analysis was carried out, and we showed that the study area is cut mainly by NW–SE-trending normal faults, which were classified according to geometric criteria (strike, throw magnitude, dip, and dip direction). We showed that these normal faults are at the origin of a large-scale trend structure (surfaces tilted toward the north-east). At a smaller scale, the normal faults create a distinct compartmentalization of the reservoirs. Then, a model of the reservoir system architecture was built by geostatistical methods. An efficient methodology was developed, to estimate the bounding faulted surfaces of the reservoir units. Emphasis was placed on (i) elaborating a methodology for variogram interpretation and modeling, whereby the importance of each variogram component is assessed in terms of probably geologic factor controlling the behavior of each structure; (ii) integrating the relevant fault characteristics, which were deduced from the previous fault classification analysis, as constraints in the kriging estimation of bounding surfaces to best reflect the geologic structure of the study area. Finally, the estimated bounding surfaces together with seismic data and variogram interpretations were used to obtain further insights into the tectonic evolution of the study area that has induced the current reservoirs configuration. 相似文献
2.
Annabel Kelly Richard W. England Peter K. H. Maguire 《Geophysical Journal International》2007,171(3):1172-1184
A regional model of the 3-D variation in seismic P -wave velocity structure in the crust of NW Europe has been compiled from wide-angle reflection/refraction profiles. Along each 2-D profile a velocity–depth function has been digitised at 5 km intervals. These 1-D velocity functions were mapped into three dimensions using ordinary kriging with weights determined to minimise the difference between digitised and interpolated values. An analysis of variograms of the digitised data suggested a radial isotropic weighting scheme was most appropriate. Horizontal dimensions of the model cells are optimised at 40 × 40 km and the vertical dimension at 1 km. The resulting model provides a higher resolution image of the 3-D variation in seismic velocity structure of the UK, Ireland and surrounding areas than existing models. The construction of the model through kriging allows the uncertainty in the velocity structure to be assessed. This uncertainty indicates the high density of data required to confidently interpolate the crustal velocity structure, and shows that for this region the velocity is poorly constrained for large areas away from the input data. 相似文献
3.
Shear and compressional velocity models of the mantle from cluster analysis of long-period waveforms
We present a new technique for the efficient measurement of the traveltimes of long period body wave phases. The technique is based on the fact that all arrivals of a particular seismic phase are remarkably similar in shape for a single event. This allows the application of cross-correlation techniques that are usually used in a regional context to measure precise global differential times. The analysis is enhanced by the inclusion of a clustering algorithm that automatically clusters waveforms by their degree of similarity. This allows the algorithm to discriminate against unusual or distorted waveforms and makes for an extremely efficient measurement technique.
This technique can be applied to any seismic phase that is observed over a reasonably large distance range. Here, we present the results of applying the algorithm to the long-period channels of all data archived at the IRIS DMC from 1976 to 2005 for the seismic phases S and P (from 23° to 100°) and SS and PP (from 50° to 170°). The resulting large data sets are inverted along with existing surface wave and updated differential traveltime measurements for new mantle models of S and P velocity. The resolution of the new model is enhanced, particularly, in the mid-mantle where SS and PP turn. We find that slow anomalies in the central Pacific and Africa extend from the core–mantle boundary to the upper mantle, but their direct connection to surface hotspots is beyond our resolution. Furthermore, we find that fast anomalies that are likely associated with subducting slabs disappear between 1700 and 2500 km, and thus are not continuous features from the upper to lower mantle despite our extensive coverage and high resolution of the mid-mantle. 相似文献
This technique can be applied to any seismic phase that is observed over a reasonably large distance range. Here, we present the results of applying the algorithm to the long-period channels of all data archived at the IRIS DMC from 1976 to 2005 for the seismic phases S and P (from 23° to 100°) and SS and PP (from 50° to 170°). The resulting large data sets are inverted along with existing surface wave and updated differential traveltime measurements for new mantle models of S and P velocity. The resolution of the new model is enhanced, particularly, in the mid-mantle where SS and PP turn. We find that slow anomalies in the central Pacific and Africa extend from the core–mantle boundary to the upper mantle, but their direct connection to surface hotspots is beyond our resolution. Furthermore, we find that fast anomalies that are likely associated with subducting slabs disappear between 1700 and 2500 km, and thus are not continuous features from the upper to lower mantle despite our extensive coverage and high resolution of the mid-mantle. 相似文献
4.
5.
Chihi Hayet Galli Alain Ravenne Christian Tesson Michel de Marsily Ghislain 《Natural Resources Research》2000,9(1):77-95
The object of this study is to build a three-dimensional (3D) geometric model of the stratigraphicunits of the margin of the Rhone River on the basis of geophysical investigations by a networkof seismic profiles at sea. The geometry of these units is described by depth charts of eachsurface identified by seismic profiling, which is done by geostatistics. The modeling starts bya statistical analysis by which we determine the parameters that enable us to calculate thevariograms of the identified surfaces. After having determined the statistical parameters, wecalculate the variograms of the variable Depth. By analyzing the behavior of the variogramwe then can deduce whether the situation is stationary and if the variable has an anisotropicbehavior. We tried the following two nonstationary methods to obtain our estimates: (a) Themethod of universal kriging if the underlying variogram was directly accessible. (b) Themethod of increments if the underlying variogram was not directly accessible. After havingmodeled the variograms of the increments and of the variable itself, we calculated the surfacesby kriging the variable Depth on a small-mesh estimation grid. The two methods then arecompared and their respective advantages and disadvantages are discussed, as well as theirfields of application. These methods are capable of being used widely in earthsciences forautomatic mapping of geometric surfaces or for variables such as a piezometricsurface or aconcentration, which are not stationary, that is, essentially, possess a gradient or a tendencyto develop systematically in space. 相似文献
6.
Thorsten W. Becker Göran Ekström Lapo Boschi John H. Woodhouse 《Geophysical Journal International》2007,171(1):451-462
We measure the degree of consistency between published models of azimuthal seismic anisotropy from surface waves, focusing on Rayleigh wave phase-velocity models. Some models agree up to wavelengths of ∼2000 km, albeit at small values of linear correlation coefficients. Others are, however, not well correlated at all, also with regard to isotropic structure. This points to differences in the underlying data sets and inversion strategies, particularly the relative 'damping' of mapped isotropic versus anisotropic anomalies. Yet, there is more agreement between published models than commonly held, encouraging further analysis. Employing a generalized spherical harmonic representation, we analyse power spectra of orientational (2Ψ) anisotropic heterogeneity from seismology. We find that the anisotropic component of some models is characterized by stronger short-wavelength power than the associated isotropic structure. This spectral signal is consistent with predictions from new geodynamic models, based on olivine texturing in mantle flow. The flow models are also successful in predicting some of the seismologically mapped patterns. We substantiate earlier findings that flow computations significantly outperform models of fast azimuths based on absolute plate velocities. Moreover, further evidence for the importance of active upwellings and downwellings as inferred from seismic tomography is presented. Deterministic estimates of expected anisotropic structure based on mantle flow computations such as ours can help guide future seismologic inversions, particularly in oceanic plate regions. We propose to consider such a priori information when addressing open questions about the averaging properties and resolution of surface and body wave based estimates of anisotropy. 相似文献
7.
Bradford E. Prather Ciaran O'Byrne Carlos Pirmez Zoltán Sylvester 《Basin Research》2017,29(3):394-416
Deciphering the role slope topography plays in partitioning sediment on siliciclastic continental slope and base‐of‐slope systems helps our understanding of slope depositional processes in significant ways: (1) by validation of large‐scale depositional process models for continental margins, (2) by validation of numerical basin‐scale stratigraphic forward models used to test and deploy source‐to‐sink (S2S) concepts and (3) by creating models for setting reservoir presence and quality expectations in frontier areas poorly constrained by wells and seismic. A global database consisting of >700 km of drilled stratigraphy provide empirical rock data lacking from most S2S studies. Analysis of calibrated seismic stratigraphic units characterised using the contextual framework laid out in this paper show that both gross depositional environments (GDEs) and sand content occur across slope profiles in systematic ways. The challenge in using these observations to quantify reservoir risk and uncertainty lies with relating the observations to depositional processes that can be used to characterise frontier basins that lack calibration. Depositional process‐based understanding encoded in 3D stratigraphic forward models (SFM) can simulate both lithologies and GDEs providing broad predictions for exploration at the scale of an entire basin or slope system. Stratigraphic forward models allow the integration of S2S understanding and provide a framework for testing sediment‐partitioning hypotheses in frontier settings. Valid S2S models must balance sediment yield from the source catchments with sinks, and be consistent with basin specific observations. The proportions of GDEs across the slope provide additional validation criteria to ensure the models are plausible. 相似文献
8.
Analysis of seismic data suggests that hydrocarbon deposits are often associated with higher than usual values of attenuation, but this is generally ignored during amplitude-versus-offset (AVO) analysis. The effect can be modelled with equivalent medium theory based on the squirt flow concept, but the excess attenuation is associated with strong velocity dispersion. Consequently, when we study reflections from the interface between such an equivalent medium and an elastic overburden we find that the reflection coefficient varies with frequency. The impact of this variation depends on the AVO behaviour at the interface; class I reflections tend to be shifted to higher frequency while class III reflections have their lower frequencies amplified. We calculate synthetic seismograms for typical models using the reflectivity method for materials with frequency dependent velocities and attenuations, and find that these effects are predicted to be detectable on stacked data. Two field data sets show frequency anomalies similar to those predicted by the analysis, and we suggest that our modelling provides a plausible explanation of the observations. 相似文献
9.
The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges. 相似文献
10.
Measurement of dispersed vitrinite reflectance in organic sediments is one of the few regional data sets used for placing bounds on the thermal history of a sedimentary basin. Reflectance data are important when access to complementary information such as high‐quality seismic data is unavailable to place bounds on subsidence history and in locations where uplift is an important part of the basin history. Attributes which make vitrinite reflectance measurements a useful data set are the relative ease of making the measurement, and the availability of archived well cores and cuttings in state, provincial, and federal facilities. In order to fully utilize vitrinite data for estimating the temperature history in a basin, physically based methods are required to calibrate an equivalent reflectance from a modelled temperature history with measured data. The most common method for calculating a numerical vitrinite reflectance from temperature history is the EASY%Ro method which we show systematically underestimates measured data. We present a new calculated reflectance model and an adjustment to EASY%Ro which makes the correlation between measured vitrinite values and calculated vitrinite values a physical relationship and more useful for constraining thermal models. We then show that calibrating the thermal history to vitrinite on a constant age date surface (e.g., top Cretaceous) instead of calibrating the thermal history in depth removes the heating rate component from the reflectance calculation and makes thermal history calibration easier to understand and more directly related to heat flow. Finally, we use bounds on the vitrinite–temperature relationships on a constant age date surface to show that significant uncertainty exists in the vitrinite data reported in most data sets. 相似文献
11.
A self-adaptive automated parametrization approach is suggested for the sequential inversion of controlled-source seismic tomography and gravity data. The velocities and interfaces are parametrized by their Haar wavelet expansion coefficients. Only those coefficients that are well constrained by the data, as measured by the number of rays that cross the corresponding wavelet function support area and their angular coverage, are inverted for, others are set to zero. This approach results in a reasonable distribution of resolution throughout the model even in cases of irregular ray coverage and does overcome the trade-off between different types of model parameters. A modified sequential inversion approach is suggested to join the traveltimes and gravity anomalies inversion. An algorithm is developed that inverts for smooth velocity and density variations inside the seismic layer, the position of its bottom interface as well as for optimal values of the velocity-to-density regression coefficients. The algorithm makes use of direct (diving), reflected and head (critically refracted) wave traveltimes. The algorithm workflow is demonstrated on a synthetic data example. 相似文献
12.
Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements 总被引：28，自引：0，他引：28
G. D. Bensen M. H. Ritzwoller M. P. Barmin A. L. Levshin F. Lin M. P. Moschetti N. M. Shapiro Y. Yang 《Geophysical Journal International》2007,169(3):1239-1260
Ambient noise tomography is a rapidly emerging field of seismological research. This paper presents the current status of ambient noise data processing as it has developed over the past several years and is intended to explain and justify this development through salient examples. The ambient noise data processing procedure divides into four principal phases: (1) single station data preparation, (2) cross-correlation and temporal stacking, (3) measurement of dispersion curves (performed with frequency–time analysis for both group and phase speeds) and (4) quality control, including error analysis and selection of the acceptable measurements. The procedures that are described herein have been designed not only to deliver reliable measurements, but to be flexible, applicable to a wide variety of observational settings, as well as being fully automated. For an automated data processing procedure, data quality control measures are particularly important to identify and reject bad measurements and compute quality assurance statistics for the accepted measurements. The principal metric on which to base a judgment of quality is stability, the robustness of the measurement to perturbations in the conditions under which it is obtained. Temporal repeatability, in particular, is a significant indicator of reliability and is elevated to a high position in our assessment, as we equate seasonal repeatability with measurement uncertainty. Proxy curves relating observed signal-to-noise ratios to average measurement uncertainties show promise to provide useful expected measurement error estimates in the absence of the long time-series needed for temporal subsetting. 相似文献
13.
Guofeng Cao Phaedon C. Kyriakidis Michael F. Goodchild 《International journal of geographical information science》2013,27(11):1773-1791
Categorical spatial data, such as land use classes and socioeconomic statistics data, are important data sources in geographical information science (GIS). The investigation of spatial patterns implied in these data can benefit many aspects of GIS research, such as classification of spatial data, spatial data mining, and spatial uncertainty modeling. However, the discrete nature of categorical data limits the application of traditional kriging methods widely used in Gaussian random fields. In this article, we present a new probabilistic method for modeling the posterior probability of class occurrence at any target location in space-given known class labels at source data locations within a neighborhood around that prediction location. In the proposed method, transition probabilities rather than indicator covariances or variograms are used as measures of spatial structure and the conditional or posterior (multi-point) probability is approximated by a weighted combination of preposterior (two-point) transition probabilities, while accounting for spatial interdependencies often ignored by existing approaches. In addition, the connections of the proposed method with probabilistic graphical models (Bayesian networks) and weights of evidence method are also discussed. The advantages of this new proposed approach are analyzed and highlighted through a case study involving the generation of spatial patterns via sequential indicator simulation. 相似文献
14.
Mohammad Reza Tarybakhsh Mehdi Assareh Mohammad Taghi Sadeghi Ali Ahmadi 《Natural Resources Research》2018,27(4):517-529
Determination of gas–oil minimum miscibility conditions is one of the important design parameters to improve the displacement efficiency of the hydrocarbon reservoir during enhanced oil recovery with gas injection. In this work, a support vector regression (SVR) model is developed using experimental data to estimate the minimum miscibility pressure (MMP) for various reservoir fluids and injection gases. Experimental MMP data taken from the reliable literature were used as input. Each data point input includes methane and intermediate components mole percent, plus fraction properties and reservoir temperature related to reservoir fluid and CO2, H2S, N2 and intermediate mole fractions, and intermediate properties of the injected gas. Experimental MMP is regarded as the model output. The database contains 135 datasets, from which 125 datasets were used for model development, and the rest were used for model evaluation. Genetic algorithm was implemented to optimize the SVR model parameters. The proposed data-driven model was verified by statistical validation data. The model results illustrate a correlation coefficient (R2) of 0.999. In addition, the SVR results demonstrate the proposed model to be a fast tool and a robust approach to map input space to output features. The SVR model was compared to popular data-driven MMP estimation models as well. This comparison presents an acceptable accuracy relative to this estimation model. Finally, the presented model was evaluated against a comprehensive theoretical model of slim tube compositional simulation on a trusted literature dataset. 相似文献
15.
We analyse earthquakes recorded at The Geysers geothermal field in California, an area where industrial activity induces seismicity. The seismicity is characterized by the seismic b -value and D , the fractal dimension of earthquake hypocentres measured from sliding windows containing 200 events. We study a group of events strongly clustered around an injection well. Over most of the time period examined we find a positive correlation between b and D . However, during the initiation of injection into a new well we find instead a negative correlation. The differences in correlation are statistically significant at the 1 σ level but only marginally so at the 2 σ level. These results provide evidence for a transient change in the seismic mechanisms operating, and may be explained by a change from conditions of slow stress loading to rapid loading as a result of the build-up of the rate of water injection into the reservoir. 相似文献
16.
S. A. Clark E. Glorstad‐Clark J. I. Faleide D. Schmid E. H. Hartz W. Fjeldskaar 《Basin Research》2014,26(4):550-566
We present tectonic models of progressive basin formation in the south‐west Barents Sea derived as part of the PETROBAR project (Petroleum‐related studies of the Barents Sea region). The basin architecture developed as a multi‐stage rift preceding the creation of the sheared/transtensional margin conjugate to NE Greenland. N‐ to NNE‐striking basins, with sediment thicknesses in places exceeding 15 km, are separated by basement highs. We use two basin analysis approaches, BMT? backstripping and TecMod?time‐forward modelling, to determine stretching factors through time along the profile PETROBAR‐07. This 550 km‐long profile derived from wide‐angle reflection/refraction seismic data acquired in 2007, coincident with deep multichannel seismic reflection data. Detailed stratigraphic analysis of the reflection profile, in concert with a dense grid of 2D profiles tied to wells, provides timing and water depth constraint for the models. Velocity analysis of the wide‐angle data provides constraint on the cumulative crustal stretching. The north‐west trending cross‐section extends from continental craton, at the Varanger Peninsula, to within 16 km of the interpreted continent–ocean boundary. Rifting along the profile was episodic, with four distinct phases of basin formation during the Carboniferous, the Late Permian–Triassic, the Late Jurassic–Early Cretaceous and the Late Cretaceous–Eocene. Collectively, the basins exhibit a general trend of younging, narrowing, and deepening oceanward, suggesting a gradual focusing of rifting prior to final breakup. Cumulative stretching factors derived from BMT and TecMod correlate well with observed crustal thinning, and the two models provide uncertainty bounds for stretching factors for the separate rift phases. In contrast to orthogonally rifted margins, stretching is relatively minor immediately prior to transform breakup, with greater stretching occurring during earlier rift phases. 相似文献
17.
Drumond David Alvarenga Rolo Roberto Mentzingen Costa João Felipe Coimbra Leite 《Natural Resources Research》2019,28(1):145-152
Experimental variograms are crucial for most geostatistical studies. In kriging, for example, the variography has a direct influence on the interpolation weights. Despite the great importance of variogram estimators in predicting geostatistical features, they are commonly influenced by outliers in the dataset. The effect of some randomly spatially distributed outliers can mask the pattern of the experimental variogram and produce a destructuration effect, implying that the true data spatial continuity cannot be reproduced. In this paper, an algorithm to detect and remove the effect of outliers in experimental variograms using the Mahalanobis distance is proposed. An example of the algorithm’s application is presented, showing that the developed technique is able to satisfactorily detect and remove outliers from a variogram.
相似文献18.
冰川冰储量不仅是冰川的重要属性,而且是核算冰川水资源及预测冰川变化的基础数据,因此准确计算冰川冰储量及其变化具有重要的理论与现实意义。目前冰川储量估算的主要方法有经验公式法、冰厚模型估算法、探地雷达法;冰川储量相对变化计算方法有实地测量法和遥感监测法。通过系统分析和讨论各计算方法的原理、现状及存在的问题,以期为冰川储量估算提供方法参考。研究表明：对于冰川冰储量计算而言,经验公式法适用于区域性或全球性的冰川储量估算;模型估算法适用于个体或小范围冰川储量估算;探地雷达法适用于人类易到达区域冰川储量的估算。对于冰川冰储量相对变化计算,实地测量法适用于对精度要求高且满足实地测量条件的单条或中小型冰川,遥感监测法适用于全球性冰储量变化估算,但需改进算法和提高数据空间分辨率。目前,随着无人机技术的逐步应用,以及冰川流速等理论模型的提出,为冰川冰储量估算方法的发展提供了新契机。 相似文献
19.
Vanessa Cerqueira Koppe João Felipe Coimbra Leite Costa Jair Carlos Koppe Fernando Gambin Gary Fallon Nick Davies 《Natural Resources Research》2007,16(4):293-303
Seismic reflection methods measure the time a seismic wave takes to travel through the ground, from the user defined source
to a series of signal monitoring sensors known as geophones. The measured times need to be depth converted to allow for integration
with other geological data. In order to convert from time to depth, an estimate of the rock volume velocity field must be
made. The velocity field estimate can be made by assignment of velocity estimates to a geological model independent of the
seismic processing. This article presents the results of using the acoustic geophysical log data extrapolated via sequential
Gaussian simulation to derive the velocity field. The uncertainties associated with the velocity estimates were significant
and provided the means to assess confidence limits for the actual depth determination. The technique is assessed by application
to a major coal deposit, approximately 2.1 m thick and 210 m deep. Considering only the uncertainty associated with estimating
the velocity field, half of the confidence interval values showed approximately 1 m of uncertainty in depth. The application
of sequential Gaussian simulation to model the 3D distribution of acoustic velocity can be extended to other geophysical log
parameters or derived estimates. 相似文献
20.
Surface-wave polarization data and global anisotropic structure 总被引：1，自引：0，他引：1
In the past few years, seismic tomography has begun to provide detailed images of seismic velocity in the Earth's interior which, for the first time, give direct observational constraints on the mechanisms of heat and mass transfer. The study of surface waves has led to quite detailed maps of upper-mantle structure, and the current global models agree reasonably well down to wavelengths of approximately 2000 km. Usually, the models contain only elastic isotropic structure, which provides an excellent fit to the data in most cases. For example, the variance reduction for minor and major arc phase data in the frequency range 7–15 mHz is typically 65–92 per cent and the data are fit to within 1–2 standard deviations. The fit to great-circle phase data, which are not subject to bias from unknown source or instrument effects, is even better. However, there is clear evidence for seismic anisotropy in various places on the globe. This study demonstrates how much (or little) the fit to the data is improved by including anisotropy in the modelling process. It also illuminates some of the trade-offs between isotropic and anisotropic structure and gives an estimate of how much bias is introduced by neglecting anisotropy. Finally, we show that the addition of polarization data has the potential for improving recovery of anisotropic structure by diminishing the trade-offs between isotropic and anisotropic effects. 相似文献