This study was conducted to develop a reservoir modelling workflow to reproduce the heterogeneous distribution of effective permeability that impacts on the performance of SAGD (Steam Assisted Gravity Drainage), the in-situ bitumen recovery technique in the Athabasca Oil Sands. Lithologic facies distribution is the main cause of the heterogeneity in bitumen reservoirs in the study area. The target formation consists of sand with mudstone facies in a fluvial-to-estuary channel system, where the mudstone interrupts fluid flow and reduces effective permeability. In this study, the lithologic facies is classified into three classes having different characteristics of effective permeability, depending on the shapes of mudstones. The reservoir modelling workflow of this study consists of two main modules; facies modelling and permeability modelling. The facies modelling provides an identification of the three lithologic facies, using a stochastic approach, which mainly control the effective permeability. The permeability modelling populates mudstone volume fraction first, then transforms it into effective permeability. A series of flow simulations applied to mini-models of the lithologic facies obtains the transformation functions of the mudstone volume fraction into the effective permeability. Seismic data contribute to the facies modelling via providing prior probability of facies, which is incorporated in the facies models by geostatistical techniques. In particular, this study employs a probabilistic neural network utilising multiple seismic attributes in facies prediction that improves the prior probability of facies. The result of using the improved prior probability in facies modelling is compared to the conventional method using a single seismic attribute to demonstrate the improvement in the facies discrimination. Using P-wave velocity in combination with density in the multiple seismic attributes is the essence of the improved facies discrimination. This paper also discusses sand matrix porosity that makes P-wave velocity differ between the different facies in the study area, where the sand matrix porosity is uniquely evaluated using log-derived porosity, P-wave velocity and photographically-predicted mudstone volume.
Kim, Seyoon;Kim, Jungtaek;Lee, Minhee;Wang, Sookyun
Economic and Environmental Geology
/
v.49
no.6
/
pp.469-477
/
2016
To evaluate the $CO_2$ storage capacity for the reservoir rock, the laboratory scale technique to measure the amount of $scCO_2$, replacing pore water of the reservior rock after the $CO_2$ injection was developed in this study. Laboratory experiments were performed to measure the $scCO_2$ displacement efficiency of the conglomerate and the sandstone in Janggi basin, which are classified as available $CO_2$ storage rocks in Korea. The high pressurized stainless steel cell containing two different walls was designed and undisturbed rock cores acquired from the deep drilling site around Janggi basin were used for the experiments. From the lab experiments, the average $scCO_2$ displacement efficiency of the conglomerate and the sandstone in Janggi basin was measured at 31.2% and 14.4%, respectively, which can be used to evaluate the feasibility of the Janggi basin as a $scCO_2$ storage site in Korea. Assuming that the effective radius of the $CO_2$ storage formations is 250 m and the average thickness of the conglomerate and the sandstone formation under 800 m in depth is 50 m each (from data of the drilling profile and the geophysical survey), the $scCO_2$ storage capacity of the reservoir rocks around the probable $scCO_2$ injection site in Janggi basin was calculated at 264,592 metric ton, demonstrating that the conglomerate and the sandstone formations in Janggi basin have a great potential for use as a pilot scale test site for the $CO_2$ storage in Korea.
Ruth, Peter J. van;Nelson, Emma J.;Hillis, Richard R.
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.
High-resolution shallow marine seismic surveys have been carried out for the resources exploration, engineering applications and Quaternary mapping. To improve the resolution of subsurface structure image, multichannel digital technique has been applied. The quality of the image depends on the vertical and horizontal resolution and signal to noise (S/N) ratio which are associated with the data acquisition parameters such as sample interval, common midpoint (CMP) interval and CMP fold. To understand the effect of the acquisition parameters, a test survey was carried out off Yeosu and the acquired data were analyzed. A 30 $in^3$ small air gun was used as a seismic source and 8 channel streamer cable with a 5 m group interval was used as a receiver. The data were digitally recorded with a shot interval of 2 s and sample interval of 0.1 ms. The acquired data were resampled with various sample intervals, CMP intervals and CMP folds. The resampled data were processed, plotted as seismic sections and compared each other. The analysis results show that thin bed structure with ${\~}1m$ thickness and ${\~}6^{\circ}$ slope can be imaged with good resolution and continuity and low noise using the acquisition parameters with a sample interval shorter than 0.2 ms, CMP interval shorter than 2.5 m and CMP fold more than 4. Because seismic resolution is associated with the acquisition parameters, the quality of the subsurface structure can be imaged successfully using suitable and optimum acquisition parameters.
To identify the potential area of gas hydrate in the Ulleung Basin, 2-D and 3-D seismic surveys using R/V Tamhae II were conducted in 2005 and 2006. Seismic survey equipment consisted of navigation system, recording system, streamer cable and air-gun source. For reliable velocity analysis in a deep sea area where water depths are mostly greater than 1,000 m and the target depth is up to about 500 msec interval below the seafloor, 3-km-long streamer and 1,035 $in^3$ tuned air-gun array were used. During the survey, a suite of quality control operations including source signature analysis, 2-D brute stack, RMS noise analysis and FK analysis were performed. The source signature was calculated to verify its conformity to quality specification and the gun dropout test was carried out to examine signature changes due to a single air gun's failure. From the online quality analysis, we could conclude that the overall data quality was very good even though some seismic data were affected by swell noise, parity error, spike noise and current rip noise. Especially, by checking the result of data quality enhancement using FK filtering and missing trace restoration technique for the 3-D seismic data inevitably contaminated with current rip noises, the acquired data were accepted and the field survey could be conducted continuously. Even in survey areas where the acquired data would be unsuitable for quality specification, the marine seismic survey efficiency could be improved by showing the possibility of noise suppression through onboard data processing.
Bang, Minkyu;Oh, Seokmin;Seol, Soon Jee;Lee, Ki Ha;Cho, Seong-Jun
Geophysics and Geophysical Exploration
/
v.21
no.3
/
pp.150-161
/
2018
Recently, unmanned aircraft EM (electromagnetic) survey based on ICT (Information and Communication Technology) has been widely utilized because of the efficiency in regional survey. We performed the theoretical study on the unmanned airship EM system developed by KIGAM (Korea Institute of Geoscience and Mineral resources) as part of the practical application of unmanned aircraft EM survey. Since this system has different configurations of transmitting and receiving loops compared to the conventional aircraft EM systems, a new technique is required for the appropriate interpretation of measured responses. Therefore, we proposed a method to calculate the EM field for the arbitrary shaped transmitter and verified its validity through the comparison with analytic solution for circular loop. In addition, to simulate the magnetic responses by three-dimensionally (3D) distributed anomalies, we have adapted our algorithm to 3D frequency-domain EM modeling algorithm based on the edge-FEM (finite element method). Though the analysis on magnetic field responses from a subsurface anomaly, it was found that the response decreases as the depth of the anomaly increases or the flight altitude increases. Also, it was confirmed that the response became smaller as the resistivity of the anomaly increases. However, a nonlinear trend of the out-of-phase component is shown depending on the depth of the anomaly and the used frequency, that makes it difficult to apply simple analysis based on the mapping of the magnitude of the responses and can cause the non-uniqueness problem in calculating the apparent resistivity. Thus, it is a prerequisite to analyze the appropriate frequency band and flight altitude considering the purpose of the survey and the site conditions when conducting a survey using the unmanned aircraft EM system.
Park, Chan Hong;Kim, Jeong U;Heo, Sik;Won, Jung Seon;Seok, Bong Chul;Yu, Hae Su
Journal of the Korean Geophysical Society
/
v.2
no.1
/
pp.27-38
/
1999
The errors between track segments or at the cross-over points of shipborne gravity were successfully reduced by applying a cross-over error adjustment technique using satellite gravity. The integration of shipborne and satellite altimeter-implied free-air gravity anomalies after the cross-over error adjustment resulted in a high resolution gravity map which contains both short and long wavelength components. The successful adjustment of the cross-over errors in the shipborne gravity using the satellite gravity suggests that the shipborne gravity can be combined with the satellite anomalies characterized by a stable and long wavelength component. The resulting free-air anomaly map is evenly harmonized with both short and long wavelength anomalies. Thus the corrected anomaly map can be better used for the geological interpretation. Free-air anomalies with more than 140 mGal in total variations generally correspond to the seafloor topographic changes in their regional patterns. A series of gravity highs are aligned from the Korea Plateau to the Oki Island, which are interpreted to be caused by seamounts or volcanic topographies. The gravity minima along the western and southern shelf edge are associated not only with the local basement morphology and thick sediment fill at the continental margin, but also possibly with the crustal edge effect known for passive continental margins. Series of NE-trending linear anomalies are possibly caused by a swarm of volcanic intrusions followed the initial opening of the Ulleung Basin. The linear high anomalies in the Ulleung Plateau are terminated by the straightly NNW-trending anomalies with a sharp gradient in its western boundary which indicates a fault-line scarp. The opposite side adjoined with the fault-line scarp shows no correlation with the fault-line scarp in geometry indicating that the block might be horizontally slided from the north. A gravity high in contrast to the deepening in seafloor toward the northeastern central Ulleung Basin is probably responsible for the thin crust and shallow seated mantle. The gravity minima along the western and southern shelf edge are associated not only with the local basement morphology and thick sediment fill at the continental margin, but also possibly with the crustal edge effect known for passive continental margins. Series of NE-trending linear anomalies are possibly caused by a swarm of volcanic intrusions followed the initial opening of the Ulleung Basin. The linear high anomalies in the Ulleung Plateau are terminated by the straightly NNW-trending anomalies with a sharp gradient in its western boundary which indicates a fault-line scarp. The opposite side adjoined with the fault-line scarp shows no correlation with the fault-line scarp in geometry indicating that the block might be horizontally slided from the north. A gravity high in contrast to the deepening in seafloor toward the northeastern central Ulleung Basin is probably suggestive of a thin crust and shallow seated mantle.
Recently, stratigraphic reservoirs are getting more attention than structural reservoirs which have mostly developed. However, recognizing stratigraphic thin gas reservoirs in a stacked section is usually difficult because of tuning effects. Moreover, if the reflections from the brine-saturated region of a thin layer have the same polarity with those from the gas-saturated region, we could not easily identify the gas reservoir with conventional data processing technique. In this study, we introduced a way to delineate the gas-saturated region in a thin layer reservoir using a spectral decomposition method. First of all, amplitude spectrum with the variation of the frequency and the incident angle was investigated for the medium which represents property of Class 3, Class 1 or Class 4 AVO response. The results show that the maximum difference in the amplitude spectra between brine and gas-saturated thin layers occurs around the peak frequency independent of the incident angle and the type of AVO responses. In addition, the amplitude spectra of the gas-saturated zone are greater than those of brine-saturated one in Class 3 and Class 4 at the peak frequency while those of phenomenon occur oppositely in Class 1. Based on the results, we applied spectral decomposition method to the stacked section in order to distinguish the gas-saturated zone from the brine-saturated zone in a thin layer reservoir. To verify our new method, we constructed a thin-layer velocity model which contains both gas and brine-saturated zones which have the same reflection polarities. As a result, in the spectral decomposed sections near the peak frequency obtained by Wigner-Ville Distribution (WVD), we could identify the difference between reflections from gas- and brinesaturated region in the thin layer reservoir, which was hardly distinguishable in the stacked section.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 2004년 10월 1일]
이용약관
제 1 장 총칙
제 1 조 (목적)
이 이용약관은 KoreaScience 홈페이지(이하 “당 사이트”)에서 제공하는 인터넷 서비스(이하 '서비스')의 가입조건 및 이용에 관한 제반 사항과 기타 필요한 사항을 구체적으로 규정함을 목적으로 합니다.
제 2 조 (용어의 정의)
① "이용자"라 함은 당 사이트에 접속하여 이 약관에 따라 당 사이트가 제공하는 서비스를 받는 회원 및 비회원을
말합니다.
② "회원"이라 함은 서비스를 이용하기 위하여 당 사이트에 개인정보를 제공하여 아이디(ID)와 비밀번호를 부여
받은 자를 말합니다.
③ "회원 아이디(ID)"라 함은 회원의 식별 및 서비스 이용을 위하여 자신이 선정한 문자 및 숫자의 조합을
말합니다.
④ "비밀번호(패스워드)"라 함은 회원이 자신의 비밀보호를 위하여 선정한 문자 및 숫자의 조합을 말합니다.
제 3 조 (이용약관의 효력 및 변경)
① 이 약관은 당 사이트에 게시하거나 기타의 방법으로 회원에게 공지함으로써 효력이 발생합니다.
② 당 사이트는 이 약관을 개정할 경우에 적용일자 및 개정사유를 명시하여 현행 약관과 함께 당 사이트의
초기화면에 그 적용일자 7일 이전부터 적용일자 전일까지 공지합니다. 다만, 회원에게 불리하게 약관내용을
변경하는 경우에는 최소한 30일 이상의 사전 유예기간을 두고 공지합니다. 이 경우 당 사이트는 개정 전
내용과 개정 후 내용을 명확하게 비교하여 이용자가 알기 쉽도록 표시합니다.
제 4 조(약관 외 준칙)
① 이 약관은 당 사이트가 제공하는 서비스에 관한 이용안내와 함께 적용됩니다.
② 이 약관에 명시되지 아니한 사항은 관계법령의 규정이 적용됩니다.
제 2 장 이용계약의 체결
제 5 조 (이용계약의 성립 등)
① 이용계약은 이용고객이 당 사이트가 정한 약관에 「동의합니다」를 선택하고, 당 사이트가 정한
온라인신청양식을 작성하여 서비스 이용을 신청한 후, 당 사이트가 이를 승낙함으로써 성립합니다.
② 제1항의 승낙은 당 사이트가 제공하는 과학기술정보검색, 맞춤정보, 서지정보 등 다른 서비스의 이용승낙을
포함합니다.
제 6 조 (회원가입)
서비스를 이용하고자 하는 고객은 당 사이트에서 정한 회원가입양식에 개인정보를 기재하여 가입을 하여야 합니다.
제 7 조 (개인정보의 보호 및 사용)
당 사이트는 관계법령이 정하는 바에 따라 회원 등록정보를 포함한 회원의 개인정보를 보호하기 위해 노력합니다. 회원 개인정보의 보호 및 사용에 대해서는 관련법령 및 당 사이트의 개인정보 보호정책이 적용됩니다.
제 8 조 (이용 신청의 승낙과 제한)
① 당 사이트는 제6조의 규정에 의한 이용신청고객에 대하여 서비스 이용을 승낙합니다.
② 당 사이트는 아래사항에 해당하는 경우에 대해서 승낙하지 아니 합니다.
- 이용계약 신청서의 내용을 허위로 기재한 경우
- 기타 규정한 제반사항을 위반하며 신청하는 경우
제 9 조 (회원 ID 부여 및 변경 등)
① 당 사이트는 이용고객에 대하여 약관에 정하는 바에 따라 자신이 선정한 회원 ID를 부여합니다.
② 회원 ID는 원칙적으로 변경이 불가하며 부득이한 사유로 인하여 변경 하고자 하는 경우에는 해당 ID를
해지하고 재가입해야 합니다.
③ 기타 회원 개인정보 관리 및 변경 등에 관한 사항은 서비스별 안내에 정하는 바에 의합니다.
제 3 장 계약 당사자의 의무
제 10 조 (KISTI의 의무)
① 당 사이트는 이용고객이 희망한 서비스 제공 개시일에 특별한 사정이 없는 한 서비스를 이용할 수 있도록
하여야 합니다.
② 당 사이트는 개인정보 보호를 위해 보안시스템을 구축하며 개인정보 보호정책을 공시하고 준수합니다.
③ 당 사이트는 회원으로부터 제기되는 의견이나 불만이 정당하다고 객관적으로 인정될 경우에는 적절한 절차를
거쳐 즉시 처리하여야 합니다. 다만, 즉시 처리가 곤란한 경우는 회원에게 그 사유와 처리일정을 통보하여야
합니다.
제 11 조 (회원의 의무)
① 이용자는 회원가입 신청 또는 회원정보 변경 시 실명으로 모든 사항을 사실에 근거하여 작성하여야 하며,
허위 또는 타인의 정보를 등록할 경우 일체의 권리를 주장할 수 없습니다.
② 당 사이트가 관계법령 및 개인정보 보호정책에 의거하여 그 책임을 지는 경우를 제외하고 회원에게 부여된
ID의 비밀번호 관리소홀, 부정사용에 의하여 발생하는 모든 결과에 대한 책임은 회원에게 있습니다.
③ 회원은 당 사이트 및 제 3자의 지적 재산권을 침해해서는 안 됩니다.
제 4 장 서비스의 이용
제 12 조 (서비스 이용 시간)
① 서비스 이용은 당 사이트의 업무상 또는 기술상 특별한 지장이 없는 한 연중무휴, 1일 24시간 운영을
원칙으로 합니다. 단, 당 사이트는 시스템 정기점검, 증설 및 교체를 위해 당 사이트가 정한 날이나 시간에
서비스를 일시 중단할 수 있으며, 예정되어 있는 작업으로 인한 서비스 일시중단은 당 사이트 홈페이지를
통해 사전에 공지합니다.
② 당 사이트는 서비스를 특정범위로 분할하여 각 범위별로 이용가능시간을 별도로 지정할 수 있습니다. 다만
이 경우 그 내용을 공지합니다.
제 13 조 (홈페이지 저작권)
① NDSL에서 제공하는 모든 저작물의 저작권은 원저작자에게 있으며, KISTI는 복제/배포/전송권을 확보하고
있습니다.
② NDSL에서 제공하는 콘텐츠를 상업적 및 기타 영리목적으로 복제/배포/전송할 경우 사전에 KISTI의 허락을
받아야 합니다.
③ NDSL에서 제공하는 콘텐츠를 보도, 비평, 교육, 연구 등을 위하여 정당한 범위 안에서 공정한 관행에
합치되게 인용할 수 있습니다.
④ NDSL에서 제공하는 콘텐츠를 무단 복제, 전송, 배포 기타 저작권법에 위반되는 방법으로 이용할 경우
저작권법 제136조에 따라 5년 이하의 징역 또는 5천만 원 이하의 벌금에 처해질 수 있습니다.
제 14 조 (유료서비스)
① 당 사이트 및 협력기관이 정한 유료서비스(원문복사 등)는 별도로 정해진 바에 따르며, 변경사항은 시행 전에
당 사이트 홈페이지를 통하여 회원에게 공지합니다.
② 유료서비스를 이용하려는 회원은 정해진 요금체계에 따라 요금을 납부해야 합니다.
제 5 장 계약 해지 및 이용 제한
제 15 조 (계약 해지)
회원이 이용계약을 해지하고자 하는 때에는 [가입해지] 메뉴를 이용해 직접 해지해야 합니다.
제 16 조 (서비스 이용제한)
① 당 사이트는 회원이 서비스 이용내용에 있어서 본 약관 제 11조 내용을 위반하거나, 다음 각 호에 해당하는
경우 서비스 이용을 제한할 수 있습니다.
- 2년 이상 서비스를 이용한 적이 없는 경우
- 기타 정상적인 서비스 운영에 방해가 될 경우
② 상기 이용제한 규정에 따라 서비스를 이용하는 회원에게 서비스 이용에 대하여 별도 공지 없이 서비스 이용의
일시정지, 이용계약 해지 할 수 있습니다.
제 17 조 (전자우편주소 수집 금지)
회원은 전자우편주소 추출기 등을 이용하여 전자우편주소를 수집 또는 제3자에게 제공할 수 없습니다.
제 6 장 손해배상 및 기타사항
제 18 조 (손해배상)
당 사이트는 무료로 제공되는 서비스와 관련하여 회원에게 어떠한 손해가 발생하더라도 당 사이트가 고의 또는 과실로 인한 손해발생을 제외하고는 이에 대하여 책임을 부담하지 아니합니다.
제 19 조 (관할 법원)
서비스 이용으로 발생한 분쟁에 대해 소송이 제기되는 경우 민사 소송법상의 관할 법원에 제기합니다.
[부 칙]
1. (시행일) 이 약관은 2016년 9월 5일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.