• Geophysics and Geophysical Exploration
• /
• v.16 no.4
• /
• pp.240-249
• /
• 2013

The seismic data quality of marine geological and engineering survey deteriorates because of the sea swell. We often conduct a marine survey when the swell height is about 1 ~ 2 m. The swell effect correction is required to enhance the horizontal continuity of seismic data and satisfy the resolution less than 1 m. We applied the swell correction to the 8 channel high-resolution airgun seismic data and 3.5 kHz subbottom profiler (SBP) data. The correct sea bottom detection is important for the swell correction. To detect the sea bottom, we used maximum amplitude of seismic signal around the expected sea bottom, and picked the first increasing point larger than threshold value related with the maximum amplitude. To find sea bottom easily in the case of the low quality data, we transformed the input data to envelope data or the cross-correlated data using the sea bottom wavelet. We averaged the picked sea bottom depths and calculated the correction values. The maximum correction of the airgun data was about 0.8 m and the maximum correction of two kinds of 3.5 kHz SBP data was 0.5 m and 2.0 m respectively. We enhanced the continuity of the subsurface layer and produced the high quality seismic section using the proper methods of swell correction.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.50-59
• /
• 2006

The risk of fault reactivation in the Gippsland Basin was calculated using the FAST (Fault Analysis Seal Technology) technique, which determines fault reactivation risk by estimating the increase in pore pressure required to cause reactivation within the present-day stress field. The stress regime in the Gippsland Basin is on the boundary between strike-slip and reverse faulting: maximum horizontal stress $({\sim}\;40.5\;Mpa/km)$ > vertical stress (21 Mpa/km) ${\sim}$ minimum horizontal stress (20 MPa/km). Pore pressure is hydrostatic above the Campanian Volcanics of the Golden Beach Subgroup. The NW-SE maximum horizontal stress orientation $(139^{\circ}N)$ determined herein is broadly consistent with previous estimates, and verifies a NW-SE maximum horizontal stress orientation in the Gippsland Basin. Fault reactivation risk in the Gippsland Basin was calculated using two fault strength scenarios; cohesionless faults $(C=0;{\mu}=0.65)$ and healed faults $(C=5.4;\;{\mu}=0.78)$. The orientations of faults with relatively high and relatively low reactivation potential are almost identical for healed and cohesionless fault strength scenarios. High-angle faults striking NE-SW are unlikely to reactivate in the current stress regime. High-angle faults oriented SSE-NNW and ENE-WSW have the highest fault reactivation risk. Additionally, low-angle faults (thrust faults) striking NE-SW have a relatively high risk of reactivation. The highest reactivation risk for optimally oriented faults corresponds to an estimated pore pressure increase (Delta-P) of 3.8 MPa $({\sim}548\;psi)$ for cohesionless faults and 15.6 MPa $({\sim}2262\;psi)$ for healed faults. The absolute values of pore pressure increase obtained from fault reactivation analysis presented in this paper are subject to large errors because of uncertainties in the geomechanical model (in situ stress and rock strength data). In particular, the maximum horizontal stress magnitude and fault strength data are poorly constrained. Therefore, fault reactivation analysis cannot be used to directly measure the maximum allowable pore pressure increase within a reservoir. We argue that fault reactivation analysis of this type can only be used for assessing the relative risk of fault reactivation and not to determine the maximum allowable pore pressure increase a fault can withstand prior to reactivation.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.60-66
• /
• 2006

Three differing sandstones, two synthetic and one field sample, have been tested ultrasonically under a range of confining pressures and pore pressures representative of in-situ reservoir pressures. These sandstones include: a synthetic sandstone with calcite intergranular cement produced using the CSIRO Calcite In-situ Precipitation Process (CIPS); a synthetic sandstone with silica intergranular cement; and a core sample from the Otway Basin Waarre Formation, Boggy Creek 1 well, from the target lithology for a trial $CO_2$ pilot project. Initial testing was carried on the cores at "room-dried" conditions, with confining pressures up to 65 MPa in steps of 5 MPa. All cores were then flooded with $CO_2$, initially in the gas phase at 6 MPa, $22^{\circ}C$, then with liquid-phase $CO_2$ at a temperature of $22^{\circ}C$ and pressures from 7 MPa to 17 MPa in steps of 5 MPa. Confining pressures varied from 10 MPa to 65 MPa. Ultrasonic waveforms for both P- and S-waves were recorded at each effective pressure increment. Velocity versus effective pressure responses were calculated from the experimental data for both P- and S-waves. Attenuations $(1/Q_p)$ were calculated from the waveform data using spectral ratio methods. Theoretical calculations of velocity as a function of effective pressure for each sandstone were made using the $CO_2$ pressure-density and $CO_2$ bulk modulus-pressure phase diagrams and Gassmann effective medium theory. Flooding the cores with gaseous phase $CO_2$ produced negligible change in velocity-effective stress relationships compared to the dry state (air saturated). Flooding with liquid-phase $CO_2$ at various pore pressures lowered velocities by approximately 8% on average compared to the air-saturated state. Attenuations increased with liquid-phase $CO_2$ flooding compared to the air-saturated case. Experimental data agreed with the Gassmann calculations at high effective pressures. The "critical" effective pressure, at which agreement with theory occurred, varied with sandstone type. Discrepancies are thought to be due to differing micro-crack populations in the microstructure of each sandstone type. The agreement with theory at high effective pressures is significant and gives some confidence in predicting seismic behaviour under field conditions when $CO_2$ is injected.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.129-138
• /
• 2006

Helicopter-borne aeromagnetic surveys at two different times separated by three years were conducted to better understand the shallow subsurface structure of the Vulcano and Lipari volcanic complex, Aeolian Islands, southern Italy, and also to monitor the volcanic activity of the area. As there was no meaningful difference between the two magnetic datasets to imply an apparent change of the volcanic activity, the datasets were merged to produce an aeromagnetic map with wider coverage than was given by a single dataset. Apparent magnetisation intensity mapping was applied to terrain-corrected magnetic anomalies, and showed local magnetisation highs in and around Fossa Cone, suggesting heterogeneity of the cone. Magnetic modelling was conducted for three of those magnetisation highs. Each model implied the presence of concealed volcanic products overlain by pyroclastic rocks from the Fossa crater. The model for the Fossa crater area suggests a buried trachytic lava flow on the southern edge of the present crater. The magnetic model at Forgia Vecchia suggests that phreatic cones can be interpreted as resulting from a concealed eruptive centre, with thick latitic lavas that fill up Fossa Caldera. However, the distribution of lavas seems to be limited to a smaller area than was expected from drilling results. This can be explained partly by alteration of the lavas by intense hydrothermal activity, as seen at geothermal areas close to Porto Levante. The magnetic model at the north-eastern Fossa Cone implies that thick lavas accumulated as another eruption centre in the early stage of the activity of Fossa. Recent geoelectric surveys showed high-resistivity zones in the areas of the last two magnetic models.

• Geophysics and Geophysical Exploration
• /
• v.9 no.1
• /
• pp.30-36
• /
• 2006

Japan's first pilot-scale $CO_2$ sequestration experiment has been conducted in Nagaoka, where 10400 t of $CO_2$ have been injected in an onshore aquifer at a depth of about 1100 m. Among various measurements conducted at the site for monitoring the injected $CO_2$, we conducted time-lapse crosswell seismic tomography between two observation wells to determine the distribution of $CO_2$ in the aquifer by the change of P-wave velocities. This paper reports the results of the crosswell seismic tomography conducted at the site. The crosswell seismic tomography measurements were carried out three times; once before the injection as a baseline survey, and twice during the injection as monitoring surveys. The velocity tomograms resulting from the monitoring surveys were compared to the baseline survey tomogram, and velocity difference tomograms were generated. The velocity difference tomograms showed that velocity had decreased in a part of the aquifer around the injection well, where the injected $CO_2$ was supposed to be distributed. We also found that the area in which velocity had decreased was expanding in the formation up-dip direction, as increasing amounts of $CO_2$ were injected. The maximum velocity reductions observed were 3.0% after 3200 t of $CO_2$ had been injected, and 3.5% after injection of 6200 t of $CO_2$. Although seismic tomography could map the area of velocity decrease due to $CO_2$ injection, we observed some contradictions with the results of time-lapse sonic logging, and with the geological condition of the cap rock. To investigate these contradictions, we conducted numerical experiments simulating the test site. As a result, we found that part of the velocity distribution displayed in the tomograms was affected by artefacts or ghosts caused by the source-receiver geometry for the crosswell tomography in this particular site. The maximum velocity decrease obtained by tomography (3.5%) was much smaller than that observed by sonic logging (more than 20%). The numerical experiment results showed that only 5.5% velocity reduction might be observed, although the model was given a 20% velocity reduction zone. Judging from this result, the actual velocity reduction can be more than 3.5%, the value we obtained from the field data reconstruction. Further studies are needed to obtain more accurate velocity values that are comparable to those obtained by sonic logging.

• Geophysics and Geophysical Exploration
• /
• v.17 no.1
• /
• pp.1-10
• /
• 2014

To reveal and classify site characteristics in densely populated areas in Chuncheon, Korea, Rayleigh-waves were recorded at 50 sites including four sites in the forest area using four 1-Hz velocity sensors and 24 4.5-Hz vertical geophones during the period of January 2011 to May 2013. Dispersion curves of the Rayleigh waves obtained by the extended spatial autocorrelation method were inverted to derive shear-wave velocity ($v_s$) models comprising 40 horizontal layers of 1-m thickness. Depths to weathered rocks ($D_b$), shear wave velocities of these basement rocks ($v_s^b$), average velocities of the overburden layer ($\bar{v}_s^s$), and the average velocity to a depth of 30 m ($v_s30$), were then derived from those models. The estimated values of $D_b$, $v_s^b$, $\bar{v}_s^s$, and $v_s30$ for 46 sites at lower altitudes were in the ranges of 5 to 29 m, 404 to 561 m/s, 208 to 375 ms/s, and 226 to 583 m/s, respectively. According to the Korean building code for seismic design, the estimated $v_s30$ indicates that the lower altitude areas in Chuncheon are classified as $S_C$ (very dense soil and soft rock) or $S_D$ (stiff soil). To determine adequate proxies for $v_s30$, we compared the computed values with land cover, lithology, topographic slope, and surface elevation at each of the measurement sites. Due to a weak correlation (r = 0.41) between $v_s30$ and elevation, the best proxy of them, applications of this proxy to Chuncheon of a relatively small area seem to be limited.

• Geophysics and Geophysical Exploration
• /
• v.8 no.2
• /
• pp.145-155
• /
• 2005

A two-dimensional (2-D) interpretation of MT data has been performed for the purpose of fracture detection for geothermal development. Remote stations have been operated in Kyushu, Japan (480 km apart) as well as in Korea (60 km and 165 km apart in 2002 and 2003 data set, respectively). Apparent resistivity and phase curves calculated by remote processing with the Japan remote data showed enough quality for 2-D inversion for the whole frequency range. Remote reference processing with Korea remote reference data also showed quite good continuity in apparent resistivity and phase curves except some noisy frequency bands; around the power frequency, 60 Hz, and around the dead band $10^{-1}Hz\;Hz\;\~1\;Hz$, where the natural EM signal is known to be very weak. Even though the subsurface showed severe three-dimensional (3-D) characteristics in the survey area so that 2-D inversion by itself could not give enough information for deep geological structures, the 2-D inversion for the 5 survey lines showed several common features. The conductive semi-consolidate mudstone layer is dipping from north to south (about 500 m depth on the south and 200 m on the north most part of the survey area). The boundary between the low (L-2) and high (H-2) resistivity anomalies can be thought as a major fault with strike $N15^{\circ}E$, passing through the sites 206, 112 and 414. The shallow (< 1 km) conductive anomalies (L-4) seem to be fracture zones having strike E-W (at site 105) and $N60^{\circ}W$ (at site 434). And there exists a conductive layer in the western and west-southern part of the survey area in the depth below $2\~3\;km$, for which further investigation is to be needed.

• Geophysics and Geophysical Exploration
• /
• v.7 no.3
• /
• pp.164-173
• /
• 2004

To delineate geothermal water movement at the Pohang geothermal development site, Self-Potential (SP) survey and monitoring were carried out during pumping tests. Before drilling, background SP data have been gathered to figure out overall potential distribution of the site. The pumping test was performed in two separate periods: 24 hours in December 2003 and 72 hours in March 2004. SP monitoring started several days before the pumping tests with a 128-channel automatic recording system. The background SP survey showed a clear positive anomaly at the northern part of the boreholes, which may be interpreted as an up-flow Bone of the deep geothermal water due to electrokinetic potential generated by hydrothermal circulation. The first and second SP monitoring during the pumping tests performed to figure out the fluid flow in the geothermal reservoir but it was not easy to see clear variations of SP due to pumping and pumping stop. Since the area is covered by some 360 m-thick tertiary sediments with very low electrical resistivity (less than 10 ohm-m), the electrokinetic potential due to deep groundwater flow resulted in being seriously attenuated on the surface. However, when we compared the variation of SP with that of groundwater level and temperature of pumping water, we could identify some areas responsible to the pumping. Dominant SP changes are observed in the south-west part of the boreholes during both the preliminary and long-term pumping periods, where 3-D magnetotelluric survey showed low-resistivity anomaly at the depth of $600m\~1,000m$. Overall analysis suggests that there exist hydraulic connection through the southwestern part to the pumping well.

• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.29-36
• /
• 2016

Recently, the studies about rock physics model (RPM) in shale reservoir are widely performed. In shale reservoir, the degree of the maturity can be estimated by kerogen and GOR (Gas-Oil Ratio). The researches on the rock physics model of shale reservoir with the amount of kerogen have been actively carried out but not with GOR. Thus, in this study, we analyzed the changes in seismic velocity and density, and AVO (Amplitude Variation with Offset) response depending on changes in GOR and the amount of kerogen. Since the shale consists of plate-like particles, it has vertical transverse isotropy (VTI). Therefore we estimated the seismic velocity and density by using Backus averaging method and analyzed AVO responses based on these estimated properties. The results of analysis showed that the changes in the velocity with the GOR variation are small but the velocity changes with the variation in kerogen amount are relatively larger. In case, GOR 180 (Litre/Litre) which is boundary between heavy oil and light oil, when volume fraction of kerogen increased from 5% to 35%, the P-wave velocity normal to the layering increased 51%. That is, it helps estimating maturity of kerogen through the velocity. Meanwhile, when rates of oil-gas mixture are large, the effect of GOR variation on the velocity change became larger. In case volume fraction of kerogen is 5%, the P-wave velocity normal to the layering was estimated $1.46km/s^2$ in heavy oil (GOR 40) but $1.36km/s^2$ in light oil (GOR 300). The AVO responses analysis showed class 4 regardless of the GOR and amount of kerogen because variation of poisson's ratio is small. Therefore, shale reservoir has possibility to have class 4.

• Geophysics and Geophysical Exploration
• /
• v.19 no.1
• /
• pp.1-10
• /
• 2016

The hypocenter distribution of microseismic events generated by hydraulic fracturing for shale gas development provides essential information for understanding characteristics of fracture network. In this study, we evaluate how inaccurate velocity models influence the inversion results of two widely used location programs, hypoellipse and hypoDD, which are developed based on an iterative linear inversion. We assume that 98 stations are densely located inside the circle with a radius of 4 km and 5 artificial hypocenter sets (S0 ~ S4) are located from the center of the network to the south with 1 km interval. Each hypocenter set contains 25 events placed on the plane. To quantify accuracies of the inversion results, we defined 6 parameters: difference between average hypocenters of assumed and inverted locations, $d_1$; ratio of assumed and inverted areas estimated by hypocenters, r; difference between dip of the reference plane and the best fitting plane for determined hypocenters, ${\theta}$; difference between strike of the reference plane and the best fitting plane for determined hypocenters, ${\phi}$; root-mean-square distance between hypocenters and the best fitting plane, $d_2$; root-mean-square error in horizontal direction on the best fitting plane, $d_3$. Synthetic travel times are calculated for the reference model having 1D layered structure and the inaccurate velocity model for the inversion is constructed by using normal distribution with standard deviations of 0.1, 0.2, and 0.3 km/s, respectively, with respect to the reference model. The parameters $d_1$, r, ${\theta}$, and $d_2$ show positive correlation with the level of velocity perturbations, but the others are not sensitive to the perturbations except S4, which is located at the outer boundary of the network. In cases of S0, S1, S2, and S3, hypoellipse and hypoDD provide similar results for $d_1$. However, for other parameters, hypoDD shows much better results and errors of locations can be reduced by about several meters regardless of the level of perturbations. In light of the purpose to understand the characteristics of hydraulic fracturing, $1{\sigma}$ error of velocity structure should be under 0.2 km/s in hypoellipse and 0.3 km/s in hypoDD.