Marcellus fracture characterization using P-wave azimuthal velocity attributes: Comparison with production and outcrop data

2015 ◽  
Vol 3 (3) ◽  
pp. SU1-SU15 ◽  
Author(s):  
Tanya L. Inks ◽  
Terry Engelder ◽  
Edward Jenner ◽  
Bruce Golob ◽  
Jacki S. Hocum ◽  
...  

Analysis of two 3D surveys, available well data, published outcrop data and subsurface information, as well as production data available from the state of Pennsylvania, demonstrates that wide-azimuth seismic is sensitive to variations in fracturing at the scale of individual pads or even individual wells. These variations in fracturing begin to explain why production varies significantly, even locally, within the Marcellus Shale gas play. Rose diagrams from quantitative fracture analysis using azimuthal seismic velocity volumes were compared with published data from Appalachian black shale outcrops and subsurface fracture models proposed in various papers to validate the results from subsurface data. These analyses provided insight into the rock fabric and the presence of systematic joints that likely affect production. There was a strong correlation between the low anisotropy and low heterogeneity of anisotropy and high estimated ultimate recovery (EUR). Additionally, interpreted fracture trend azimuths differed between areas of larger gas EUR and areas of smaller gas EUR as defined by decline curve analysis. Some perforations were likely to perform much better than others along the borehole, based on observed heterogeneity in the seismic profiles and map view.

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.


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.


2009 ◽  
Author(s):  
Zhirang Zhang ◽  
Mingjin Zhao ◽  
Hongluo Wang ◽  
Shaoguo Yang

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.


2021 ◽  
Vol 9 ◽  
Author(s):  
Hidayat Hidayat ◽  
Andri Dian Nugraha ◽  
Awali Priyono ◽  
Marjiyono Marjiyono ◽  
Januar H. Setiawan ◽  
...  

The Banyumas Basin is a tertiary sedimentary basin located in southern Central Java, Indonesia. Due to the presence of volcanic deposits, 2-D seismic reflection methods cannot provide a good estimation of the sediment thickness and the subsurface geology structure in this area. In this study, the passive seismic tomography (PST) method was applied to image the 3-D subsurface Vp, Vs, and Vp/Vs ratio. We used 70 seismograph borehole stations with a recording duration of 177 days. A total of 354 events with 9, 370 P and 9, 368 S phases were used as input for tomographic inversion. The checkshot data of a 4, 400-meter deep exploration well (Jati-1) located within the seismic network were used to constrain the shallow crustal layer of the initial 1-D velocity model. The model resolution of the tomographic inversions was assessed using the checkerboard resolution test (CRT), the diagonal resolution element (DRE), and the derivative weight sum (DWS). Using the obtained Vp, Vs, and Vp/Vs ratio, we were able to sharpen details of the geological structures within the basin from previous geological studies, and a fault could be well-imaged at a depth of 4 km. We interpreted this as the main dextral strike-slip fault that controls the pull apart process of the Banyumas Basin. The thickness of the sediment layers, as well as its layering, were also could be well determined. We found prominent features of the velocity contrast that aligned very well with the boundary between the Halang and Rambatan formations as observed in the Jati-1 well data. Furthermore, an anticline structure, which is a potential structural trap for the petroleum system in the Banyumas Basin, was also well imaged. This was made possible due to the dense borehole seismographic stations which were deployed in the study area.


Geophysics ◽  
2018 ◽  
Vol 83 (2) ◽  
pp. KS1-KS10 ◽  
Author(s):  
Zhishuai Zhang ◽  
James W. Rector ◽  
Michael J. Nava

We have studied microseismic data acquired from a geophone array deployed in the horizontal section of a well drilled in the Marcellus Shale near Susquehanna County, Pennsylvania. Head waves were used to improve event location accuracy as a substitution for the traditional P-wave polarization method. We identified that resonances due to poor geophone-to-borehole coupling hinder arrival-time picking and contaminate the microseismic data spectrum. The traditional method had substantially greater uncertainty in our data due to the large uncertainty in P-wave polarization direction estimation. We also identified the existence of prominent head waves in some of the data. These head waves are refractions from the interface between the Marcellus Shale and the underlying Onondaga Formation. The source location accuracy of the microseismic events can be significantly improved by using the P-, S-wave direct arrival times and the head wave arrival times. Based on the improvement, we have developed a new acquisition geometry and strategy that uses head waves to improve event location accuracy and reduce acquisition cost in situations such as the one encountered in our study.


2018 ◽  
Vol 18 (11) ◽  
pp. 06018030
Author(s):  
Xiao-Ping Zhou ◽  
Xiao-Cheng Huang ◽  
Jin-Xin Li

Sign in / Sign up

Export Citation Format

Share Document