JP2022013627A - 岩石破砕によって展開する岩石の割れ目ネットワークのトポロジー構造を描写する方法 - Google Patents
岩石破砕によって展開する岩石の割れ目ネットワークのトポロジー構造を描写する方法 Download PDFInfo
- Publication number
- JP2022013627A JP2022013627A JP2021016639A JP2021016639A JP2022013627A JP 2022013627 A JP2022013627 A JP 2022013627A JP 2021016639 A JP2021016639 A JP 2021016639A JP 2021016639 A JP2021016639 A JP 2021016639A JP 2022013627 A JP2022013627 A JP 2022013627A
- Authority
- JP
- Japan
- Prior art keywords
- crack
- rock
- steps
- network
- describing
- 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.)
- Granted
Links
- 239000011435 rock Substances 0.000 title claims abstract description 72
- 238000000034 method Methods 0.000 title claims abstract description 47
- 230000008569 process Effects 0.000 claims abstract description 17
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000006073 displacement reaction Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 13
- 230000008859 change Effects 0.000 claims description 11
- 238000011065 in-situ storage Methods 0.000 claims description 10
- 230000035945 sensitivity Effects 0.000 claims description 9
- 230000009466 transformation Effects 0.000 claims description 7
- 238000011161 development Methods 0.000 abstract description 6
- 238000004422 calculation algorithm Methods 0.000 description 10
- 238000001514 detection method Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0062—Crack or flaws
- G01N2203/0066—Propagation of crack
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/006—Crack, flaws, fracture or rupture
- G01N2203/0067—Fracture or rupture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/0641—Indicating or recording means; Sensing means using optical, X-ray, ultraviolet, infrared or similar detectors
- G01N2203/0647—Image analysis
Abstract
【課題】岩石破砕割れ目ネットワークのトポロジー構造を描写する方法を提供する。【解決手段】一実施例による描写方法は、画像処理装置を用いて変形時刻の異なる岩石の原位置画像を取得するステップと、各画像に基づいて、原位置のサブピクセルクラックの空間的広がりを取得するステップと、を含む。この解決手段によれば、構造離散化の形態の違いによって結果が異なることが避けられるために、画像境界条件を簡略化する必要がなく、かつ固有の割れ目情報を含まないので、岩石割れ目ネットワークを効果的に描写し、岩石割れ目ネットワークの展開の全過程を得ることができる。【選択図】図１
Description
本発明は、岩石破砕検出の技術分野に関し、具体的には、岩石破砕割れ目ネットワークのトポロジー構造を描写する方法に関する。
岩石破砕は、岩盤工学では一般的な現象であり、岩石破砕の進化法則の研究、および岩石破砕時間やその位置の予測は、基礎的な工学の構築や天然資源の効率的な開発、深部地球における断層活性化の理解と評価に不可欠である。岩石の変形・破壊時に非連続的な超低周波音が発生し、掘削中の周囲岩石の騒音をモニタニングした結果、岩石のバーストには顕著な時間的・空間的特性があることが分かる。
岩石破砕の時間的・空間的影響を研究するための有効な手段の一つとして、岩石割れ目ネットワークの展開が挙げられる。岩石割れ目ネットワークの描写は、主に数値シミュレーション、割れ目の幾何学的特徴法（長さ、幅、形状などを含む）、統計的手法（形状、長さなどの異なる割れ目を含む）が用いられている。
数値シミュレーション手法では、岩石割れ目ネットワークの展開の全過程を得ることができるが、数値シミュレーションの過程では境界条件を簡略化する必要があり、また、構造離散化の形態によって得られる結果が異なる。割れ目の幾何学的特徴法や統計的手法では、固有の割れ目情報しか含まれておらず、割れ目の膨張傾向を記述することができない。
どのように効果的に岩石割れ目ネットワークを描写するかは、現在、当業者が解決すべき問題である。
上記技術的課題を解決するために、または少なくとも部分的に上記技術的課題を解決するために、本発明は、岩石破砕割れ目ネットワークのトポロジー構造を描写する方法を提供する。
岩石破砕割れ目ネットワークのトポロジー構造を描写する方法は、
画像処理装置を用いて、変形時刻の異なる岩石の原位置画像を取得するステップと、
前記各画像に基づいて、原位置のサブピクセルクラックの空間的広がりを取得するステップと、を含む。
画像処理装置を用いて、変形時刻の異なる岩石の原位置画像を取得するステップと、
前記各画像に基づいて、原位置のサブピクセルクラックの空間的広がりを取得するステップと、を含む。
選択的に、上述した原位置のサブピクセルクラックの空間的広がりを取得することは、
岩石原位置の変位場を構築するステップと、
サブピクセルのクラック分布を描写するステップと、
割れ目ネットワークのトポロジー構造を構築するステップと、を含む。
岩石原位置の変位場を構築するステップと、
サブピクセルのクラック分布を描写するステップと、
割れ目ネットワークのトポロジー構造を構築するステップと、を含む。
選択的に、上述した岩石原位置の変位場を構築するステップでは、２つの前記画像の類似度を測定するデジタルボリューム計算処理を含む。
選択的に、前記デジタルボリューム計算処理は、
参照画像と変形画像を入力するステップと、
初期変形関数を設定するステップと、
割れ目分布サブセットと検索サブセットを分割するステップと、
各サブセットの変更行列を計算するステップと、
各サブセットの相関係数を計算するステップと、
相関係数の最大値と座標を見つけるステップと、
各サブセットの変位、角度を計算するステップと、
エラー要件を満たしているか否かを判定するステップと、
変形関数を求めるステップと、を含む。
参照画像と変形画像を入力するステップと、
初期変形関数を設定するステップと、
割れ目分布サブセットと検索サブセットを分割するステップと、
各サブセットの変更行列を計算するステップと、
各サブセットの相関係数を計算するステップと、
相関係数の最大値と座標を見つけるステップと、
各サブセットの変位、角度を計算するステップと、
エラー要件を満たしているか否かを判定するステップと、
変形関数を求めるステップと、を含む。
選択的に、上述したエラー要件を満たしているか否かを判定するステップでは、
エラー要件を満たしていない場合、各サブセットの変更行列を再計算し、上述した各サブセットの変更行列を計算するステップから、他のステップを順次行い、
エラー要件を満たしている場合、変形関数を求める。
エラー要件を満たしていない場合、各サブセットの変更行列を再計算し、上述した各サブセットの変更行列を計算するステップから、他のステップを順次行い、
エラー要件を満たしている場合、変形関数を求める。
選択的に、上述したサブピクセルのクラック分布を描写するステップでは、空間のスカラーを取得することが含まれ、具体的には、
クラック面がｘ軸のクラックに垂直になり、ｘ軸に沿った変位を関数ｕとする場合、クラック位置では、
を確立し、ここで、Ｋは空間のスカラーであり、Ｕは変位ベクトル場であり、前記スカラーを利用して、各位置で発生する破砕の傾向を評価する。
クラック面がｘ軸のクラックに垂直になり、ｘ軸に沿った変位を関数ｕとする場合、クラック位置では、
選択的に、上述した割れ目ネットワークのトポロジー構造を構築することは、
岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立するステップと、
異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、前記割れ目ネットワークに基づく隣接行列を確立するステップと、
割れ目頂点の感度を記述するステップと、を含む。
岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立するステップと、
異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、前記割れ目ネットワークに基づく隣接行列を確立するステップと、
割れ目頂点の感度を記述するステップと、を含む。
選択的に、上述した岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立する方法は、
であり、ここで、Ｇは図形であり、Ｖは前記岩石割れ目の頂点であり、Ｅは前記各頂点間に形成された辺であり、
前記接続関係では、前記割れ目のぞれぞれは、前記辺に対応し、
上述した割れ目頂点の感度は、次の段階で発生する破砕の確率である。
前記接続関係では、前記割れ目のぞれぞれは、前記辺に対応し、
上述した割れ目頂点の感度は、次の段階で発生する破砕の確率である。
上述した異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、前記割れ目ネットワークに基づく隣接行列を確立する方法は、
選択的に、前記割れ目ネットワークでは、前記割れ目の展開を記述し、前記割れ目の２つの頂点が接続可能であるか否かを決定し、前記割れ目の２つの頂点確率を描写する方法は、
であり、ここで、Ｃは割れ目ネットワークノードｉのクラスタリング係数であり、Ｅは頂点ｉを通過する辺の数であり、ｋはノードｉの感度であり、
特定の頂点ｉを通過する割れ目が１つしかない場合、ｋｉ＝１、
特定の頂点ｉを通過する割れ目がない場合、ｋｉ＝０、Ｅｉ＝０、Ｃｉ＝０、
クラスタリングシステムは、岩石荷重の違いによって異なり、前記クラスタリングシステムは、応力－クラスタリング係数曲線を確立し、割れ目の一点接続性を描写し、割れ目ノード度Ｋのクラスタリング係数を利用して岩石割れ目全体の接続性を定量的に描写するように構成され、割れ目クラスタリング係数を計算する方法は、
特定の頂点ｉを通過する割れ目が１つしかない場合、ｋｉ＝１、
特定の頂点ｉを通過する割れ目がない場合、ｋｉ＝０、Ｅｉ＝０、Ｃｉ＝０、
クラスタリングシステムは、岩石荷重の違いによって異なり、前記クラスタリングシステムは、応力－クラスタリング係数曲線を確立し、割れ目の一点接続性を描写し、割れ目ノード度Ｋのクラスタリング係数を利用して岩石割れ目全体の接続性を定量的に描写するように構成され、割れ目クラスタリング係数を計算する方法は、
本発明では、構造離散化の形態の違いによって結果が異なることを避けるために、画像境界条件を簡略化する必要がなく、かつ固有の割れ目情報を含まないので、岩石割れ目ネットワークを効果的に描写し、岩石割れ目ネットワークの展開の全過程を得ることができる。
本発明の上記目的、特徴および利点をより明らかにするために、以下、図面および実施例を組み合わせて本発明のいくつかの実施例をより詳細に説明する。説明される実施例は、本発明の実施例の一部に過ぎず、その全てではないことを理解すべきである。以下の実施例は、本発明を解釈するためのものに過ぎず、本発明を限定するためのものではない。説明される本発明の実施例に基づいて、当業者によって得られる他の全ての実施例は、本発明の保護の範囲内に属する。特に指定されていない場合、実施例で使用される技術的解決手段は、当業者が熟知している従来の手段または先行技術である。
本文では、「第１」や「第２」などの関係用語は、あるエンティティや操作を他のエンティティや操作と区別するためにのみ使用され、これらのエンティティや操作の間の実際の関係や順序を必ずしも必要としたり、暗示したりするものではないことに留意されたい。「接続」、「連結」などの用語は、広義に理解されるべきであり、例えば、固定接続されてもよく、着脱可能に接続または一体に接続されてもよいが、機械的に接続されてもよく、電気的に接続されてもよいが、直接的に接続されてもよく、中間媒体を介して間接的に接続されてもよい。「含む」、「包含」または他の任意の変形は、非排他的な包含をカバーすることを意図しており、一連の要素を含む過程、方法、物品や装置は、それらの要素だけでなく、明示的に列挙されていない他の要素、またはこの過程、方法、物品や装置に固有の要素をも含む。さらなる制限がない限り、「……を含む」という記述によって限定される要素は、前記要素を含む過程、方法、物品や装置に追加した他の同一の要素を排除しない。
本発明によって提供される一実施例は、
岩石破砕割れ目ネットワークのトポロジー構造を描写する方法であって、
画像処理装置を用いて、変形時刻の異なる岩石の原位置画像を取得するステップと、
前記原位置画像の各々に基づいて、原位置のサブピクセルクラックの空間的広がりを取得するステップと、を含む。
岩石破砕割れ目ネットワークのトポロジー構造を描写する方法であって、
画像処理装置を用いて、変形時刻の異なる岩石の原位置画像を取得するステップと、
前記原位置画像の各々に基づいて、原位置のサブピクセルクラックの空間的広がりを取得するステップと、を含む。
本実施例では、構造離散化の形態の違いによって結果が異なることを避けるために、画像境界条件を簡略化する必要がなく、かつ固有の割れ目情報を含まないので、岩石割れ目ネットワークを効果的に描写し、岩石割れ目ネットワーク展開の全過程を得ることができる。
画像処理装置を用いて、変形時刻の異なる岩石原位置画像を取得する。その場合、画像処理装置は、例えば高エネルギー加速器ＣＴ試験機である。ＣＴ画像は、離散的な岩石グレースケール分布を記録しており、どのような統計的手法で計算しても、その計算結果は離散的な画素を単位とし、岩石破砕変形および微小割れ目の展開も画素でしか調べることができず、岩石変形および破砕に重要なサブピクセルの情報を失うことになる。そこで、補間法を用いてＣＴ画像に対して最初の前処理を行い、ほぼ連続した岩石グレースケール分布を得て、連続した画像を再構成する。また、最大エントロピー検出アルゴリズムとクラスタリング画像アルゴリズムを用いて、反転した連続に対して二次処理を行う。
最大エントロピー検出アルゴリズムは、グレースケール分布が周辺環境と大きく異なる領域を強調するのに適しており、計算効率が高く、高速であるという利点を有する。しかし、最大エントロピー検出アルゴリズムは騒音の影響を受けやすいため、エッジおよび／または形状がぼやけているＣＴ画像に対しては、より良好な認識効果が得られるクラスタリング画像アルゴリズムを用いて二次処理を行う。最大エントロピー検出アルゴリズムとクラスタリング画像アルゴリズムに基づいて、ＣＴの画像解像度を向上させることができる。なお、最大エントロピー検出アルゴリズムとクラスタリング画像アルゴリズムは、先行技術におけるアルゴリズムであるため、詳細な説明を省略する。
本発明によって提供される他の実施例では、
原位置のサブピクセルクラックの空間的広がりを取得することは、
岩石原位置の変位場を構築するステップと、
サブピクセルのクラック分布を描写するステップと、
割れ目ネットワークのトポロジー構造を構築するステップと、を含む。
原位置のサブピクセルクラックの空間的広がりを取得することは、
岩石原位置の変位場を構築するステップと、
サブピクセルのクラック分布を描写するステップと、
割れ目ネットワークのトポロジー構造を構築するステップと、を含む。
本発明によって提供される他の実施例では、岩石原位置の変位場を構築するステップで、２つの画像の類似度を測定するデジタルボリューム計算処理を含む。
本実施例では、異なる段階で取得したデータに基づいて、変位場とひずみ場の変化パラメータを判定する。
例えば図１に示すように、本発明によって提供される他の実施例おいて、デジタルボリューム計算処理は、
参照画像と変形画像を入力するステップと、
初期変形関数を設定するステップと、
割れ目分布サブセットと検索サブセットを分割するステップと、
各サブセットの変更行列を計算するステップと、
各サブセットの相関係数を計算するステップと、
相関係数の最大値と座標を見つけるステップと、
各サブセットの変位、角度を計算するステップと、
エラー要件を満たしているか否かを判定するステップと、
変形関数を求めるステップと、を含む。
参照画像と変形画像を入力するステップと、
初期変形関数を設定するステップと、
割れ目分布サブセットと検索サブセットを分割するステップと、
各サブセットの変更行列を計算するステップと、
各サブセットの相関係数を計算するステップと、
相関係数の最大値と座標を見つけるステップと、
各サブセットの変位、角度を計算するステップと、
エラー要件を満たしているか否かを判定するステップと、
変形関数を求めるステップと、を含む。
本実施例のデジタルボリューム計算処理では、計算結果が安定しており、類似度の判定が困難になったり、複数の類似点の発生により相関係数の最適値の判定が困難になったりすることが回避され、さらに、このデジタルボリューム計算処理では、干渉に強く、計算結果を信頼できなくなるような騒音を多く発生させることが回避される。
デジタルボリューム計算処理の各詳細なステップは、先行技術における計算ステップを使用する。
本発明によって提供される他の実施例において、エラー要件を満たしているか否かを判定するステップでは、
エラー要件を満たしていない場合、各サブセットの変更行列を再計算し、各サブセットの変更行列を計算するステップから、他のステップを順次行い、
エラー要件を満たしている場合、変形関数を求める。
エラー要件を満たしていない場合、各サブセットの変更行列を再計算し、各サブセットの変更行列を計算するステップから、他のステップを順次行い、
エラー要件を満たしている場合、変形関数を求める。
本実施例では、さらに、より安定した計算結果を実現し、干渉に対する耐性を向上させる。
本発明によって提供される他の実施例において、サブピクセルのクラック分布を描写するステップでは、空間のスカラーを取得することが含まれ、具体的には、
クラック面がｘ軸のクラックに垂直になり、ｘ軸に沿った変位を関数ｕとする場合、クラック位置では、
を確立し、ここで、Ｋは空間のスカラーであり、Ｕは変位ベクトル場であり、スカラーを利用して、各位置で発生する破砕の傾向を評価する。
クラック面がｘ軸のクラックに垂直になり、ｘ軸に沿った変位を関数ｕとする場合、クラック位置では、
本実施例では、取得したスカラーに基づいて、各位置で発生する破砕の傾向を評価する。
本発明によって提供される他の実施例において、割れ目ネットワークのトポロジー構造を構築するステップでは、
岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立するステップと、
異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、割れ目ネットワークに基づく隣接行列を確立するステップと、
割れ目頂点の感度を記述するステップと、を含む。
岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立するステップと、
異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、割れ目ネットワークに基づく隣接行列を確立するステップと、
割れ目頂点の感度を記述するステップと、を含む。
本実施例では、岩石割れ目の頂点の発達状況を定量的に特徴付けることができる。
本発明によって提供される他の実施例において、岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立する方法は、
であり、ここで、Ｇは図形であり、Ｖは岩石割れ目の頂点であり、Ｅは各頂点の間に形成された辺であり、
接続関係では、割れ目のそれぞれは辺に対応し、
割れ目頂点の感度は、次の段階で発生する破砕の確率である。
接続関係では、割れ目のそれぞれは辺に対応し、
割れ目頂点の感度は、次の段階で発生する破砕の確率である。
本発明によって提供される他の実施例では、異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、割れ目ネットワークに基づく隣接行列を確立する方法は、
本実施例では、岩石微割れ目の展開を定量的に特徴付けることができる。
本発明によって提供される他の実施例では、割れ目ネットワークにおいて、割れ目の展開を記述し、割れ目の２つの頂点が接続可能であるか否かを決定し、割れ目の２つの頂点確率を描写する方法は、
であり、ここで、Ｃは割れ目ネットワークノードｉのクラスタリング係数であり、Ｅは頂点ｉを通過する辺の数であり、ｋはノードｉの感度であり、
特定の頂点ｉを通過する割れ目が１つしかない場合、ｋｉ＝１、
特定の頂点ｉを通過する割れ目がない場合、ｋｉ＝０、Ｅｉ＝０、Ｃｉ＝０、
クラスタリングシステムは、岩石荷重の違いによって異なり、応力－クラスタリング係数曲線を確立し、割れ目の一点接続性を描写し、割れ目ノード度Ｋのクラスタリング係数を利用して岩石割れ目全体の接続性を定量的に描写するように構成され、割れ目クラスタリング係数を計算する方法は、
Ｃ（ｋ）が圧力に従って変化した曲線の傾きの最大値は岩石の破砕時間である。
特定の頂点ｉを通過する割れ目が１つしかない場合、ｋｉ＝１、
特定の頂点ｉを通過する割れ目がない場合、ｋｉ＝０、Ｅｉ＝０、Ｃｉ＝０、
クラスタリングシステムは、岩石荷重の違いによって異なり、応力－クラスタリング係数曲線を確立し、割れ目の一点接続性を描写し、割れ目ノード度Ｋのクラスタリング係数を利用して岩石割れ目全体の接続性を定量的に描写するように構成され、割れ目クラスタリング係数を計算する方法は、
本実施例では、岩石割れ目展開のＣ（ｋ）が力率法の形態を有し、岩石力学における応力－ひずみ曲線と１対１の対応関係を有する応力－Ｃ（Ｋ）の変化曲線を確立することができ、応力－ひずみ曲線に応答して岩石の内部破砕メカニズムを明らかにするために使用される。
最後に、上記の実施例は、本発明の技術的解決手段を例示するためにのみ使用され、それらを限定するものではないことに留意すべきである。上記実施例を参照して本発明を詳細に説明したが、当業者であれば、上記の各実施例に記載された技術的解決手段を変更し、またはそれらの技術的特徴の一部を等価的に置き換えることができ、本発明の精神から逸脱しないことに基づいて行われた変更や置き換えは、本発明の保護の範囲に含まれる。
Claims (10)
- 画像処理装置を用いて、変形時刻の異なる岩石の原位置画像を取得するステップと、
各前記画像に基づいて、原位置のサブピクセルクラックの空間的広がりを取得するステップと、を含むことを特徴とする岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。 - 上述した原位置のサブピクセルクラックの空間的広がりを取得することは、
岩石原位置の変位場を構築するステップと、
サブピクセルのクラック分布を描写するステップと、
割れ目ネットワークのトポロジー構造を構築するステップと、を含むことを特徴とする請求項１に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。 - 上述した岩石原位置の変位場を構築するステップでは、２つの前記画像の類似度を測定するデジタルボリューム計算処理を含むことを特徴とする請求項２に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。
- 前記デジタルボリューム計算処理は、
参照画像と変形画像を入力するステップと、
初期変形関数を設定するステップと、
割れ目分布サブセットと検索サブセットを分割するステップと、
各サブセットの変更行列を計算するステップと、
各サブセットの相関係数を計算するステップと、
相関係数の最大値と座標を見つけるステップと、
各サブセットの変位、角度を計算するステップと、
エラー要件を満たしているか否かを判定するステップと、
変形関数を求めるステップと、を含むことを特徴とする請求項３に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。 - 上述したエラー要件を満たしているか否かを判定するステップでは、
エラー要件を満たしていない場合、各サブセットの変更行列を再計算し、上述した各サブセットの変更行列を計算するステップから、他のステップを順次行い、
エラー要件を満たしている場合、変形関数を求めることを特徴とする請求項４に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。 - 上述したサブピクセルのクラック分布を描写するステップでは、空間のスカラーを取得することが含まれ、具体的には、
クラック面がｘ軸のクラックに垂直になり、ｘ軸に沿った変位を関数ｕとする場合、クラック位置では、
を確立し、ここで、Ｋは空間のスカラーであり、Ｕは変位ベクトル場であり、前記スカラーを利用して、各位置で発生する破砕の傾向を評価することを特徴とする請求項２に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。
- 上述した割れ目ネットワークのトポロジー構造を構築することは、
岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立するステップと、
異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、前記割れ目ネットワークに基づく隣接行列を確立するステップと、
割れ目頂点の感度を記述するステップと、を含むことを特徴とする請求項２に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。 - 上述した岩石割れ目の頂点と接続関係に基づいて異なる段階の割れ目ネットワークを確立する方法は、
であり、ここで、Ｇは図形であり、Ｖは前記岩石割れ目の頂点であり、Ｅは前記各頂点間に形成された辺であり、
前記接続関係では、前記割れ目のそれぞれは、前記辺に対応し、
前記割れ目頂点の感度は、次の段階で発生する破砕の確率であることを特徴とする請求項７に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。 - 前記異なる破砕段階の岩石中の割れ目を異なる割れ目ネットワーク構造として確立し、前記割れ目ネットワークに基づく隣接行列を確立する方法は、
- 前記割れ目ネットワークでは、前記割れ目の展開を記述し、前記割れ目の２つの頂点が接続可能であるか否かを決定し、前記割れ目の２つの頂点確率を描写する方法は、
であり、ここで、Ｃは割れ目ネットワークノードｉのクラスタリング係数であり、Ｅは頂点ｉを通過する辺の数であり、ｋはノードｉの感度であり、
特定の頂点ｉを通過する割れ目が１つしかない場合、ｋｉ＝１、
特定の頂点ｉを通過する割れ目がない場合、ｋｉ＝０、Ｅｉ＝０、Ｃｉ＝０、
クラスタリングシステムは、岩石荷重の違いによって異なり、前記クラスタリングシステムは、応力－クラスタリング係数曲線を確立し、割れ目の一点接続性を描写し、割れ目ノード度Ｋのクラスタリング係数を利用して岩石割れ目全体の接続性を定量的に描写するように構成され、割れ目クラスタリング係数を計算する方法は、
Ｃ（ｋ）が圧力に従って変化した曲線の傾きの最大値は岩石の破砕時間であることを特徴とする請求項９に記載の岩石破砕割れ目ネットワークのトポロジー構造を描写する方法。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010617168.2 | 2020-06-30 | ||
| CN202010617168.2A CN111751211B (zh) | 2020-06-30 | 2020-06-30 | 一种岩石压裂缝网拓扑结构刻画方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP6926352B1 JP6926352B1 (ja) | 2021-08-25 |
| JP2022013627A true JP2022013627A (ja) | 2022-01-18 |
Family
ID=72678314
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2021016639A Active JP6926352B1 (ja) | 2020-06-30 | 2021-02-04 | 岩石破砕によって展開する岩石の割れ目ネットワークのトポロジー構造を描写する方法 |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP6926352B1 (ja) |
| CN (1) | CN111751211B (ja) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113899746A (zh) * | 2021-09-30 | 2022-01-07 | 江苏纹动测控科技有限公司 | 一种基于dic的钢结构疲劳裂纹扩展形态测量方法 |
| CN115220100B (zh) * | 2022-06-15 | 2023-09-01 | 西安恒泰艾普能源发展有限公司 | 一种碳酸盐岩水窜通道的分析方法及系统 |
| CN116665079B (zh) * | 2023-06-21 | 2024-02-23 | 南京师范大学 | 复杂场景下顾及空间关系的岩体结构特征识别方法与系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003268778A (ja) * | 2002-03-15 | 2003-09-25 | Railway Technical Res Inst | 剥落型落石および岩石崩壊を対象とした斜面の安定性評価方法 |
| JP2016532867A (ja) * | 2013-08-06 | 2016-10-20 | ビーピー・コーポレーション・ノース・アメリカ・インコーポレーテッド | 疑似応力およびひずみ条件の下での岩石物理特性の画像ベース直接数値シミュレーション |
| WO2019009007A1 (ja) * | 2017-07-06 | 2019-01-10 | 富士フイルム株式会社 | 機械学習データ生成方法、機械学習データ生成プログラム、機械学習データ生成システム、サーバ装置、画像処理方法、画像処理プログラム、及び画像処理システム |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4292095B2 (ja) * | 2004-02-26 | 2009-07-08 | 計測検査株式会社 | 微細ひび割れ幅検量方法 |
| CN101256156B (zh) * | 2008-04-09 | 2011-06-08 | 西安电子科技大学 | 平板裂缝天线裂缝精密测量方法 |
| CN101957178B (zh) * | 2009-07-17 | 2012-05-23 | 上海同岩土木工程科技有限公司 | 一种隧道衬砌裂缝测量方法及其测量装置 |
| JP5645730B2 (ja) * | 2011-03-25 | 2014-12-24 | 公益財団法人鉄道総合技術研究所 | コンクリート表面の閉合ひび割れ検出方法 |
| CN107256557B (zh) * | 2017-05-03 | 2020-05-22 | 华南理工大学 | 一种误差可控的细分曲面图像矢量化方法 |
| CN110763697B (zh) * | 2019-08-26 | 2020-11-06 | 清华大学 | 一种使用飞行器检测工程结构表面裂缝的方法 |
-
2020
- 2020-06-30 CN CN202010617168.2A patent/CN111751211B/zh active Active
-
2021
- 2021-02-04 JP JP2021016639A patent/JP6926352B1/ja active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003268778A (ja) * | 2002-03-15 | 2003-09-25 | Railway Technical Res Inst | 剥落型落石および岩石崩壊を対象とした斜面の安定性評価方法 |
| JP2016532867A (ja) * | 2013-08-06 | 2016-10-20 | ビーピー・コーポレーション・ノース・アメリカ・インコーポレーテッド | 疑似応力およびひずみ条件の下での岩石物理特性の画像ベース直接数値シミュレーション |
| WO2019009007A1 (ja) * | 2017-07-06 | 2019-01-10 | 富士フイルム株式会社 | 機械学習データ生成方法、機械学習データ生成プログラム、機械学習データ生成システム、サーバ装置、画像処理方法、画像処理プログラム、及び画像処理システム |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111751211A (zh) | 2020-10-09 |
| JP6926352B1 (ja) | 2021-08-25 |
| CN111751211B (zh) | 2021-05-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6926352B1 (ja) | 岩石破砕によって展開する岩石の割れ目ネットワークのトポロジー構造を描写する方法 | |
| Truong‐Hong et al. | Combining an angle criterion with voxelization and the flying voxel method in reconstructing building models from LiDAR data | |
| Nishikawa et al. | Concrete crack detection by multiple sequential image filtering | |
| US7933749B2 (en) | System and method for computing minimum distances between two point clouds | |
| Li | Markov chain random fields for estimation of categorical variables | |
| US20070291031A1 (en) | Three dimensional geometric data correction | |
| US20020012451A1 (en) | Method for target detection and identification by using proximity pixel information | |
| US9959670B2 (en) | Method for rendering terrain | |
| Gain et al. | Preventing self-intersection under free-form deformation | |
| JP6380051B2 (ja) | 有限要素演算プログラム、有限要素演算装置および有限要素演算方法 | |
| CN115131560A (zh) | 基于全局特征学习和局部特征判别聚合的点云分割方法 | |
| Kovitz et al. | Spatial statistics of clustered data | |
| CN109947597A (zh) | 一种网络流量数据恢复方法及系统 | |
| CN110793438B (zh) | 基于模糊熵和人工鱼群算法的低速冲击位置辨识方法 | |
| Silva et al. | Multiview range image registration using the surface interpenetration measure | |
| CN111754500B (zh) | 一种岩石压裂缝网拓扑结构刻画系统 | |
| Krysl et al. | An efficient linear‐precision partition of unity basis for unstructured meshless methods | |
| JP2004295255A (ja) | 干渉判定システム及び干渉判定プログラム | |
| Suárez et al. | Four-triangles adaptive algorithms for RTIN terrain meshes | |
| Rahman et al. | Protein similarity analysis by wavelet decomposition of cellular automata images | |
| JP2011221930A (ja) | 三次元地形データ高精度化装置、三次元地形データ高精度化方法及びプログラム | |
| CN110910499B (zh) | 基于Revit软件的地质环境载体断层模型的构建方法及装置 | |
| CN116305461B (zh) | 结构响应计算方法、装置、电子设备及存储介质 | |
| Meidow et al. | Algebraic reasoning for the enhancement of data-driven building reconstructions | |
| Janeja et al. | FS/sup 3: a random walk based free-form spatial scan statistic for anomalous window detection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210204 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20210204 |
|
| A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20210312 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210330 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20210524 |
|
| 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: 20210803 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210804 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 6926352 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |