3D magnetic inversion, based on the assemblage of 2D forward modeling and inversion as a practical technique to reflect the a priori information, was conducted to investigate the spatial distribution features of black-shale related and pyrometamorphic uranium deposit and several lithological units of Ogcheon Super Group in an area of Geumsan. By using the 3D visualization technique with suitable susceptibility interval and horizontal slice map, the spatial distribution of magnetic susceptibility corresponded to the black shale related uranium bearing lithological units, Black Slate member was well coincided with a information of uranium deposit. Also, even though it is indirect indicator for the detetction of uranium deposits interbedded in Gray Hornfels member, spatial susceptibility distribution which shows the south-east magnetic linearment corresponding to the Majeon-ri formation and Dark Gray Slate were matched well. From this investigation, we inferred that maximum depth extension which Black Slate member can be separately recognized with respect to adjacent Dark Gray Slate with strong magnetic susceptibility anomaly is about 150m with reference elevation level of 306m. In addition, Majeon-ri formation located south of Black Slate member revels relatively high magnetic susceptibility range but shows high spatial susceptibility fluctation. And, as an intrusive rocks, Jurassic Biotite Granite shows relatively low magnetic susceptibility characteristics. On the contrary, Cretaceous granite distributed in soutthern part of the study area shows the relatively high susceptibility distribution.
In this case study, we present the various and consistent processing techniques for the reasonable interpretation of aeromagnetic data. In the processing stage, we especially focused on the three major respects. First, in the low latitude area, severe artifacts are occurred as a result of reduction to the pole technique. To overcome this problem, variable alternative methods were investigated. From the comparison of each technique, we concluded that energy balancing method gives more fruitful result. Second, because of limited a priori information, it is nearly impossible to employ detailed geological survey due to wide and thick spreading of soils in the survey area. So we especially investigated the new techniques such as extracting slope, curvature and aspect information mainly used in GIS field as well as conventional methods. Finally, by using the Euler deconvolution, we extracted the depth information on the magnetic anomalous body. From the synthetic analysis between depth information and previous discussed results, the detailed future survey area was proposed. We think that a series of processing techniques discussed in this study may perform an important role in the domestic and abroad resource development project as a useful guideline.
The simultaneous-source full waveform inversion improves the applicability of full waveform inversion by reducing the computational cost. Since this technique adopts simultaneous multi-source for forward modeling, unwanted events remain in the residual seismograms when the receiver geometry of field acquisition is different from that of numerical modeling. As a result, these events impede the convergence of the full waveform inversion. In particular, the streamer-type data with limited offsets is the most difficult data to apply the simultaneous-source technique. To overcome this problem, the global-correlation-based objective function was suggested and it was successfully applied to the simultaneous-source full waveform inversion in time domain. However, this method distorts residual wavefields due to the modified objective function and has a negative influence on the inversion result. In addition, this method has not been applied to the frequency-domain simultaneous-source full waveform inversion. In this paper, we apply a timedamping function to the observed and modeled data, which are used to compute global correlation, to minimize the distortion of residual wavefields. Since the damped wavefields optimize the performance of the global correlation, it mitigates the distortion of the residual wavefields and improves the inversion result. Our algorithm incorporates the globalcorrelation-based full waveform inversion into the frequency domain by back-propagating the time-domain residual wavefields in the frequency domain. Through the numerical examples using the streamer-type data, we show that our inversion algorithm better describes the velocity structure than the conventional global correlation approach does.
Kim, Rinhui;Chang, Sung-Joon;Mai, Martin;Zahran, Hani
Geophysics and Geophysical Exploration
/
v.22
no.1
/
pp.21-28
/
2019
We perform a 3D tomographic inversion using surface wave dispersion curves to obtain S-velocity model and radial anisotropy beneath Saudi Arabia. The Arabian Peninsula is geologically and topographically divided into a shield and a platform. We used event data with magnitudes larger than 5.5 and epicentral distances shorter than $40^{\circ}$ during 2008 ~ 2014 from the Saudi Geological Survey. We obtained dispersion curves by using the multiple filtering technique after preprocessing the event data. We constructed SH- and SV-velocity models and consequently radial anisotropy model at 10 ~ 60 km depths by inverting Love and Rayleigh group velocity dispersion curves with period ranges of 5 ~ 140 s, respectively. We observe high-velocity anomalies beneath the Arabian shield at 10 ~ 30 km depths and low-velocity anomalies beneath the Arabian platform at 10 km depth in the SV-velocity model. This discrepancy may be caused by the difference between the Arabian shield and the Arabian platform, that is, the Arabian shield was formed in Proterozoic thereby old and cold, while the Arabian platform is covered by predominant Paleozoic, Mesozoic, and Cenozoic sedimentary layers. Also we obtained radial anisotropy by estimating the differences between SH- and SV-velocity models. Positive anisotropy is observed, which may be generated by lateral tension due to the slab pull of subducting slabs along the Zagros belt.
The complex resistivity method is an exploration technique that can obtain various characteristic information of underground media by measuring resistivity and phase in the frequency domain, and its utilization has recently increased. In this paper, a three-dimensional inversion algorithm for the CR data was developed to increase the utilization of this method. The Poisson equation, which can be applied when the electromagnetic coupling effect is ignored, was applied to the modeling, and the inversion algorithm was developed by modifying the existing algorithm by adopting comlex variables. In order to increase the stability of the inversion, a technique was introduced to automatically adjust the Lagrangian multiplier according to the ratio of the error vector and the model update vector. Furthermore, to compensate for the loss of data due to noisy phase data, a two-step inversion method that conducts inversion iterations using only resistivity data in the beginning and both of resistivity and phase data in the second half was developed. As a result of the experiment for the synthetic data, stable inversion results were obtained, and the validity to real data was also confirmed by applying the developed 3D inversion algorithm to the analysis of field data acquired near a hydrothermal mine.
Kim, Dae-sik;Shin, Jungkyun;Ha, Jiho;Kang, Nyeon Keon;Oh, Ju-Won
Geophysics and Geophysical Exploration
/
v.25
no.3
/
pp.109-119
/
2022
The computational efficiency of reverse time migration (RTM) based on numerical modeling is not secured due to the high-frequency band of several hundred Hz or higher for data acquired through a three-dimensional (3D) ultra-high-resolution (UHR) seismic survey. Therefore, this study develops an RTM program to derive high-quality 3D geological structures using UHR seismic data. In the traditional 3D RTM program, an excitation amplitude technique that stores only the maximum amplitude of the source wavefield and a domain-limiting technique that minimizes the modeling area where the source and receivers are located were used to significantly reduce memory usage and calculation time. The program developed through this study successfully derived a 3D migration image with a horizontal grid size of 1 m for the 3D UHR seismic survey data obtained from the Korea Institute of Geoscience and Mineral Resources in 2019, and geological analysis was conducted.
The physics-informed neural network (PINN) has been proposed to overcome the limitations of various numerical methods used to solve partial differential equations (PDEs) and the drawbacks of purely data-driven machine learning. The PINN directly applies PDEs to the construction of the loss function, introducing physical constraints to machine learning training. This technique can also be applied to wave equation modeling. However, to solve the wave equation using the PINN, second-order differentiations with respect to input data must be performed during neural network training, and the resulting wavefields contain complex dynamical phenomena, requiring careful strategies. This tutorial elucidates the fundamental concepts of the PINN and discusses considerations for wave equation modeling using the PINN approach. These considerations include spatial coordinate normalization, the selection of activation functions, and strategies for incorporating physics loss. Our experimental results demonstrated that normalizing the spatial coordinates of the training data leads to a more accurate reflection of initial conditions in neural network training for wave equation modeling. Furthermore, the characteristics of various functions were compared to select an appropriate activation function for wavefield prediction using neural networks. These comparisons focused on their differentiation with respect to input data and their convergence properties. Finally, the results of two scenarios for incorporating physics loss into the loss function during neural network training were compared. Through numerical experiments, a curriculum-based learning strategy, applying physics loss after the initial training steps, was more effective than utilizing physics loss from the early training steps. In addition, the effectiveness of the PINN technique was confirmed by comparing these results with those of training without any use of physics loss.
A field technique for assessing the transmissive fractures in an aquifer was applied to a fractured rock formation in Youngchun area Korea. Geological mapping and detailed acoustic borehole teleview(BHTV) logging were performed to obtain information about the fractures. The study area consists predominantly of two types of fractures. The fracture sets of low angle partings such as bedding and sheeting plains have strike N70-80$^{\circ}$W, 25$^{\circ}$-30$^{\circ}$SW and N3S$^{\circ}$W, 12$^{\circ}$NE, respectively. In areas of high fractures, on the other hand, the major fracture sets show strike N80$^{\circ}$W and dip 70$^{\circ}$-85$^{\circ}$SW, N10$^{\circ}$E.85$^{\circ}$SE in sedimentry rocks, N40-50$^{\circ}$E.85$^{\circ}$SE/85$^{\circ}$NE, N70$^{\circ}$E.80$^{\circ}$SE, and N7$^{\circ}$-75$^{\circ}$W.80$^{\circ}$SW in granites and volcanic rocks. Injection tests have been performed to identify discrete production zones and quantify the vertical distribution of hydraulic conductivity. The calculated hydraulic conductivities range from 3.363E-10 to 2.731E-6, showing that the difference between maximum and minimum value is four order of magnitude. Dominant section in hydraulic conductivity is extensively fractured. Geophysical logging was carried out to clarify characterization of the distribution of fracture zones. Transmissive fractures were evaluated through the comparison of the results obtained by each method. The temperature logs appeared to be a good indicator that can distinguish a high transmissive fractures from a common fractures in hydraulic conductivity. In numerous cases, evidence of fluid movement was amplified in the temperature gradient log. The fracture sets of N70-80$^{\circ}$W.60-85$^{\circ}$NE/SW N75-80$^{\circ}$W.25-30$^{\circ}$SW, N50-64$^{\circ}$W.60-85$^{\circ}$NE, N35-45$^{\circ}$E.65-75$^{\circ}$SE, and N65-72$^{\circ}$E.80$^{\circ}$SE/60$^{\circ}$NW were idenfied as a distinct transmissive fractures through the results of each tests.
This study is the case of groundwater development based on the geometrical fracture model of target area established only through geological fracture mapping technique. A fracture mapping of $9km^2$, eastern Munsan, has been conducted to determine geological and hydrological factors for new water well placement in the Gyeonggi gneiss complex. Geophysical exploration was not applicable because of small restricted area and dense underground utilities at the site. Form line mapping on the basis of foliation orientation and rock type revealed a synform of NS fold axis bearing to the south. An EW geological cross-section passed through the site area shows a F2 synform as a double-wall ice cream spoon shape. Three regional faults of $N20^{\circ}E,\;N30^{\circ}W$, and NS have been dragged into the site to help understand extensional fault paths. The $N20^{\circ}E$ fault with dextral sense is geometrically interpreted as a western fault of two flexural conjugate type-P shear faults in the F2 synformal fold. The NE cross-section reveals that a possible groundwater belt in the western limb of super-posed fold area is formed as a trigonal prism within 100 m depth of the intersectional space between the $N20^{\circ}E$ fault plane and the weakly sheared plane of transposed foliation. Another possible fault for water resource strikes $N40^{\circ}E$. Recommended sites for new water well placement are along the $N20^{\circ}E\;and\;N40^{\circ}E$ faults. As a result of fracture mapping, 145 ton/day of water can be produced at one well along the $N20^{\circ}E$ fault line. Exploration of groundwater in the area is succeeded only using with geological fracture mapping and interpretation of geological cross-section, without any geophysical survey. Intersection of fault generated with the F2 synformal fold and foliation supply space of groundwater reserver.
Vertical electrical sounding and 2D electrical resistivity survey were applied for evaluating the characteristics of alluvial layers at a groundwater artificial recharge site. The fine particles in alluvial layer, main target layer of groundwater artificial recharge, may cause clogging phenomena. In this case, electrical resistivity method is an effective technique to verify the spatial distribution of low-resistivity layers, such as saturated silts and clays. On the other hand, much attention should be paid to interpret the resistivity data in unconsolidated layers, because thick clayey overburden sometimes produces a masking effect on underlying interbedded resistive sands and gravels. Considering these points, we designed 35 points arranged in a grid form for vertical electrical sounding and 10 lines for 2D electrical resistivity survey, and concentrated our effort on enhancing the vertical and horizontal resolution of resistivity images. According to the results, 15 meters thick layers consisting of sands and gravels are located in 30 meters below ground. And the spatial distribution of silts and clays are mapped, which may cause clogging. Consequently, this approach can contribute to design and determine the location and depth of injection and observation wells for groundwater artificial recharge.
본 웹사이트에 게시된 이메일 주소가 전자우편 수집 프로그램이나
그 밖의 기술적 장치를 이용하여 무단으로 수집되는 것을 거부하며,
이를 위반시 정보통신망법에 의해 형사 처벌됨을 유념하시기 바랍니다.
[게시일 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일부터 적용되며, 종전 약관은 본 약관으로 대체되며, 개정된 약관의 적용일 이전 가입자도 개정된 약관의 적용을 받습니다.