JP6872762B2 - 堆積物の岩石コアにおけるイベント層の識別方法 - Google Patents
堆積物の岩石コアにおけるイベント層の識別方法 Download PDFInfo
- Publication number
- JP6872762B2 JP6872762B2 JP2020142431A JP2020142431A JP6872762B2 JP 6872762 B2 JP6872762 B2 JP 6872762B2 JP 2020142431 A JP2020142431 A JP 2020142431A JP 2020142431 A JP2020142431 A JP 2020142431A JP 6872762 B2 JP6872762 B2 JP 6872762B2
- Authority
- JP
- Japan
- Prior art keywords
- layer
- event
- rock core
- ratio
- depth
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000011435 rock Substances 0.000 title claims description 86
- 239000011575 calcium Substances 0.000 claims description 57
- 239000013049 sediment Substances 0.000 claims description 48
- 238000004458 analytical method Methods 0.000 claims description 42
- 238000000034 method Methods 0.000 claims description 40
- 230000008859 change Effects 0.000 claims description 21
- 238000010606 normalization Methods 0.000 claims description 17
- 229910052791 calcium Inorganic materials 0.000 claims description 11
- 235000013619 trace mineral Nutrition 0.000 claims description 6
- 239000011573 trace mineral Substances 0.000 claims description 6
- 229910052701 rubidium Inorganic materials 0.000 claims description 5
- 229910052712 strontium Inorganic materials 0.000 claims description 5
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 claims description 3
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 3
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 claims description 3
- 229910052782 aluminium Inorganic materials 0.000 claims description 3
- IGLNJRXAVVLDKE-UHFFFAOYSA-N rubidium atom Chemical compound [Rb] IGLNJRXAVVLDKE-UHFFFAOYSA-N 0.000 claims description 3
- 229910052710 silicon Inorganic materials 0.000 claims description 3
- 239000010703 silicon Substances 0.000 claims description 3
- CIOAGBVUUVVLOB-UHFFFAOYSA-N strontium atom Chemical compound [Sr] CIOAGBVUUVVLOB-UHFFFAOYSA-N 0.000 claims description 3
- 229910052726 zirconium Inorganic materials 0.000 claims description 3
- QCWXUUIWCKQGHC-UHFFFAOYSA-N Zirconium Chemical compound [Zr] QCWXUUIWCKQGHC-UHFFFAOYSA-N 0.000 claims description 2
- 239000002245 particle Substances 0.000 description 7
- 230000002159 abnormal effect Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011946 reduction process Methods 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 238000013517 stratification Methods 0.000 description 2
- 238000006424 Flood reaction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 235000010755 mineral Nutrition 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- VSZWPYCFIRKVQL-UHFFFAOYSA-N selanylidenegallium;selenium Chemical compound [Se].[Se]=[Ga].[Se]=[Ga] VSZWPYCFIRKVQL-UHFFFAOYSA-N 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N23/00—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00
- G01N23/02—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material
- G01N23/04—Investigating or analysing materials by the use of wave or particle radiation, e.g. X-rays or neutrons, not covered by groups G01N3/00 – G01N17/00, G01N21/00 or G01N22/00 by transmitting the radiation through the material and forming images of the material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/03—Investigating materials by wave or particle radiation by transmission
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/10—Different kinds of radiation or particles
- G01N2223/101—Different kinds of radiation or particles electromagnetic radiation
- G01N2223/1016—X-ray
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/30—Accessories, mechanical or electrical features
- G01N2223/33—Accessories, mechanical or electrical features scanning, i.e. relative motion for measurement of successive object-parts
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2223/00—Investigating materials by wave or particle radiation
- G01N2223/60—Specific applications or type of materials
- G01N2223/616—Specific applications or type of materials earth materials
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Description
本発明は、海洋堆積学におけるイベント層識別及び古気候研究分野に属し、具体的には海域堆積物の岩石コアにおけるイベント層の識別方法に関するものである。
異常気象イベント（例えば台風や洪水など）は、異常な降雨の発生を招き、これによって海洋中に輸送される物質が突然に増加し、最終的に海洋中でこれらの物質の堆積が生じることで海洋堆積物が形成される。従って、海洋堆積物中の異常堆積層を抽出し、且つ室内での実験解析結果と結合することで、異常気象イベントの発生履歴を復元することができ、将来の世界各地における異常気象イベントの発生を予測するうえで重要な科学的根拠が提供される。また、岩石コアにおけるイベント層の存在は、岩石コアの年代較正の正確性にも影響するため、岩石コアにおけるイベント層の識別は、海洋堆積学研究の発展にとって極めて重要である。
現在、海洋堆積物の岩石コアからイベント層を抽出する方法は主に、１）堆積物の岩石コアを２つに分割し、岩石コアの堆積層全体における異質性を肉眼で観察し、イベント層を判断する方法と、２）岩石コアを等間隔に従ってサンプリングした後、実験室内で粒度と元素を解析し、測定結果によってイベント層を確定する方法がある。但し、これら２種類の堆積物の岩石コアにおけるイベント層の識別方法には、欠点もある程度存在する。１）岩石コアの堆積層の異質性を肉眼的に観察すると、人によって違いが生じ、イベント層の正確な位置を確定することができない。２）等間隔でサンプリングして室内解析を行う場合、試料の縮分間隔が限定され、極めて小さなイベント層を正確に識別することができない。また、当該方法で実験結果を得るには期間が長く、時間と手間がかかる。
よって、海洋における古気候及び古環境研究の発展に有力な技術支援を提供するために、堆積物の岩石コアにおけるイベント層を迅速且つ精確に識別する方法の構築が期待されている。
本発明は、従来技術の堆積物の岩石コアにおけるイベント層を識別する方法に存在する精度の悪さなどの欠点に鑑み、岩石コア中の微小なイベント層を識別できるだけでなく、識別効率を大幅に高め、イベント層識別の正確性を大きく向上させた、堆積物の岩石コアにおけるイベント層の識別方法を提出する。
本発明は、以下の技術案を採用して実現したものである。堆積物の岩石コアにおけるイベント層の識別方法であって、以下の工程を含む。
工程Ａ、堆積物の岩石コアを分割し、Ｘ線により堆積物の岩石コアの元素に対する走査を行い、堆積物の岩石コアに含まれる元素の強度データを得る。
工程Ｂ、得られた元素の強度データを基に、堆積層における元素の賦存特性に従って堆積層の解析指数を確定し、解析指数は元素の強度データによって表す。
工程Ｃ、工程Ｂで確定した堆積層の解析指数を基に、その元素の強度データに対して解析を行うことでイベント層の識別が実現するが、それには以下を含む。
（Ｃ１）解析指数の平均値
及び標準偏差δを確定する。
（Ｃ２）解析指数を正規化した後の正規化標準偏差ＳＤを基に、イベント層を識別する閾値を確定する。正規化後の標準偏差ＳＤは、解析指数の元素強度データに対して正規化処理を行って得たものであり、閾値とは、堆積物の岩石コアにおける正常堆積層とイベント層を区分する臨界値のことをいう。
（Ｃ３）正規化標準偏差ＳＤの値に従って閾値を確定し、正規化後の解析指数の深度別変化曲線に対して、確定した閾値に基づく深度に平行な直線を作成し、直線と正規化後の解析指数である元素強度の深度別変化曲線とが交差する。これにより、「直線右側の点の数量が５個より多い層位、且つ解析指数中で合致する層位をイベント層とする」に従って、イベント層の識別を実現する。
さらに、工程Ａは具体的に以下の方式により実現する。
実験室において得られた堆積物の岩石コアを中心軸に沿って均等に分割し、そのうちの半分の岩石コアの表面を平坦化処理した後、コアスキャナーを用いて走査解析を行い、堆積物の岩石コア中の常量元素と微量元素の強度データを得る。常量元素には、ケイ素Ｓｉ、カルシウムＣａ、アルミニウムＡｌが含まれ、微量元素には、ジルコニウムＺｒ、ルビジウムＲｂ、ストロンチウムＳｒが含まれる。
さらに、工程Ｂにおいて、堆積層の解析指数はＳｉ／Ｔｉ比率及びＣａの元素強度である。
さらに、工程Ｂおいて、堆積層の解析指数はＺｒ／Ｒｂ比率及びＳｒの元素強度である。
さらに、工程Ｃ１において、Ｓｉ／Ｔｉ比率及びＣａの元素強度データの平均値
及び標準偏差δをそれぞれ確定するときには、以下の方式を採用する。
（１）最初にＳｉ／Ｔｉ比率及びＣａの元素強度データの正規化処理を行い、且つ正規化後のＳｉ／Ｔｉ及びＣａ元素の正規化後の正規化標準偏差ＳＤをそれぞれ計算する。
（２）次に箱ひげ図及び四分位数の解析方法により、上述のＳｉ／Ｔｉ比率及びＣａの元素強度データに対してそれぞれ処理を行い、且つＳｉ／Ｔｉ比率及びＣａの元素強度データの四分位点における数値をそれぞれ確定する。
（３）四分位数を上回った外れ値を除去した後、Ｓｉ／Ｔｉ比率及びＣａの元素強度データの平均値
及び標準偏差δをそれぞれ計算する。
さらに、工程Ｃ２において、正規化標準偏差ＳＤを基に、岩石コア堆積層を識別する閾値を確定する。
１）ＳＤ≧０．１の場合、
値を閾値として設定する。
２）ＳＤ＜０．１の場合、
値を閾値として設定する。
閾値とは、堆積物の岩石コアにおける正常堆積層とイベント層を区分する臨界値のことをいう。
さらに、工程Ｃ３において、イベント層の確定は、具体的には以下の方式を採用する。
１）ＳＤ≧０．１の場合、正規化後のＳｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線において、
及び、
の２つの臨界値により深度に平行な直線を作成し、直線と曲線が交差する部分の深度層位が区間内にある層位は正常堆積層と定義し、
より小さく
より大きい区間内にある層位はいずれもイベント層と定義する。
より小さい区間内の層位は適切な堆積学的解釈が得られないため、
より大きい区間内の層位をイベント層と定義する。
２）ＳＤ＜０．１の場合、正規化後のＳｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線中、
及び
の２つの臨界値により深度に平行な垂線を作成し、垂線と曲線が交差する部分の深度層位が区間
内にある層位は正常堆積層と定義し、
より小さく
より大きい区間内にある層位はいずれもイベント層と定義する。
より小さい区間内の層位は適切な堆積学的解釈が得られないため、
より大きい区間内の層位をイベント層と定義する。
３）正規化標準偏差ＳＤの値に基づき、正規化後のＳｉ／Ｔｉ比率及びＣａ元素の深度別変化曲線に対して、選定した閾値により深度に平行な直線を作成し、直線右側に位置する点の数量が５個以上の層位をイベント層とし、且つ直線と正規化後のＳｉ／Ｔｉ比率及びＣａ元素の深度別変化曲線の交差部における曲線中の一つ目の点をイベント層の起点とし、交差部における曲線中の最後の点をイベント層の終点とし、これにより岩石コア中の第１イベント層を得る。第１イベント層以下をこのように類推して、岩石コア中の第ｎイベント層を識別する。
従来技術と比較すると、本発明は以下の優位点及び好ましい効果を有する。
本発明は、無損失且つ高分解能の岩石コア元素走査技術により、得られた海域堆積物の岩石コアに対して走査解析を行い、岩石コア中の元素の地球化学的な強度データを迅速に取得し、イベント層の識別方法と結合することで、岩石コア中のイベント層に対する高効率且つ精確な識別を実現することが可能であり、さらに本案は、岩石コア中の微小なイベント層を識別できるだけでなく、識別効率を大幅に高め、イベント層識別の正確性を大きく向上させており、海洋における古気候及び古環境研究の発展に有力な技術支援を提供し、高い実用価値を具備している。
本発明の上記目的及び利点をより明解にするため、以下で図面に基づき本発明の具体的な実施例について詳細に説明する。
＜実施例＞
＜実施例＞
堆積物の岩石コアにおけるイベント層の識別方法であって、本実施例では海域堆積物の岩石コアを例に紹介するが、図１に示す通り、以下の工程を含む。
工程１、堆積物の岩石コアを分割し、Ｘ線により堆積物の岩石コアの元素に対する走査を行い、堆積物の岩石コアに含まれる元素の強度データを得る。具体的には以下の通りである。
実験室において得られた堆積物の岩石コアを中心軸に沿って均等に分割し、そのうちの半分の岩石コアの表面を平坦化処理した後、コアスキャナーを用いて走査解析を行い、堆積物の岩石コア中の常量元素と微量元素の強度データを得る。常量元素には、ケイ素（Ｓｉ）、カルシウム（Ｃａ）、アルミニウム（Ａｌ）などが含まれ、微量元素には、ジルコニウム（Ｚｒ）、ルビジウム（Ｒｂ）、ストロンチウム（Ｓｒ）などが含まれる。
一般的に、取得する海域堆積物の岩石コア試料は、直径７５ｍｍ又は１１０ｍｍの円柱体であり、長さは通常数十センチメートルから数メートルの間とされる。そのうち、本実施例において採用するコアスキャナーは限定されず、ＣＯＲＴＥＸシリーズのＸ−ｒａｙ Ｃｏｒｅｓｃａｎｎｅｒ又はＩＴＲＡＸシリーズのＸ−ｒａｙ Ｃｏｒｅｓｃａｎｎｅｒでもよい。堆積層中の微小なイベント層を精確に識別できるようにするため、コアスキャナーの走査ステップ幅は１ｍｍに設定する。ステップ幅が小さいほど、精度はより高くなる。即ち、長さ２００ｃｍの岩石コアから２０００列の元素強度データを得ることができる。
工程２、得られた元素の強度データを基に、堆積層の解析指数を確定する。
得られた複数種類の元素強度のデータを基にした実験、モデリング、解析により、且つ堆積層における元素の賦存特性を考慮することにより、Ｓｉ／Ｔｉ及びＣａの元素強度又はその他の元素指数（例えばＺｒ／Ｒｂ及びＳｒの元素強度）を解析指数として選択する。一般的に、Ｘ−ｒａｙ Ｃｏｒｅｓｃａｎｎｅｒによる走査解析で得られるＺｒ、Ｒｂ及びＳｒの元素強度は弱めであり、識別精度を保証するため、本実施例ではＳｉ／Ｔｉ及びＣａの元素強度を堆積物の岩石コアにおけるイベント層の解析指数として選択している。
海洋堆積物中のＳｉは主に堆積物中の粗粒成分中に賦存されており、Ｔｉは主に細粒堆積物中に賦存されているため、Ｓｉ／Ｔｉ比率は堆積層中の粗粒・細粒鉱物含有量の差異を良好に示し得る。Ｃａは主に堆積物中のデトライタス（生物破片）中に存在するため、Ｃａ含有量は比較的粗いデトライタスの堆積物における相対的寄与を示し得る。
工程３、工程２で確定した堆積層の解析指数を基に、その元素の強度データに対して解析を行うことでイベント層の識別が実現するが、具体的には以下を含む。
（１）Ｓｉ／Ｔｉ比率及びＣａの元素強度データの平均値
及び標準偏差δをそれぞれ計算する。
イベント層を識別する前に、最初にＳｉ／Ｔｉ比率及びＣａの元素強度データの正規化処理を行い、且つ正規化後のＳｉ／Ｔｉ比率及びＣａ元素の標準偏差ＳＤをそれぞれ計算する。次に、箱ひげ図及び四分位数の解析方法により、上述のＳｉ／Ｔｉ比率及びＣａの元素強度データに対してそれぞれ処理を行い、当該方法に従ってＳｉ／Ｔｉ比率及びＣａの元素強度データの四分位点における数値をそれぞれ確定する。その後、四分位数を上回った外れ値を除去した後、Ｓｉ／Ｔｉ比率及びＣａの元素強度データの平均値
及び標準偏差δをそれぞれ計算する。また、本実施例中、上述の正規化処理方法、箱ひげ図及び四分位数の解析方法はいずれも当業者公知の技術であり、ここでは詳述しない。
（２）正規化標準偏差ＳＤを基に、岩石コア堆積層を識別する閾値を確定する。
１）ＳＤ≧０．１の場合、
値を閾値として設定する。
２）ＳＤ＜０．１の場合、
値を閾値として設定する。
上述の１）及び２）中、閾値とは、岩石コア中の正常堆積層とイベント層を区分する臨界値のことをいう。
（３）設定した岩石コア堆積層を識別する閾値により、イベント層を確定する。
１）ＳＤ≧０．１の場合、本実施例では正規化後のＳｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線中、閾値により深度に平行な直線を作成し、直線と曲線が交差する部分の深度層位が区間
内にある層位は正常堆積層と定義し、この区間以外の堆積層はイベント層と定義するが、
より小さい区間内の層位は適切な堆積学的解釈が得られないため、本発明では
より大きい区間内の層位をイベント層と定義する。
２）ＳＤ＜０．１の場合、本発明では正規化後のＳｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線中、閾値により深度に平行な垂線を作成し、垂線と曲線が交差する部分の深度層位が区間
内にある層位は正常堆積層と定義し、この区間以外の堆積層はイベント層と定義するが、
より小さい区間内の層位は適切な堆積学的解釈が得られないため、本発明では
より大きい区間内の層位をイベント層と定義する。
３）正規化標準偏差ＳＤの値に基づき、正規化後のＳｉ／Ｔｉ比率及びＣａ元素の深度別変化曲線に対して、適切な閾値を選定して（閾値を確定する原則に従い、ＳＤ≧０．１の場合、
値を閾値として設定し、ＳＤ＜０．１の場合、
値を閾値として設定する）深度に平行な直線を作成し、直線と正規化後のＳｉ／Ｔｉ比率及びＣａ元素の深度別変化曲線とが交差する。直線右側に位置する点の数量が５個以上の層位をイベント層とし、且つ直線右側の交差部における曲線中の一つ目の点をイベント層の起点とし、交差部における曲線中の最後の点をイベント層の終点とし、これにより岩石コア中の第１イベント層を得る。第１イベント層以下において、直線右側に位置する点の数量が５個以上の層位を第２イベント層とし、このように類推して、岩石コア中の第ｎイベント層を識別し、岩石中のＳｉ／Ｔｉ比率及びＣａ元素の曲線中のイベント層をそれぞれ識別することができるが、堆積層をより正確に識別するため、本発明では正規化後のＳｉ／Ｔｉ比率及びＣａ元素の深度別に変化する２本の曲線中、互いに合致するイベント層をイベント層と定義している。
以下は、実際の応用事例を例に本案の優位性について詳述する。
事案１：中国長江の三角州の水中から岩石コアＡを選んで採取した。
工程１：岩石コアの長さは１００ｃｍであり、実験室内で、得られた直径７５ｍｍの堆積物の岩石コアＡをコアカッターで中心軸に沿って２つに分割し、解析準備としてそのうちの半分にトリミング・バフ処理を行った。
工程２：実験室内でコアスキャナーを用いてバフ処理後の半分の岩石コアの走査解析を行い、岩石コア中の元素強度の深度別に変化する強度変化曲線を得た。スキャナーの動作条件は、以下の通り設定した。電圧５０ｋＶ、電流３０ｍＡ、露光時間８００ｍｓ、スキャン幅１ｍｍ。
工程３：工程２で得た岩石コアの元素強度に対して処理を行い、Ｓｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線を得るとともに、箱ひげ図及び四分位数の解析方法により、Ｓｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線の三分位点及び四分位点における数値を確定した。次に、上述の閾値を確定する原則に従い、計算して得られたＳｉ／Ｔｉ比率及びＣａ元素強度のＳＤ値はそれぞれ０．０３と０．０５であったため、
（０．２６）値をＳｉ／Ｔｉ曲線の閾値として選択し、
（０．１８）をＣａ元素強度曲線の閾値として選択した。
工程４：図２を参照して、Ｓｉ／Ｔｉ及びＣａ曲線中で０．２６と０．１８の垂線がそれぞれ作成され、垂線と２本の曲線が交差しており、「垂線右側に位置する点の数量が５個以上の層位をイベント層とし、且つ垂線右側の交差部における曲線中の一つ目の点をイベント層の起点とし、交差部における曲線中の最後の点をイベント層の終点とし、これにより岩石コア中の第１イベント層を得る。第１イベント層以下において、垂線右側に位置する点の数量が５個以上の層位を第２イベント層とし、且つ垂線右側の交差部における曲線中の一つ目の点をイベント層の起点とし、交差部における曲線中の最後の点をイベント層の終点とし、これにより岩石コア中の第２イベント層を得る。このように類推して、岩石コア中の第ｎイベント層を識別することができる」という判定方法に従って、岩石中のＳｉ／Ｔｉ曲線中で１０個のイベント層が識別され、Ｃａ元素曲線中で９個のイベント層が識別され、２本の曲線中で互いに合致したイベント層はそれぞれ[１]、[２]、[３]、[４]、[５]、[６]、[７]、[８]の合計８か所であった。一番右側は堆積学において常套の平均粒径に基づく指数であり、合計８個のイベント層が識別された。
２種類の方法の比較からは、従来の堆積学における識別方法を用いた場合、幾つかの微小なイベント層は識別が得られていないことが分かる（イベント層[７]と[８]など）。平均粒径の結果は手動による縮分で得たものであり、縮分過程中で人為的な誤差が存在し、識別されたイベント層の位置が不正確になっている。また、堆積物粒度の実験過程においても人為的な誤差要因が存在し、識別されたイベント層に誤りが存在するおそれがある。例えば、平均粒径に基づく指数では、イベント層１と２が識別されているが、Ｓｉ／Ｔｉ及びＣａ元素強度の２つの指数中では識別されていない。この点も、本発明の方法の高い効果性と正確性をさらに示すものである。
事例２：中国東シナ海の大陸棚内から岩石コアＢを選んで採取した。
工程１’：岩石コアの長さは１００ｃｍであり、実験室内で、得られた直径７５ｍｍの堆積物の岩石コアＢをコアカッターで中心軸に沿って２つに分割した。解析準備として、そのうちの半分にトリミング・バフ処理を行った。
工程２’：実験室内でコアスキャナーを用いてバフ処理後の半分の岩石コアの走査解析を行い、岩石コア中の元素の深度別に変化する強度変化曲線を得た。スキャナーの動作条件は以下の通り設定した。電圧５０ｋＶ、電流３０ｍＡ、露光時間８００ｍｓ、スキャン幅１ｍｍ。
工程３’：工程２’中で得た岩石コアの元素強度に対して処理を行い、Ｓｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線を得るとともに、箱ひげ図及び四分位数の解析方法により、Ｓｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線の三分位点及び四分位点における数値を確定した。次に、上述の閾値を確定する原則に従い、計算して得られたＳｉ／Ｔｉ比率及びＣａ元素強度のＳＤ値はそれぞれ０．０６と０．０７であったため、
を閾値として選択し、計算を経て０．２６と０．６６をＳｉ／Ｔｉ及びＣａ元素強度曲線の閾値として確定した。
工程４’：図３を参照して、Ｓｉ／Ｔｉ及びＣａ曲線中で０．２６と０．６６の垂線がそれぞれ作成され、垂線と２本の曲線が交差しており、「垂線右側に位置する点の数量が５個以上の層位をイベント層とし、且つ垂線右側の交差部における曲線中の一つ目の点をイベント層の起点とし、交差部における曲線中の最後の点をイベント層の終点とし、これにより岩石コア中の第１イベント層を得る。第１イベント層以下において、垂線右側に位置する点の数量が５個以上の層位を第２イベント層とし、且つ垂線右側の交差部における曲線中の一つ目の点をイベント層の起点とし、交差部における曲線中の最後の点をイベント層の終点とし、これにより岩石コア中の第２イベント層を得る。このように類推して、岩石コア中の第ｎイベント層を識別することができる」という判定方法に従って、岩石中のＳｉ／Ｔｉ曲線中で５個のイベント層が識別され、Ｃａ元素曲線中で５個のイベント層が識別され、２本の曲線中で互いに合致したイベント層はそれぞれ[１]、[２]、[３]、[４]、[５]の合計５か所であった。一番右側は堆積学において常套の平均粒径に基づく指数であり、合計４個のイベント層が識別された。
２種類の方法の比較からは、従来の堆積学における識別方法を用いた場合、幾つかの微小なイベント層は識別が得られていないことが分かる（イベント層[４]など）。平均粒径の結果は手動による縮分で得たものであり、縮分過程で人為的な誤差が存在し、識別されたイベント層の位置が不正確になっている。例えば、２種類の方法が識別したイベント層[２]、[３]及び[４]の位置がやや大きく異なっている。従来の堆積学における識別方法を用いて識別されたイベント層は大きめであり、且つ確定されたイベント層の位置も不正確であるため、後続の解析過程において誤認識をもたらす可能性がある。
上述は本発明の好ましい実施例に過ぎず、本発明が他の形態を取るのを限定するものではなく、当業者であれば上述で開示した技術内容を利用し、変更又は改造を加えて等価変更した同等の実施例を他の分野に応用することができるが、そのいずれも本発明の技術案の内容を逸脱せず、本発明の技術的本質に基づき上述の実施例に対して行う何らかの簡単な修正、等価変更及び改造はいずれも本発明の技術案の保護範囲に属する。
Claims (5)
- 堆積物の岩石コアを分割し、Ｘ線により堆積物の岩石コアの元素に対する走査を行い、堆積物の岩石コアに含まれる元素の強度データを得る工程Ａと、
得られた元素の強度データを基に、堆積層における元素の賦存特性に従って堆積層の解析指数を確定し、前記解析指数はＳｉ／Ｔｉ比率及びＣａの元素強度を含む工程Ｂと、
前記工程Ｂで確定した堆積層の解析指数を基に、前記元素の強度データに対して解析を行うことでイベント層の識別を行う工程Ｃと、を含み、前記工程Ｃは、
（Ｃ１）解析指数の平均値
及び標準偏差δを確定する；
（Ｃ２）解析指数を正規化した後の正規化標準偏差ＳＤを基に、イベント層を識別する閾値を確定し、前記正規化後の標準偏差ＳＤは、解析指数の元素強度データに対して正規化処理を行って得たものであり、前記閾値とは、堆積物の岩石コアにおける正常堆積層とイベント層を区分する臨界値のことをいう；
（Ｃ３）イベント層を確定するが、具体的には以下の方式を採用する：
１）ＳＤ≧０．１の場合、正規化後のＳｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線中、
及び
の２つの臨界値により深度に平行な直線を作成し、前記直線と曲線が交差する部分の深度層位が区間
内にある層位は正常堆積層と定義し、
より大きい区間内の層位をイベント層と定義する；
２）ＳＤ＜０．１の場合、正規化後のＳｉ／Ｔｉ比率及びＣａ元素強度の深度別変化曲線中、
及び
の２つの臨界値により深度に平行な垂線を作成し、前記垂線と曲線が交差する部分の深度層位が区間
内にある層位は正常堆積層と定義し、
より大きい区間内の層位をイベント層と定義する；
３）正規化標準偏差ＳＤの値に基づき、正規化後のＳｉ／Ｔｉ比率及びＣａ元素の深度別変化曲線に対して、選定した閾値により深度に平行な直線を作成し、正規化後の元素強度の変化曲線中のデータ点の値が閾値を上回る位置にある点の数量が５個以上の層位をイベント層とし、且つ前記直線と正規化後のＳｉ／Ｔｉ比率及びＣａ元素の深度別変化曲線の交差部における曲線中の一つ目の点をイベント層の起点とし、交差部における曲線中の最後の点をイベント層の終点とし、これにより岩石コア中の第１イベント層を得て、第１イベント層以下を類推し、岩石コア中の第ｎイベント層を識別する、という工程を含むことを特徴とする、堆積物の岩石コアにおけるイベント層の識別方法。 - 前記工程Ａは、具体的には、
実験室において得られた堆積物の岩石コアを中心軸に沿って均等に分割し、そのうちの半分の岩石コアの表面を平坦化処理した後、コアスキャナーを用いて走査解析を行、堆積物の岩石コア中の常量元素と微量元素の強度データを得るが、前記常量元素には、ケイ素Ｓｉ、カルシウムＣａ、アルミニウムＡｌが含まれ、前記微量元素には、ジルコニウムＺｒ、ルビジウムＲｂ、ストロンチウムＳｒが含まれる、という方式で行われることを特徴とする、請求項１に記載の堆積物の岩石コアにおけるイベント層の識別方法。 - 前記工程Ｂにおいて、前記堆積層の解析指数はＺｒ／Ｒｂ比率及びＳｒの元素強度であることを特徴とする、請求項１又は請求項２に記載の堆積物の岩石コアにおけるイベント層の識別方法。
- 前記工程Ｃ１において、Ｓｉ／Ｔｉ比率及びＣａの元素強度データの平均値
及び標準偏差δをそれぞれ確定するときには、
（１）最初にＳｉ／Ｔｉ比率及びＣａの元素強度データの正規化処理を行い、且つ正規化後のＳｉ／Ｔｉ及びＣａ元素の正規化後の正規化標準偏差ＳＤをそれぞれ計算する；
（２）次に、箱ひげ図及び四分位数の解析方法により、前記Ｓｉ／Ｔｉ比率及びＣａの元素強度データに対してそれぞれ処理を行い、且つＳｉ／Ｔｉ比率及びＣａの元素強度データの四分位点における数値をそれぞれ確定する；
（３）四分位数を上回った外れ値を除去した後、Ｓｉ／Ｔｉ比率及びＣａの元素強度データの平均値
及び標準偏差δをそれぞれ計算する、という方式を採用することを特徴とする、請求項１に記載の堆積物の岩石コアにおけるイベント層の識別方法。 - 前記工程Ｃ２中、正規化標準偏差ＳＤを基に岩石コア堆積層を識別する閾値を確定するときには、
１）ＳＤ≧０．１の場合、
値を閾値として設定し、
２）ＳＤ＜０．１の場合、
値を閾値として設定し、
前記閾値とは、堆積物の岩石コアにおける正常堆積層とイベント層を区分する臨界値をいうことを特徴とする、請求項４に記載の堆積物の岩石コアにおけるイベント層の識別方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910923092.3A CN110726740B (zh) | 2019-09-27 | 2019-09-27 | 一种沉积物岩心中事件沉积层的识别方法 |
| CN201910923092.3 | 2019-09-27 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2021056212A JP2021056212A (ja) | 2021-04-08 |
| JP6872762B2 true JP6872762B2 (ja) | 2021-05-19 |
Family
ID=69218495
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020142431A Active JP6872762B2 (ja) | 2019-09-27 | 2020-08-26 | 堆積物の岩石コアにおけるイベント層の識別方法 |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP6872762B2 (ja) |
| CN (1) | CN110726740B (ja) |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106204591B (zh) * | 2016-07-11 | 2018-12-18 | 青岛海洋地质研究所 | 一种海域沉积物岩心中古环境信息载体的取样方法 |
| CN106885815B (zh) * | 2017-03-22 | 2020-01-31 | 中国石油大学(华东) | 一种基于砂岩中元素xrf法判别沉积环境的方法 |
-
2019
- 2019-09-27 CN CN201910923092.3A patent/CN110726740B/zh active Active
-
2020
- 2020-08-26 JP JP2020142431A patent/JP6872762B2/ja active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110726740A (zh) | 2020-01-24 |
| CN110726740B (zh) | 2020-06-30 |
| JP2021056212A (ja) | 2021-04-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ojala et al. | Characteristics of sedimentary varve chronologies–a review | |
| CN103206208B (zh) | 一种微观不同赋存状态剩余油宏观量化的方法 | |
| KR101707440B1 (ko) | 가도리니움 화합물을 반응시킨 암석 시료의 전자현미경 이미지를 활용한 공극률 측정방법 | |
| GB2524810A (en) | Method of analysing a drill core sample | |
| Kazak et al. | Machine-learning-assisted segmentation of focused ion beam-scanning electron microscopy images with artifacts for improved void-space characterization of tight reservoir rocks | |
| Reedy et al. | High-resolution micro-CT with 3D image analysis for porosity characterization of historic bricks | |
| CN1198217A (zh) | 自动确定地层中层面的方法 | |
| CN109580687A (zh) | 一种识别砂岩型铀成矿目的层物源的综合方法 | |
| CN115113299B (zh) | 智能标志地层定位划分方法、装置、设备及存储介质 | |
| CN111060922A (zh) | 基于机载激光雷达点云空间分布特征的树木点云提取方法 | |
| JP6872762B2 (ja) | 堆積物の岩石コアにおけるイベント層の識別方法 | |
| Kolendrianou et al. | The Palaeolithic cave of Kalamakia (Mani Peninsula), Greece: new insights on the palaeoenvironment using microvertebrates and mesowear analysis of ruminant teeth | |
| Hale et al. | Elemental signatures in the vertebral cartilage of the round stingray, Urobatis halleri, from Seal Beach, California | |
| CN109709623B (zh) | 非破坏性的金矿化作用探测方法 | |
| Cecil et al. | Feasibility of field-portable XRF to identify obsidian sources in Central Peten, Guatemala | |
| US20230154208A1 (en) | Method of Detecting at Least One Geological Constituent of a Rock Sample | |
| CN108805158B (zh) | 一种致密油储层成岩相划分方法 | |
| Choi et al. | Microstructural and crystallographic evolution of palaeognath (Aves) eggshells | |
| CN113075731B (zh) | 深层储层连续性井筒数字建模方法及装置 | |
| CN112950551A (zh) | 一种海洋沉积物岩心年代快速标定方法 | |
| CN109444189B (zh) | 利用数字岩石分析技术开展复杂地层对比与定量评价方法 | |
| Tripathi et al. | Exporting petrophysical properties of sandstones from thin section image analysis | |
| CN114428063B (zh) | 污染层光谱数据识别方法、装置、电子设备和存储介质 | |
| Petersen et al. | Provenance of Jurassic sedimentary rocks of south-central Quesnellia, British Columbia: implications for paleogeography | |
| CN112612935B (zh) | 一种基于自推理模型的完整测井数据获取方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200826 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20200826 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201015 |
|
| A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20201116 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20201201 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210225 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20210316 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210412 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6872762 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |