Approximate traveltime inversion in downhole microseismic monitoring     

Sergey Yaskevich  Anton A. Duchkov  Yuriy Ivanov 《Geophysical Prospecting》2020,68(3):918-925

In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media.  相似文献   

Saltwater Upconing Due to Cyclic Pumping by Horizontal Wells in Freshwater Lenses          下载免费PDF全文

Pieter S. Pauw  Sjoerd E.A.T.M. van der Zee  Anton Leijnse  Gualbert H.P. Oude Essink 《Ground water》2016,54(4):521-531

This article deals with the quantification of saltwater upconing below horizontal wells in freshwater lenses using analytical solutions as a computationally fast alternative to numerical simulations. Comparisons between analytical calculations and numerical simulations are presented regarding three aspects: (1) cyclic pumping; (2) dispersion; and (3) finite horizontal wells in a finite domain (a freshwater lens). Various hydrogeological conditions and pumping regimes within a dry half year are considered. The results show that the influence of elastic and phreatic storage (which are not taken into account in the analytical solutions) on the upconing of the interface is minimal. Furthermore, the analytical calculations based on the interface approach compare well with numerical simulations as long as the dimensionless interface upconing is below 1/3, which is in line with previous studies on steady pumping. Superimposing an analytical solution for mixing by dispersion below the well over an analytical solution based on the interface approach is appropriate in case the vertical flow velocity around the interface is nearly constant but should not be used for estimating the salinity of the pumped groundwater. The analytical calculations of interface upconing below a finite horizontal well compare well with the numerical simulations in case the distance between the horizontal well and the initial interface does not vary significantly along the well and in case the natural fluctuation of the freshwater lens is small. In order to maintain a low level of salinity in the well during a dry half year, the dimensionless analytically calculated interface upconing should stay below 0.25.  相似文献   

Isotopic Composition of Waters in the Dem?novská dolina Valley and Its Underground Hydrologic System During Winter and Spring of 2010/2011          下载免费PDF全文

Peter MALIK  Dagmar HAVIAROVA  Anton AUXT  Zuzana GROLMUSOVA  Juraj MICHALKO 《地下水科学与工程》2013,1(2):14-23

In the Dem?novská dolina Cave system (Slovakia) and its vicinity, 32 sampling places for regular observation (in 2-months interval) of δ18O and δ2H in water were established. This monitoring included precipitation waters, waters in the surface streams, waters of the underground hydrological system as well as the dripping seepage waters of the cave system. Altitudinal extent of the area was from 800 m a.s.l. (lowermost cave entrance) to 2024 m a.s.l. (Chopok Mt. on the top of the crystalline range). Initial results show some similarities but also differences within the analyzed water types. For precipitation, a high variability of isotopic composition was confirmed, from quite depleted up to more enriched waters (δ18O from -16.8‰ up to -5.7‰; δ2H from -121.6‰ to -32.7‰). During the recharge process and groundwater/surface water formation, precipitation water is homogenized, what is reflected in much more stable isotope content. The most depleted (δ18O ≈ -11.7‰ to -10.8‰; δ2H ≈ -78.9‰ to -73.4‰) were the waters of surface streams, running from the northern slopes of the Nízke Tatry Mts., formed by crystalline rocks, alochtonous to the under?ground hydrological system. Smaller autochtonous surface water streams (formed in the side valleys of the main karstic canyon) are slightly enriched (heavier, as δ18O ≈ -11.4‰ to -10.6‰; δ2H ≈ -78.3‰ to -71.5‰), what reflects lower altitudes of their watersheds. Interesting is the distribution of the isotope content of the underground streams in the cave system. The most depleted are the underground streams directly (visibly) communicating with surface waters (δ18O≈-11.33±0.13‰; δ2H≈-76.88±1.68‰). Extent of the relationship of underground streams to the autochtonous seepage waters (slow circulation through the fissures) is manifested by respective portion of iso?topically enriched waters–as the underground streams show different isotope composition. The combination of the alochtonous water components (from surface streams reaching the karstic area from the adjacent crystalline via swallow holes) and autochtonous water components (recharged directly in karstified limestones) is visible especially on the subsurface stream of Dem?novka. The most isotopically enriched (heaviest) of all water types are dripping seepage waters (δ18O ≈ -10.4‰ to -9.4‰; δ2H ≈ -71.6‰ to -65.0‰).  相似文献   

