Directional Acoustic-Velocity Distribution in a Petroleum Reservoir

1975 ◽  
Vol 15 (03) ◽  
pp. 234-246
Author(s):  
L.Z. Shuck ◽  
K-H Frohne ◽  
G.E. Rennick ◽  
D.M. Evans

Abstract Experimental conditions, theoretical considerations, and results of adirectional acoustic-velocity survey in the plane of formation of an oilreservoir in Pennsylvania are given. More than 60 blasting-cap acoustic sourceswere discharged opposite the Bradford Third sand in 10 different wells locatedwithin a 1,200- by 1,350-ft rectangular array of 41 wells. Each source wasmonitored simultaneously by 11 different hydrophones strategically locatedwithin the 41-well array. Both acoustic sources and sensors were located at theformation elevation, about 2,000 ft below the surface, so that the in-situvelocities determined were in the plane of the sandstone formation. Well-to-well acoustic velocities were calculated and are presented in tabularand graphical form. The velocity data were further analyzed statistically toyield over-all reservoir velocity trends. Results of the study show that thereservoir is acoustically heterogeneous, with local variations in P-wavevelocity of as much as 20 percent. The over-all P-wave velocity trend for thereservoir was greater in the east-west direction than in the north-southdirection by about 700 ft/sec. The mean acoustic velocity for all wells in alldirections was 15,621 ft/ sec. The acoustic-velocity distribution is also correlated to directionalpermeability data obtained by oriented core analyses. Introduction The directional variation of geophysical properties affecting wave propagationhas long been known. A brief examination of the history and evolutionaryprocesses of sedimentary rock should leave no doubt about the nature of theformations. However, little attention has been given to directional propertiesin the evaluation and recovery of natural resources. Evaluation of thedirectional acoustic (seismic) velocity distribution in reservoirs is perhapsthe most promising means of evaluating directional reservoir properties. Themagnitude of directional variation can be predicted to a degree by consideringthe lithology and the geological features of an area. However, there are alwaysuncertainties in the study of wave propagation through the earth; numerousvariables influencing wave-propagation velocities have been determined. This research was initiated originally to aid in acoustically mapping thegrowth of hydraulically induced fractures in time and space, and to assist indelineating existing fracture systems. The reservoir investigated had been thesubject of other research programs, and considerable information was available.Using the available information, however, large directional variations inacoustic velocities over small areas would not necessarily have beenanticipated. As will be shown, the velocity variations were substantial andwere correlatable, to a degree, with other geophysical and reservoirproperties. Since the late 1940's and early 1950's, a tremendous amount of work has beendone in active and passive seismic monitoring and in the general area of wavepropagation through geological materials. Most of the work falls into one offour general categories:active seismic exploration,passive earthquakeor nuclear undergroundexplosion monitoring,active oil and gas welllogging, androck noise and slope-stability studies. Obviously, eachcategory directly or indirectly involves the directional nature of propagationvelocities. The first two categories, and usually the fourth, are associatedwith very low frequency monitoring systems, such as geophones usually operatingin the 0- to 50-hz range. Frequencies used in Category 3 are usually in the20-khz range since the waves typically are transmitted only 3 ft verticallyfrom a transmitter to a receiver.

Geophysics ◽  
2016 ◽  
Vol 81 (4) ◽  
pp. D293-D301 ◽  
Author(s):  
Fuyong Yan ◽  
De-hua Han ◽  
Qiuliang Yao ◽  
Xue-Lian Chen

We have made various experimental investigations on rock-physics properties of the halite salt coming from a salt dome in the U.S. Gulf Coast Basin. The effect of crystal defects and intercrystal cracks on the P-wave velocity of the halite salt sample can be mitigated after high-pressure annealing. The temperature effect on seismic velocities of halite salt is dominant relative to the stress effect. Azimuthal anisotropy is not observed on the halite salt samples. We have observed that the velocity variations in different directions were mostly caused by the crystal-scale heterogeneity. For the salt structures primarily made of creep-deformed halite, the effect of seismic velocity anisotropy might be negligible for seismic exploration. No significant dispersion of seismic velocities was observed from the low-frequency measurement. Our measurements supply basic information for seismic velocity model building that may help to improve seismic imaging of salt structures.


2020 ◽  
Vol 25 (3) ◽  
pp. 415-423
Author(s):  
Ahmed Lachhab ◽  
El Mehdi Benyassine ◽  
Mohamed Rouai ◽  
Abdelilah Dekayir ◽  
Jean C. Parisot ◽  
...  

The tailings of Zeida's abandoned mine are found near the city of Midelt, in the middle of the high Moulouya watershed between the Middle and the High Atlas of Morocco. The tailings occupy an area of about 100 ha and are stored either in large mining pit lakes with clay-marl substratum or directly on a heavily fractured granite bedrock. The high contents of lead and arsenic in these tailings have transformed them into sources of pollution that disperse by wind, runoff, and seepage to the aquifer through faults and fractures. In this work, the main goal is to identify the pathways of contaminated water with heavy metals and arsenic to the local aquifers, water ponds, and Moulouya River. For this reason, geophysical surveys including electrical resistivity tomography (ERT), seismic refraction tomography (SRT) and very low-frequency electromagnetic (VLF-EM) methods were carried out over the tailings, and directly on the substratum outside the tailings. The result obtained from combining these methods has shown that pollutants were funneled through fractures, faults, and subsurface paleochannels and contaminated the hydrological system connecting groundwater, ponds, and the river. The ERT profiles have successfully shown the location of fractures, some of which extend throughout the upper formation to depths reaching the granite. The ERT was not successful in identifying fractures directly beneath the tailings due to their low resistivity which inhibits electrical current from propagating deeper. The seismic refraction surveys have provided valuable details on the local geology, and clearly identified the thickness of the tailings and explicitly marked the boundary between the Triassic formation and the granite. It also aided in the identification of paleochannels. The tailings materials were easily identified by both their low resistivity and low P-wave velocity values. Also, both resistivity and seismic velocity values rapidly increased beneath the tailings due to the compaction of the material and lack of moisture and have proven to be effective in identifying the upper limit of the granite. Faults were found to lie along the bottom of paleochannels, which suggest that the locations of these channels were caused by these same faults. The VLF-EM surveys have shown tilt angle anomalies over fractured areas which were also evinced by low resistivity area in ERT profiles. Finally, this study showed that the three geophysical methods were complementary and in good agreement in revealing the pathways of contamination from the tailings to the local aquifer, nearby ponds and Moulouya River.


2020 ◽  
Vol 14 (3) ◽  
pp. 7082-7093
Author(s):  
Jahirwan Ut Jasron ◽  
Sudjito Soeparmani ◽  
Lilis Yuliati ◽  
Djarot B. Darmadi

The hydrodynamic performance of oscillating water column (OWC) depends on the depth of the water, the size of the water column and its arrangement, which affects the oscillation of the water surface in the column. An experimental method was conducted by testing 4 water depths with wave periods of 1-3 s. All data recorded by the sensor is then processed and presented in graphical form. The research focused on analyzing the difference in wave power absorption capabilities of the three geometric types of OWC based on arrangements of water columns. The OWC devices designed as single water column, the double water column in a series arrangement which was perpendicular to the direction of wave propagation, and double water column in which the arrangement of columns was parallel to the direction of wave propagation. This paper discussed several factors affecting the amount of power absorbed by the device. The factors are the ratio of water depth in its relation to wavelength (kh) and the inlet openings ratio (c/h) of the devices. The test results show that if the water depth increases in the range of kh 0.7 to 0.9, then the performance of the double chamber oscillating water column (DCOWC) device is better than the single chamber oscillating water column (SCOWC) device with maximum efficiency for the parallel arrangement 22,4%, series arrangement 20.8% and single column 20.7%. However, when referring to c/h, the maximum energy absorption efficiency for a single column is 27.7%, double column series arrangement is 23.2%, and double column parallel arrangement is 29.5%. Based on the results of the analysis, DCOWC devices in parallel arrangement showed the ability to absorb better wave power in a broader range of wave frequencies. The best wave of power absorption in the three testing models occurred in the wave period T = 1.3 seconds.


Author(s):  
Fumiaki Nagashima ◽  
Hiroshi Kawase

Summary P-wave velocity (Vp) is an important parameter for constructing seismic velocity models of the subsurface structures by using microtremors and earthquake ground motions or any other geophysical exploration data. In order to reflect the ground survey information in Japan to the Vp structure, we investigated the relationships among Vs, Vp, and depth by using PS-logging data at all K-NET and KiK-net sites. Vp values are concentrated at around 500 m/s and 1,500 m/s when Vs is lower than 1,000 m/s, where these concentrated areas show two distinctive characteristics of unsaturated and saturated soil, respectively. Many Vp values in the layer shallower than 4 m are around 500 m/s, which suggests the dominance of unsaturated soil, while many Vp values in the layer deeper than 4 m are larger than 1,500 m/s, which suggests the dominance of saturated soil there. We also investigated those relationships for different soil types at K-NET sites. Although each soil type has its own depth range, all soil types show similar relationships among Vs, Vp, and depth. Then, considering the depth profile of Vp, we divided the dataset into two by the depth, which is shallower or deeper than 4 m, and calculated the geometrical mean of Vp and the geometrical standard deviation in every Vs bins of 200 m/s. Finally, we obtained the regression curves for the average and standard deviation of Vp estimated from Vs to get the Vp conversion functions from Vs, which can be applied to a wide Vs range. We also obtained the regression curves for two datasets with Vp lower and higher than 1,200 m/s. These regression curves can be applied when the groundwater level is known. In addition, we obtained the regression curves for density from Vs or Vp. An example of the application for those relationships in the velocity inversion is shown.


2021 ◽  
Vol 9 (1) ◽  
Author(s):  
Britta Wawerzinek ◽  
Hermann Buness ◽  
Hartwig von Hartmann ◽  
David C. Tanner

AbstractThere are many successful geothermal projects that exploit the Upper Jurassic aquifer at 2–3 km depth in the German Molasse Basin. However, up to now, only P-wave seismic exploration has been carried out. In an experiment in the Greater Munich area, we recorded S-waves that were generated by the conventional P-wave seismic survey, using 3C receivers. From this, we built a 3D volume of P- to S-converted (PS) waves using the asymptotic conversion point approach. By combining the P-volume and the resulting PS-seismic volume, we were able to derive the spatial distribution of the vp/vs ratio of both the Molasse overburden and the Upper Jurassic reservoir. We found that the vp/vs ratios for the Molasse units range from 2.0 to 2.3 with a median of 2.15, which is much higher than previously assumed. This raises the depth of hypocenters of induced earthquakes in surrounding geothermal wells. The vp/vs ratios found in the Upper Jurassic vary laterally between 1.5 and 2.2. Since no boreholes are available for verification, we test our results against an independently derived facies classification of the conventional 3D seismic volume and found it correlates well. Furthermore, we see that low vp/vs ratios correlate with high vp and vs velocities. We interpret the latter as dolomitized rocks, which are connected with enhanced permeability in the reservoir. We conclude that 3C registration of conventional P-wave surveys is worthwhile.


Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1838-1842 ◽  
Author(s):  
C. M. Schmeissner ◽  
K. T. Spikes ◽  
D. W. Steeples

Ultrashallow seismic reflection surveys require dense spatial sampling during data acquisition, which increases their cost. In previous efforts to find ways to reduce these costs, we connected geophones rigidly to pieces of channel iron attached to a farm implement. This method allowed us to plant the geophones in the ground quickly and automatically. The rigidly interconnected geophones used in these earlier studies detected first‐arrival energy along with minor interfering seismic modes, but they did not detect seismic reflections. To examine further the feasibility of developing rigid geophone emplacement systems to detect seismic reflections, we experimented with four pieces of channel iron, each 2.7 m long and 10 cm wide. Each segment was equipped with 18 geophones rigidly attached to the channel iron at 15‐cm intervals, and the spikes attached to all 18 geophones were pushed into the ground simultaneously. The geophones detected both refracted and reflected energy; however, no significant signal distortion or interference attributable to the rigid coupling of the geophones to the channel iron was observed in the data. The interfering seismic modes mentioned from the previous experiments were not detected, nor was any P‐wave propagation noted within the channel iron. These results show promise for automating and reducing the cost of ultrashallow seismic reflection and refraction surveys.


