한국 남동해역 내대륙붕 이토대 퇴적물의 물리적 성질

김길영,김대철,서영교,박수철,최진혁,김정창,Kim, Gil-Young,Kim, Dae-Choul,Seo, Young-Kyo,Park, Soo-Chul,Choi, Jin-Hyuk,Kim, Jeong-Chang 한국해양학회 1999 바다 Vol.4 No.4

한국 남동해역 내대륙붕 이토대(mudbelt) 퇴적물의 물리적 성질을 연구하기 위하여 이토대지역 총 14개 정점에서 피스톤시추기를 이용하여 해저 퇴적물을 채취하였다. 시추퇴적물은 실험실에서 각 코어별 일정한 깊이로 퇴적물의 조직, 물리적 성질 및 음파전달속도를 측정하였다. 표층퇴적물의 조직 중 평균입도는 연구지역의 남쪽지역인 울산 앞바다에서 북쪽으로 갈수록 등수심선과 거의 평행하게 감소하는 경향을 보인다. 물리적 성질인 습윤전밀도와 공극율 그라고 속도 등의 분포 형태도 평균입도 분포와 유사하게 변화하는 경향을 뚜렷하게 보여주어 퇴적물의 입도가 물성에 큰 영향을 끼치고 있음을 알 수 있다. 이러한 퇴적물의 조직 및 물성의 분포특성은 낙동강 및 한반도 주변에서 공급된 세립질 퇴직물이 연안류나 대마난류 등에 의해 해안선을 따라 북동쪽으로 이동하면서 분급화되어 북쪽으로 갈수록 점진적으로 조립질에서 세립질로 변해가는 퇴적과정을 보여주는 것으로 해석된다. 퇴적물의 조직, 물성 및 속도등의 수직적인 변화도 최적후의 다짐작용보다는 퇴적물의 조직에 의해 영향을 받고 있는 것으로 나타났다. 조직 및 물성간의 상관관계를 보면 비교모델로 시용한 북태평양 대륙붕 및 사면퇴적물은 물론 남해대륙붕 퇴적물과 전반적인 경향은 유사하나 절대값의 차이가 있다. 즉 같은 조직 및 물성값에 대해 본 연구지역의 속도가 더 높다. 이러한 결과는 퇴적물의 기원에 기인한 광물조성 및 퇴적환경의 차이와 본 연구지역에 비교적 흔한 가스함유퇴적층의 존재로 인한 가스의 방출로 퇴적물의 공극비가 감소하여 높은 속도감을 보이는 것으로 사료된다. Physical properties of mudbelt sediments in the southeastern inner shelf of Korea are studied from 14 cores. Physical properties, compressional wave velocity, and sediment texture for core sediments are analyzed. The major source of sediment in the study area is the Nakdong River. Fine-grained sediments from the river are transported northeastward by coastal circulation and the Tsushima Current, resulting in a gradual northeastward increase in porosity and a decrease in wet bulk density and velocity. The trend matches well with the bathymetry. The mean grain size appears to be the most important variable to determine the physical properties and velocity. The variations of physical properties with burial depth are dependent more strongly on sediment texture than compaction and/or consolidation. Correlations between the physical properties and the sediment texture show slight deviations from those of the continental terrace sediment in the North Pacific and inner shelf sediment in the South Sea of Korea. The velocity is higher than that of the North Pacific and the South Sea sediments between these areas. This is probably due to differences in sedimentary, environment and mineral compositions. The higher sediment velocity in the study area may also be attributed to the escape of gas from pore space which decreases void ratio.

Acoustic Properties of Gassy Sediments: Preliminary Result of Jinhae Bay, Korea

김길영,김대철,여정윤,유동근,Kim, Gil-Young,Kim, Dae-Choul,Yeo, Jung-Yoon,Yoo, Dong-Geun The Acoustical Society of Korea 2007 韓國音響學會誌 Vol.26 No.e1

Compressional wave velocity and shear wave velocity were measured for gassy sediments collected from Jinhae Bay, Korea. To distinguish inhomogeneities of gassy sediments, Computed Tomography (CT) was carried out for gassy sediment using CT Scanner. The cored sediments are composed of homogeneous and soft mud (greater than $8{\Phi}$ in mean grain size) containing clay content more than 50%. In depth interval of gassy sediments, compressional wave velocity is significantly decreased from 1480m/s to 1360m/s, indicating that the gas greatly affects compressional wave velocity due to a gas and/or degassing cracks. Shear wave velocity shows a slight increasing pattern from ${\sim}55\;m/s$ in the upper part of the core to ${\sim}58\;m/s$ at 320 cm depth, and then decreases to ${\sim}54\;m/s$ in the lower part of the core containing a small amount of gas. But shear wave velocity in the gassy sediments is slightly greater than that of non-gassy sediments in the upper part of the core. Thus, the Vp/Vs ratio is decreased (from 30 to 25) in gas charged zone. The Vp/Vs ratio is well correlated with shear wave velocity, but no correlation with compressional wave velocity. This suggests that low concentrations of gas have little affects on shear wave velocity. By CT images, the gas in the sediments is mostly concentrated around inner edge of core liner due to a long duration after sediment collection.

거제도 동쪽 해역 표층 퇴적물의 지음향모델

김길영,김대철,신보경,서영교,이광훈,KIM, GIL YOUNG,KIM, DAE CHOUL,SHIN, BO KYOUNG,SEO, YOYUNG KYO,LEE, GWANG HOON 한국해양학회 2005 바다 Vol.10 No.2

Sediment texture, physical (porosity, water content, bulk density, grain density, and shear strength), and geoacoustic properties (compressional wave velocity and attenuation) were measured on eighteen core samples collected from the shelf off eastern Geoje Island, the South Sea of Korea. Based on these properties, the study area is divided into three different sub-areas: (1) Area I affected directly by the Nakdong River discharge; (2) Area II covered by the southern branch of the Nakdong River discharge; and (3) Area III dominated by relict sediment. Mean grain size, velocity, and bulk density decrease from Area $I(7.4\Phi,\;1528m/s,\;1.6g/cm^3,\;respectively)$ to Area $II(8.1\Phi,\;1485m/s,\;and\;1.5g/cm^3)$, and then increase rather rapidly in Area $III(1.4\Phi,\;1664m/s,\;and\;2.2g/cm^3)$. Porosity, on the other hand, exhibits an opposite trend, increasing from Area $I(64.5\%)$ to Area$II(73.9\%)$ and then decreasing significantly in Area $III(32.9\%)$ From the results measured and calculated, we suggest a specified geoacoustic model in the study area. 대한해협과 거제도 동쪽 사이에 분포하고 있는 표층퇴적물의 지음향모델을 구현하기 위하여 18개 정점에서 채취한 시추 코어를 대상으로 퇴적물 조직, 물리적 성질(공극률, 함수율, 전밀도, 입자밀도, 전단응력) 및 지음향 특성(음파전달속도, 음파감쇠)을 측정하였다. 평균입도, 음파전달속도 및 전밀도는 연구지역의 외해로 갈수록 증가하며, 공극률은 감소하는 경향을 보였다. 퇴적물 특성을 통해 본 연구지역을 3개의 구역(구역 I: 낙동강으로부터 직접 퇴적물이 퇴적되는 지역:구역 II:낙동강 부유퇴적물이 남쪽 방향으로 이동 및 퇴적된 니질 퇴적물 지역; 구역 III: 조립질 잔류퇴적물이 퇴적되어 있는 지역)으로 구분 할 수 있었다. 평균입도, 음파전달속도 및 전밀도는 구역 I(각각 $7.4\Phi,\;1528m/s,\;1.6g/cm^3$)에서 구역 II(각각 $8.1\Phi,\;1485m/s,\;1.5g/cm^3$)로 줄어들다가 구역 $III(1.4\Phi,\;1664m/s,\;2.2g/cm^3)$으로 급격히 증가하는 경향을 보인다 반면, 공극률은 구역 $I(64.5\%)$에서 구역 $II(73.9\%)$로 증간하다 구역 $III(32.9\%)$으로 감소하는 반대의 경향을 보인다. 따라서 본 연구에서는 이러한 구역별 분석 및 측정된 자료와 계산된 음향특성 값을 이용하여 연구지역을 구역별로 나누어진 세분화된 지음향모델을 제시하고자 한다.

동해 울릉분지 가스하이드레이트 함유 퇴적물의 음파전달속도 특성