Single-crystal elasticity of stishovite: New experimental data obtained using high-frequency resonant ultrasound spectroscopy and a Gingham check structure model     

Akira Yoneda  Titus Cooray  Anton Shatskiy 《Physics of the Earth and Planetary Interiors》2012

Single-crystal elasticity of stishovite was examined using a new experimental technique and an empirical macroscopic model. Employing high-frequency resonant ultrasound spectroscopy, single-crystal elastic constants of stishovite were determined: C₁₁ = 443(3), C₃₃ = 781(4), C₁₂ = 193(2), C₂₃ = 199(2), C₄₄ = 256(2), and C₃₃ = 316(2) GPa. The frequency range of the resonant ultrasound spectroscopy techniques was 6–20 MHz, which is much lower than the ～10 GHz range of the Brillouin scattering technique. Of the elastic constants, the shear elastic constants C₄₄ and C₆₆ are consistent with the average value of the previously mentioned Brillouin scattering. Conversely, the four elastic constants, C₁₁, C₃₃, C₁₂, and C₂₃, slightly deviate outside the range of previous Brillouin scattering results. The present results, except those for C₁₂, are consistent with recent lattice dynamic analysis of inelastic X-ray scattering data. The adiabatic bulk modulus was calculated as 298 GPa, which is smaller and more consistent with the result of compression experiments than any other Brillouin scattering results (301–312 GPa). The present result shows greater P-wave velocity anisotropy (24.7%) than any preceding work. To understand the unique elastic properties of stishovite, the Gingham check model was proposed and examined. The result shows that the octahedron of 6-coordinated Si in stishovite crystal has stiffness comparable to that of diamond.  相似文献   

Prestack shot-gather depth migration by a rigid flow of Gaussian wave packets     

Anton A. Duchkov  Fredrik Andersson  Rikard Ojala 《Studia Geophysica et Geodaetica》2012,56(1):83-106

We introduce a new method for prestack depth migration of seismic common-shot gathers. The computational procedure follows standard steps of the reverse-time migration, i.e., downward continuation of the source and the receiver wavefields, followed by application of an imaging condition (e.g. zero-lag cross-correlation of these fields). In our method we first find a sparse data representation with a small number of Gaussian wave packets. We then approximate the downward wavefield propagation (for the source and the receiver fields) by a rigid flow of these wave packets along seismic rays. In this case, the wave packets are simply translated and rotated according to the ray geometry. One advantage of using Gaussian wave packets is that analytic formulas can be used for translation, rotation, and the application of the cross-correlation imaging condition. Moreover, they allow more sparse representations than competing methods. Finally we formulate a computationally and memory efficient migration procedure, as only few rays have to be traced, and since it is cheap to compute the cross-correlation for the intersecting rays.  相似文献   

Fractal characteristics of the ULF emissions along a meridian profile, based on the 210 MM stations data     

Anton Varlamov  Natalia Smirnova  Masashi Hayakawa  Kiyohumi Yumoto 《Acta Geophysica》2012,60(3):928-941

Fractal analysis of magnetic records (1 Hz sampling rate) of 5 stations (Guam, Moshiri, Paratunka, Magadan, and Chokurdakh) located along the 210 magnetic meridian (210 MM) has been performed using the Higuchi method. The period of 22 months (October 1992 to July 1994) that embodies the date of the strong Guam earthquake of 8 August1993 has been considered. A comparison of the ULF emissions scaling parameters (spectral exponents β and fractal dimensions D) obtained at different latitudes has been made. Dependence of β and D on the Kp index of geomagnetic activity has been analyzed for each of the 24 local time intervals. It is revealed that D decreases ( β increases) with increasing geomagnetic activity at all stations, but the rates of decrease (increase) are different at different stations and in different time intervals. It is shown that the evening, night and early morning hours are preferable to study magnetospheric effects, whereas the noon hours are the most suitable for the analysis of lithospheric effects. A possibility of using the data of the 210 MM stations as reference materials for the Guam seismically active area is discussed.  相似文献   

Review of vibroseis data acquisition and processing for better amplitudes: adjusting the sweep and deconvolving for the time-derivative of the true groundforce   总被引：1，自引：0，他引：1  

Anton Ziolkowski 《Geophysical Prospecting》2010,58(1):41-54

The goal of vibroseis data acquisition and processing is to produce seismic reflection data with a known spatially-invariant wavelet, preferably zero phase, such that any variations in the data can be attributed to variations in geology. In current practice the vibrator control system is designed to make the estimated groundforce equal to the sweep and the resulting particle velocity data are cross-correlated with the sweep. Since the downgoing far-field particle velocity signal is proportional to the time-derivative of the groundforce, it makes more sense to cross-correlate with the time-derivative of the sweep. It also follows that the ideal amplitude spectrum of the groundforce should be inversely proportional to frequency. Because of non-linearities in the vibrator, bending of the baseplate and variable coupling of the baseplate to the ground, the true groundforce is not equal to the pre-determined sweep and varies not only from vibrator point to vibrator point but also from sweep to sweep at each vibrator point. To achieve the goal of a spatially-invariant wavelet, these variations should be removed by signature deconvolution, converting the wavelet to a much shorter zero-phase wavelet but with the same bandwidth and signal-to-noise ratio as the original data. This can be done only if the true groundforce is known. The principle may be applied to an array of vibrators by employing pulse coding techniques and separating responses to individual vibrators in the frequency domain. Various approaches to improve the estimate of the true groundforce have been proposed or are under development; current methods are at best approximate.  相似文献   

Modelling ice conditions in the easternmost Gulf of Finland in the Baltic Sea     

Vladimir Ryabchenko  Anton Dvornikov  Jari Haapala  Kai Myrberg 《Continental Shelf Research》2010

To model ice conditions in the eastern Gulf of Finland, a high-resolution three-dimensional hydrodynamic model is coupled with the advanced sea-ice model HELMI (Haapala et al., 2005). To test the model in extreme situations, the ice pattern in the eastern Gulf of Finland was simulated for a mild ice winter (2007–2008) and for a moderate one (2003–2004). The reference runs were performed on the assumption that the ice in the model domain is fast ice if the sea depth is less than 10 m. Using this assumption, the ice thickness averaged over the Neva Bay (the easternmost part of the Gulf of Finland) is overestimated by the model for almost the entire wintertime in the mild winter and during the ice formation and melting periods in the moderate winter, as compared with the thickness reported in ice charts.  相似文献   

A semantic annotation tool for hydrologic sciences     

Michael Piasecki  Bora Beran 《Earth Science Informatics》2009,2(3):157-168

  相似文献   

Green space planning and land use: An assessment of urban regional and green structure planning in Greater Copenhagen     

Ole H. Caspersen MSc  PhD.  Cecil C. Konijnendijk  Anton S. Olafsson 《Geografisk tidskrift / udgivet af Bestyrelsen for Det Kongelige danske geografiske selskab》2013,113(2):7-20

Geografisk Tidsskrift, Danish Journal of Geography 106(2):7–20, 2006

In its Regional Plan of 2005, The Greater Copenhagen Authority (abbreviated as “HUR” in Danish) places special emphasis on the future recreational values associated with the regional green structure. In this paper, the development of the urban green structure in Greater Copenhagen is elaborated upon, focusing on land use changes and the effectiveness of regionally coordinated planning measures. EU MOLAND data are used to analyse the development of the region's green structure during the period 1954 to 1998. Analysis of two “green wedges” within the green structure illustrates that the development of the green recreational areas is the result of both formal and more informal planning initiatives. Development has shown equal phases that correspond to the applied regional planning measures and the general economic conditions. However, local preferences in the involved municipalities likewise have played an important role and have resulted in different urbanisation pressure within the two wedges. Land use has transformed from an agricultural to a primarily recreational landscape. In some areas, however, urbanisation pressure has resulted in residential and green industrial areas instead of the planned recreational land use. Based on its historical development it is concluded that future preservation and development of the green structure in Greater Copenhagen requires regional planning measures to be incorporated into municipal plans. In this way the increasingly independent municipalities will comply to objectives of the Regional Plan 2005.  相似文献