Author(s):  
Suresh Dande ◽  
◽  
Robert R. Stewart ◽  
Nikolay Dyaur ◽  
◽  
...  

Laboratory physical models play an important role in understanding rock properties and wave propagation, both theoretically and at the field scale. In some cases, 3D-printing technology can be adopted to construct complex rock models faster, more inexpensively, and with more specific features than previous model-building techniques. In this study, we use 3D-printed rock models to assist in understanding the effects of various fluids (air, water, engine oil, crude oil, and glycerol) on the models’ elastic properties. We first used a 3D-printed, 1-in. cube-shaped layered model. This model was created with a 6% primary porosity and a bulk density of 0.98 g/cc with VTI anisotropy. We next employed a similar cube but with horizontal inclusions embedded in the layered background, which contributed to its total 24% porosity (including primary porosity). For air to liquid saturation, P-velocities increased for all liquids in both models, with the highest increase being with glycerol (57%) and an approximately 45% increase for other fluids in the inclusion model. For the inclusion model (dry and saturated), we observed a greater difference between two orthogonally polarized S-wave velocities (Vs1 and Vs2) than between two P-wave velocities (VP0 and VP90). We attribute this to the S2-wave (polarized normal to both the layering and the plane of horizontal inclusions), which appears more sensitive to horizontal inclusions than the P-wave. For the inclusion model, Thomsen’s P-wave anisotropic parameter (ɛ) decreased from 26% for the air case to 4% for the water-saturated cube and to 1% for glycerol saturation. The small difference between the bulk modulus of the frame and the pore fluid significantly reduces the velocity anisotropy of the medium, making it almost isotropic. We compared our experimental results with theory and found that predictions using Schoenberg’s linear slip theory combined with Gassmann’s anisotropic equation were closer to actual measurements than Hudson’s isotropic calculations. This work provides insights into the usefulness of 3D-printed models to understand elastic rock properties and wave propagation under various fluid saturations.


2014 ◽  
Vol 988 ◽  
pp. 502-507 ◽  
Author(s):  
Shao Bo Chai ◽  
Jian Chun Li ◽  
Hai Bo Li ◽  
Ya Qun Liu

According to the displacement discontinuity method and the conservation of momentum at the wave fronts, analysis for cylindrical P-wave propagation across a linear elastic rock joint is carried out. Considering the energy variation for wave propagation in one medium, the wave propagation equation was derived and expressed in an iterative form. The transmission and reflection coefficients are then obtained from the equation. By verification, the results agree very well with those from the existing results.


Geophysics ◽  
2005 ◽  
Vol 70 (4) ◽  
pp. R45-R56 ◽  
Author(s):  
Lars Nielsen ◽  
Hans Thybo ◽  
Martin Glendrup

Seismic wide-angle data were recorded to more than 300-km offset from powerful airgun sources during the MONA LISA experiments in 1993 and 1995 to determine the seismic-velocity structure of the crust and uppermost mantle along three lines in the southeastern North Sea with a total length of 850 km. We use the first arrivals observed out to an offset of 90 km to obtain high-resolution models of the velocity structure of the sedimentary layers and the upper part of the crystalline crust. Seismic tomographic traveltime inversion reveals 2–8-km-thick Paleozoic sedimentary sequences with P-wave velocities of 4.5–5.2 km/s. These sedimentary rocks are situated below a Mesozoic-Cenozoic sequence with variable thickness: ∼2–3 km on the basement highs, ∼2–4 km in the Horn Graben and the North German Basin, and ∼6–7 km in the Central Graben. The thicknesses of the Paleozoic sedimentary sequences are ∼3–5 km in the Central Graben, more than 4 km in the Horn Graben, up to ∼4 km on the basement highs, and up to 8 km in the North German Basin. The Paleozoic strata are clearly separated from the shallower and younger sequences with velocities of ∼1.8–3.8 km/s and the deeper crystalline crust with velocities of more than 5.8–6.0 km/s in the tomographic P-wave velocity model. Resolution tests show that the existence of the Paleozoic sediments is well constrained by the data. Hence, our wide-angle seismic models document the presence of Paleozoic sediments throughout the southeastern North Sea, both in the graben structures and in deep basins on the basement highs.


Sign in / Sign up

Export Citation Format

Share Document