김길영,유동근,류병재,Kim, Gil-Young,Yoo, Dong-Geun,Ryu, Byong-Jae 한국음향학회 2009 韓國音響學會誌 Vol.28 No.5

This study investigates the difference of sound velocity (compressional wave velocity) between gas hydrate-bearing sediments and nongas hydrate-bearing sediments in the Ulleung Basin, East Sea. We use a dataset measured from one site in the central part of the Ulleung Basin. Sound velocity for gas hydrate-bearing sediment shows the range from 1600 m/s to 2200 m/s. However, the value for nongas hydrate-bearing sediment is mostly around 1500 m/s, being less than 1400 m/s below 140 m subbottom depth. This trend is probably due to the presence of free gas below BSR (Bottom Simulating Reflector). Gas hydrate-bearing sediments show high value (maximum 150 Ohm-m) of resistivity. The physical properties between gas hydrate-bearing sediment and nongas hydrate-bearing sediment are characterized by the different patterns due to the presence of gas hydrate in comparison with those of marine unconsolidated sediments. Therefore, in order to investigate acoustic and physical properties for gas hydrate-bearing sediments, the study for the occurrence type and the amount of gas hydrates should be conducted simultaneously. 본 연구는 지구물리검층을 통해서 동해 울릉분지의 가스하이드레이트 함유지층과 비함유지층에서의 음파전달속도 (종파)의 특성차이를 규명하였다. 지구물리 검층 자료는 울릉분지의 중앙지역에서 취득하였다. 음파전달속도는 가스하이드레이트 함유지층의 경우 최대 속도값이 약 2200 m/s의 높은 값에서 1600 m/s의 낮은 값을 갖는다. 반면에 비함유지층에서는 약 1500 m/s 내외에서 1400 m/s보다 더 낮은 값을 갖기도 한다. 일반 해양퇴적물에서의 값 (약 $1500{\sim}1600$ m/s)보다 높은 값을 보이는 것은 가스하이드레이트가 퇴적물내에 함유되어 있기 때문이다. 가스하이드레이트 비함유구간중 해수의 속도보다 낮은 값을 보이는 구간은 해저면에서 약 140 m 이하의 해저모방반사면 (Bottom Simulating Reflector)의 하부층에 존재하는 자유가스 때문이다. 전기비저항값도 가스하이드레이트 함유구간에서 최대 150 Ohm-m까지 높게 나타나 비함유구 간에서의 값과 차이가 크다. 각 물성간의 상관관계를 보면 가스하이드레이트 함유지층과 비함유지층의 구분이 명확하게 나타나고 있으며 그 상관관계도 일반적인 해양 미고결퇴적물에서 보이는 값과는 상이한 양상을 보여주고 있다. 따라서 가스하이드레이트가 존재하는 지층의 물성과 음향특성을 해석할 때에는 가스하이드레이트의 존재 형태나 함량 등에 대한 연구가 병행되어야 할 것으로 생각된다.

해양 과학시추 50년 (1968-2018): 한국의 성과 및 미래 방향

김길영,KIM, GIL YOUNG 한국해양학회 2019 바다 Vol.24 No.1

2018년은 해양 과학시추가 시작된 지 50년이 되는 해이다. 그럼에도 불구하고 우리는 지구의 대양저보다 달의 표면을 더 많이 알고 있을 정도로 지구내부에 대한 정보가 많지 않다. 대양저에 관한 연구는 해양 과학시추로부터 얻어진 시료를 통해서 알 수 있다. 이러한 심해의 시료획득은 50년 전인 1968년 8월 11일 미연방차원에서 지원된 심해저시추계획(DSDP: Deep Sea Drilling Project)에서 글로마 챌린저(Glomar Challenger)호를 이용한 멕시코 만 시추로부터 시작되었다. 이후 해저지각시추프로그램(ODP: Ocean Drilling Program), 통합해저지각시추프로그램(old IODP: Integrated Ocean Drilling Program), 그리고 국제해양시추탐사프로그램 (new IODP: International Ocean Discovery Program)으로 이어져 오고 있다. 해양 과학시추로부터 얻어진 가장 큰 성과는 두 가지 기술적인 성과와 다양한 과학적인 성과로 나눠진다. 첫 번째 기술적인 성과는 시추선이 시추위치를 벗어나지 않고 연속적으로 시추코어를 획득할 수 있도록 위치를 유지시켜주는 동적위치유지 시스템(dynamic positioning system)이다. 다른 하나는 시추동안 드릴 비트를 교체한 후 동일한 시추공에 드릴비트가 쉽게 투입될 수 있도록 해주는 재투입 콘(re-entry cone)의 개발이다. 이러한 기술적인 혁신 외에도 다양한 과학적 성과 즉 판구조론 증명, 지구의 역사 규명, 그리고 심해 퇴적물 내 생명체의 발견 등이 있다. 2013년 10월 시작된 국제해양시추탐사프로그램(new IODP)은 2023년까지 계속될 예정이고, 2023년 이후 다음 단계를 위해 참여 회원국들은 새로운 과학계획 수립과 더불어 미래의 해양 과학시추 50년을 준비하고 있다. 우리나라도 이러한 국제적인 동향에 발맞추어 회원국가로서 지속적인 참여와 다음단계를 위한 준비가 필요한 시점이다. The year 2018 is the $50^{th}$ anniversary of scientific ocean drilling. Nevertheless, we know more about the surface of the moon than the Earth's ocean floor. In other words, there are still no much informations about the Earth interior. Much of what we do know has come from the scientific ocean drilling, providing the systematic collection of core samples from the deep seabed. This revolutionary process began 50 years ago, when the drilling vessel Glomar Challenger sailed into the Gulf of Mexico on August 11, 1968 on the first expedition of the federally funded Deep Sea Drilling Project (DSDP). DSDP followed successively by Ocean Drilling Program (ODP), Integrated Ocean Drilling Program (old IODP), and International Ocean Discovery Program (new IODP). Concerning on the results of scientific ocean drilling, there are two technological innovations and various scientific research results. The one is a dynamic positioning system, enables the drilling vessel to stay fixed in place while drilling and recovering cores in the deep water. Another is the finding of re-entry cone to replace drill bit during the drilling. In addition to technological innovation, there are important scientific results such as confirmation of plate tectonics, reconstruction of earth's history, and finding of life within sediments. New IODP has begun in October, 2013 and will continue till 2023. IODP member countries are preparing for the IODP science plan beyond 2023 and future 50 years of scientific ocean drilling. We as IODP member also need to participate in keeping with the international trend.

초음파 열분해법를 이용한 ZnO 성장

김길영,정연식,변동진,최원국,Kim, Gil-Young,Jung, Yeon-Sik,Byun, Dong-Jin,Choi, Won-Kook 한국세라믹학회 2005 한국세라믹학회지 Vol.42 No.4

단결정 사파이어 (0001) 기판 위에 저가의 초산아연(Zinc Acetate Dehydrate; ZAH) 전구체를 이용하여 초음파 열분해법과 Ar 가스를 이용한 ZnO 박막을 성장시켰다. Thermogravimetry-Differential Scanning Calorimetry(TG-DSC) 초산아연의 열분해 과정을 조사하여 $380^{\circ}C$ 이상에서 ZnO로 분해되는 것을 확인하였다. $380-700^{\circ}C$에서 증착된 ZnO 박막은 모두 ZnO (002), (101) 결정면으로 부터의 회절피크를 보여주고 있었으며, $400^{\circ}C$ 박막의 경우 c-압축 스트레인 ${\Sigma}Z=0.2\%$, 압축 응력 $\sigma=-0.907\;GPa$이 작용하고 있음을 알 수 있었다. 전자 현미경을 이용한 미세 구조의 관찰을 통하여 $380-600^{\circ}C$에서는 초산아연과 ZnO 초미세 입자가 혼합된 aggregate 형태의 결정립을 형성하고 있었으며, nanoblade 형태의 미세구조를 보였다. 한편 $700^{\circ}C$에서 증착된 박막내의 결정립은 찌그러진 육방정계의 형태를 취하고 있으며, 10-25nm 정도의 부결정림 초미세 ZnO 입자로 이루어져 있음을 알 수 있었다. 초미세 입자의 형성을 임의 핵형성 기구(random nucleation mechanism)로 설명하였고, photoluminescence(PL) 측정을 통하여 광 특성을 조사하였다. ZnO was deposited on sapphire single crystal substrate by an ultrasonic pyrolysis of Zinc Acetate Dehydrate (ZAH) with carrying Ar gas. Through Thermogravimetry-Differential Scanning Calorimetry(TG-DSC), zinc acetate dihydrate was identified to be dissolved into ZnO above $380^{\circ}C$. ZnO deposited at $380-700^{\circ}C$ showed polycrystalline structures with ZnO (101) and ZnO (002) diffraction peaks like bulk ZnO in XRD, and from which c-axis strain ${\Sigma}Z=0.2\%$ and compressive biaxial stress$\sigma=-0.907\;GPa$ was obtained for the ZnO deposited $400^{\circ}C$. Scanning electron microscope revealed that microstructures of the ZnO were dependent on the deposition temperature. ZnO grown below temperature $600^{\circ}C$ were aggregate consisting of zinc acetate and ZnO particles shaped with nanoblades. On the other hand the grain of the ZnO deposited at $700^{\circ}C$ showed a distorted hexagonal shape and was composed of many ultrafine ZnO powers of 10-25 nm in size. The formation of these ulrafine nm scale ZnO powers was explained by the model of random nucleation mechanism. The optical property of the ZnO was analyzed by the photoluminescence (PL) measurement.

IODP Exp. 346 동해 과학시추와 부산항 Port Call

김길영(Gil Young Kim),장세원(Se Won Chang),현상민(Sangmin Hyun) 대한지질학회 2013 대한지질학회 학술대회 Vol.2013 No.10

청소년 탄력성이 자살충동에 미치는 영향과 자살예방에 관한 연구

김길영(Kim, Gil-Young),조성제(Cho, Sung-Je) 한국산학기술학회 2013 한국산학기술학회논문지 Vol.14 No.5

본 연구는 청소년 탄력성과 자살충동과의 상관성을 명확하게 파악하여, 외적보호요인을 강화하는 보다 효과 적인 청소년 자살예방 대안 제시를 목적으로 한다. 연구대상은 경기도에 소재하는 일반고 고등학교 1학년 남녀 학생 243명이며 2013년 3월 5일부터 1주일간 설문조사를 실시하였다. 분석방법은 빈도분석, 상관분석, 다중회귀분석, t-test, 일원변량분석 등을 실시하였고, 실증분석은 모두 유의수준 5%에서 검증하였다. 연구결과는 자살충동과 탄력성은 부의 상관관계가 있었으며, 특히 가정차원의 상호작용이 높을수록, 또래차원과 학교차원의 돌봄과 기대가 높을수록 자살충 동이 낮았다. 본 연구의 시사점은 청소년 탄력성이 자살충동에 미치는 영향에 대한 연구결과를 바탕으로 자살예방정 책마련을 위한 자료가 될 것으로 사료된다. This study clearly identify the correlation between adolescents resiliency and suicidal ideation and more effective to strengthen the external protective factors for adolescents suicide prevention alternatives presented. Research subjects were 243 boy and girl first year students in a high school located in Kyunggi Province; Questionnaires were distributed for one week from March 5, 2013. Frequency analysis, correlation analysis, multiple regression analysis, t-test, one-way ANOVA, etc. were executed. All empirical tests were verified with 5% of significance level. From the result of analysis, suicidal ideation has a negative correlation with the adolescent resilience, and especially, when interaction in family dimension and the care, expectation level of school and peer dimension become higher, suicidal ideation was reduced. The implications of this study suggest that the material prepared for suicide prevention policy based on the results of a study on the impact of adolescents resiliency to suicidal ideation

종이 헬리콥터 낙하해석모델의 통계적 교정 및 검증

김길영(Gil Young Kim),유성범(Sung Bum Yoo),김동영(Dong Young Kim),김동성(Dong Seong Kim),최주호(Joo Ho Choi) 대한기계학회 2015 大韓機械學會論文集A Vol.39 No.8

수학적 해석모델은 물리적 현상을 파악하고 실험비용을 절감하는데 활발하게 사용되지만 편의를 위한 단순화 또는 파라미터가 가지고 있는 불확실성에 의해 해석모델에 의한 예측결과는 실제현상과 차이가 발생한다. 본 연구에서는 이러한 문제에 대해 통계적 기법을 이용하여 해석모델의 불확실성을 반영한 교정 및 검증 방법을 종이 헬리콥터를 통해 제시한다. 먼저, 같은 제원의 세 가지 종이 헬리콥터로 실시한 실험 데이터를 각 그룹으로 형성하여 두 가지 낙하해석모델에서 미지의 입력 파라미터인 항력계수를 교정하는데 사용했다. 그리고 확률분포로 예측된 낙하시간을 실험 데이터 분포와 비교하여 해석 모델을 검증하였다. 이 때, Markov Chain Monte Carlo 기법을 활용하여 항력계수의 불확실성을 정량화하였다. 또한 종이 헬리콥터의 그룹별 데이터에 대해 분산분석(Analysis of Variance)를 이용하여 제작오차와 실험오차의 관계를 비교하였고, 각 그룹이 모두 동일한 대상으로 간주해도 됨을 증명하였다. Mathematical models are actively used to reduce the experimental expenses required to understand physical phenomena. However, they are different from real phenomena because of assumptions or uncertain parameters. In this study, we present a calibration and validation method using a paper helicopter and statistical methods to quantify the uncertainty. The data from the experiment using three nominally identical paper helicopters consist of different groups, and are used to calibrate the drag coefficient, which is an unknown input parameter in both analytical models. We predict the predicted fall time data using probability distributions. We validate the analysis models by comparing the predicted distribution and the experimental data distribution. Moreover, we quantify the uncertainty using the Markov Chain Monte Carlo method. In addition, we compare the manufacturing error and experimental error obtained from the fall-time data using Analysis of Variance. As a result, all of the paper helicopters are treated as one identical model.

해설 : 미고결 퇴적층내 가스하이드레이트 포화도 계산

김길영 ( Gil Young Kim ) 한국지구물리·물리탐사학회 2012 지구물리와 물리탐사 Vol.15 No.2

이 논문은 퇴적물내 포함되어 있는 가스하이드레이트의 포화도를 계산하는 여러 방법에 대하여 토론하고자 한다. 가스하이드레이트의 포화도를 계산하는 방법은 물리검층 자료를 이용하는 방법과 코어자료(압력코어 포함)를 이용하는 방법, 그리고 탄성파 탐사자료로부터 얻을 수 있는 속도 자료를 이용하는 방법 등 크게 세가지 방법으로 나눌 수 있다. 물리검층 자료중 전기비저항 자료를 이용하는 방법의 경우 Archie 식을 주로 이용하는데 이 경우 각각의 변수 값을 정확하게 정의하는 게 중요하다. 또한 가스하이드레이트의 산출형태도 포화도 계산에 큰 영향을 주기 때문에 주의해야 한다. 코어자료를 이용하는 경우 공극수의 염소량을 측정하는 방법과 압력코어를 취득할 경우 이를 이용하는 방법이 있다. 지금까지 발표된 정량적이고 가장 정확한 가스하이드레이트 포화도값을 구할 수 있는 방법이 압력코어를 이용하는 것이다. 그러나 이는 비용과 시간이 많이 소요되기 때문에 연속적인 자료를 얻기가 어렵다는 단점이 있다. 지금까지 발표된 가스하이드레이트 포화도 값을 비교해 보면 전기비저항 값을 이용한 경우가 가장 높은 값을 압력코어를 이용하여 측정한 경우가 가장 낮은 값을 보여주는 경향이 있다. 그러나 이러한 값이 모든 경우에 있어서 절대적인 경향을 보여준다고 볼 수는 없다. 그러므로 가스하이드레이트의 포화도를 정확하게 계산하기 위해서는 여러 가지 방법을 이용하여 계산해야 하며 이를 비교하여 가장 적절한 값을 사용해야 할 것이다. The purpose of this paper is to review several different methods calculating gas hydrate saturations. There are three methods using downhole log data, core data (including pressure core), and seismic velocity data. Archie`s equation using electrical resistivity of downhole log data is widely used for saturation calculation. In this case, Archie`s parameters should be defined accurately. And the occurrence types of gas hydrate significantly affect to saturation calculation. Thus saturation calculation should be carefully conducted. The methods using chlorinity and pressure core data are directly calculated from core sample. So far, the saturation calculated from pressure core gives accurate and quantitative values. But this method is needed much more time and cost. Thus acquisition of the continuous data with sediment depth is realistically hard. The recent several results show that the saturation calculated from resistivity data is the highest values, while the value calculated from pressure core is the lowest. But this trend is not always absolutely. Thus, to estimate accurate gas hydrate saturation, the values calculated from several methods should be compared.