KR102177448B1 - X-선 ct 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법 - Google Patents

X-선 ct 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법 Download PDF

Info

Publication number
KR102177448B1
KR102177448B1 KR1020190098876A KR20190098876A KR102177448B1 KR 102177448 B1 KR102177448 B1 KR 102177448B1 KR 1020190098876 A KR1020190098876 A KR 1020190098876A KR 20190098876 A KR20190098876 A KR 20190098876A KR 102177448 B1 KR102177448 B1 KR 102177448B1
Authority
KR
South Korea
Prior art keywords
volume
voxel
gray
ray
pore size
Prior art date
Application number
KR1020190098876A
Other languages
English (en)
Other versions
KR102177448B9 (ko
Inventor
윤태섭
강동훈
Original Assignee
연세대학교 산학협력단
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 연세대학교 산학협력단 filed Critical 연세대학교 산학협력단
Priority to KR1020190098876A priority Critical patent/KR102177448B1/ko
Application granted granted Critical
Publication of KR102177448B1 publication Critical patent/KR102177448B1/ko
Publication of KR102177448B9 publication Critical patent/KR102177448B9/ko

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

본 발명은 컴퓨터를 포함하는 연산처리수단에 의하여 실행되는 프로그램 형태로 이루어지는 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법에 관한 것으로서, X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법이다. 본 발명은 재료의 내부구조 정보인 X-선 CT 이미지와 시편의 공극 크기 정보를 결합하여 정확하고 일관적인 삼상 분할 결과를 제공하는 효과가 있다.

Description

X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법{Evaluation Method of Three-Dimensional Fluid Flow in the Porous or Fractured Media by X-ray CT images and Pore Size Distribution}
본 발명은 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법에 관한 것이다. 구체적으로는 X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법에 관한 것이다.
다공성 재료와 균열 재료는 대부분 연결된 공극 구조로 포함하고 있으며 이로 인해 여러 가지 유동 현상들을 야기하게 된다.
전산 유체 역학 방법과 CT 장비들의 발달로 최근 다공성 재료 및 균열 재료에서의 유동을 공극 규모(pore scale)에서 수치적으로 해석하여 평가하는 방법들이 제안되고 있다(도 1 참조).
일반적으로 다공성 재료 또는 균열 재료를 대상으로 X-선(X-ray) CT 이미징을 사용한다. X-선 CT 이미징으로부터 재료 내부의 이미지를 획득할 수 있으며, 이때 각 픽셀(pixel)이 가지는 밝기 값(또는 CT number)는 해당 지점에서의 밀도와 유효 원자 번호의 함수로 표현된다.
일반적으로 매질이 비슷한 유효 원자 번호들을 가지는 경우 X-선 흡수계수(CT number)와 밀도를 비례관계로 근사(approximate)할 수 있다.
따라서, 많은 경우 이진 분할(Binary segmentation) 방법을 통해 재료 내부의 X-ray CT 이미지를 고체와 공극 이미지로 분할한다. 그리고 해당 공극 이미지로부터 유동 해석에 필요한 도메인 메쉬를 생성하여 전통적인 전산 유체 역학 방법들 (유한 요소 해석, 유한 차분 해석, 유한 체적 해석 등) 또는 격자 볼츠만 방법을 이용하여 편미분 방정식인 나비에 스토크스 방정식 (Navier-Stokes equation) 또는 스토크스 방정식 (Stokes equation)을 풀게 된다. 그 결과로서 재료 내부에서 발생되는 유동을 예측하게 되는 것이다.
참고로, 도 2는 X-선 CT 이미징을 통한 일반적인 유동 해석 과정으로서의 격자 볼츠만 방법을 나타낸다.
한편, 연속방정식, 나비에-스토크스 방정식 및 스토크스 방정식의 예시는 다음과 같다.
[연속방정식]
Figure 112019082999326-pat00001
[나비에-스토크스 방정식]
Figure 112019082999326-pat00002
[스토크스 방정식]
Figure 112019082999326-pat00003
여기서 u는 국부유속(local velocity), ρ는 밀도(density), μ는 절대점성계수(absolute viscosity), 그리고 p는 국부압력(local pressure)이다.
그리고 이때 이진 분할은 X-선 CT 이미지에서 밝기 값이 특정임계값(threshold value)보다 작은 픽셀들을 공극으로 큰 픽셀들을 고체로 보는 전역 문턱치 처리 방법이 대표적이다(도 3 참조).
하지만 X-선 CT의 해상도와 비슷하거나 작은 공극들은 고체와 공극이 하나의 픽셀 위에 혼재하게 되어, X-선 흡수계수(CT number)가 X-선 CT 이미지의 밝기 값 히스토그램에서 명확하게 기대되는 공극과 고체 사이의 애매한 위치에 존재하게 된다.
이진 분할 방법 종류와 히스토그램의 형태에 따라 해당 픽셀이 공극 또는 고체로 평가되어 공극 또는 고체를 과평가 하게 되는 문제점이 발생된다.
이러한 픽셀을 회색 픽셀 또는 회색 복셀이라 칭하며, 종래의 이진 분할과 이를 이용한 유동해석의 한계점은 명확한 것이다.
그리고 CT 해상도가 감소하여 회색 복셀이 증가할수록 계산되는 유동의 결과는 실제 실험 값과 맞지 않게 되는 문제점이 있다(도 4 참조).
따라서, 이러한 한계점을 극복하기 위해 회색 복셀을 고려하는 유동 해석이 필요하게 되었고, 스토크스-브링크만 방정식(Stokes-Brinkman equation)을 이용한 해석 방법이 제안되었다(도 5 참조).
스토크스-브링크만 방정식에서는 해석시 회색 복셀을 지나게 되는 유체에 일정한 저항력을 부여하는 것이 가능하다.
Figure 112019082999326-pat00004
위 식에서, 두 번째 항의 μe는 유효점성계수이고, 세 번째 항이 회색 복셀에서의 점성 저항력을 나타내며 분모의 km는 미시투수계수이며 기존 연구들(Soulaine et al., 2016)에서는 일반적으로 다음과 같은 Kozeny-Carman 식을 사용하였다.
Figure 112019082999326-pat00005
여기서 d는 특성 길이이고
Figure 112019082999326-pat00006
은 공극률이다. 위의 Kozeny-Carman 식에서 특성 길이는 원래 흙 입자들의 입경이나 회색 복셀을 구성하는 입자들의 입경은 구할 수 없으므로, 기존의 연구들에서는 CT 이미지의 해상도 즉 복셀의 한 변의 길이인 1dx 값을 사용해 왔다.
하지만 이 방법은 점성 저항력 항의 계수인 미시투수계수가 dx값에 의존하게 되어 일관성과 정확성에 큰 문제점을 야기하였다.
예를 들어, 미시투수계수는 해당 지점을 구성하는 매질의 특성 계수로 동일한 구조를 가진다면 복셀의 크기와 관계없이 일정해야 한다. 하지만 CT 해상도에 따라 일정하지 않은 문제점이 있는 것이다.
따라서 스토크스-브링크만 방정식을 사용하더라도 여전히 (1)공극 복셀, 회색 복셀, 고체 복셀을 분할하는 방법과 (2)회색 복셀에서의 저항력 계산에 필요한 미시투수계수(micro permeability)를 추정하는 방법이 부정확하고 일관적이지 못한 문제점이 제기되었다.
(문헌 1) 대한민국 공개특허공보 제10-2016-0069379호(2016.06.16)
본 발명에 따른 X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법은 다음과 같은 해결과제를 가진다.
첫째, 본 발명은 X-선 CT 이미지와 공극 크기 분포 정보를 결합하여 명확한 공극 복셀, 회색 복셀, 명확한 고체 복셀의 3상 분리 방법(ternary segmentation)을 제안하고자 한다.
둘째, 회색 복셀에서의 투수계수(voxel permeability)의 계산에서 필요한 특성길이인 공극 크기를 시편의 공극 크기 분포로부터 계산하는 방법을 제안하고자 한다.
본 발명의 해결과제는 이상에서 언급한 것들에 한정되지 않으며, 언급되지 아니한 다른 해결과제들은 아래의 기재로부터 당업자에게 명확하게 이해되어질 수 있을 것이다.
본 발명은 컴퓨터를 포함하는 연산처리수단에 의하여 실행되는 프로그램 형태로 이루어지는 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법에 관한 것으로서, X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법이다.
본 발명은 데이터 준비부(10)에서, 대상 재료의 X-선 CT 이미지를 확보하고, 관심체적(VOI)을 추출하고, 관심체적(VOI)의 부피(V)를 계산하고, 시편 기공 측정의 시험 데이터를 준비하는 S1 단계; 제1 부피 계산부(20)에서, 관심체적 내의 공극의 부피(Vp)와, 고체의 부피(Vs)를 계산하는 S2 단계; 분리직경값 설정부(30)에서, CT 이미지 해상도를 고려한 분리 직경(ds) 값을 설정하는 S3 단계; 제2 부피 계산부(40)에서, 상기 분리 직경(ds) 값을 기준으로 관심체적의 명확한 공극의 부피(Vap)와 회색 공극의 부피(Vgp)를 산출하는 S4 단계; 제1 X-선 흡수계수 산출부(50)에서, X-선 흡수계수 히스토그램에서 낮은 X-선 흡수계수를 갖는 복셀들의 부피들을 더해가며, 명확한 공극의 부피(Vap)와 가장 근접한 X-선 흡수계수(Ip)를 구하는 S5 단계; 제2 X-선 흡수계수 산출부(60)에서, X-선 흡수계수 히스토그램에서 명확한 고체 복셀의 가장 작은 X-선 흡수계수(Is)를 산출하는 S6 단계; 맵핑부(70)에서, 관심체적(VOI)에서 X-선 흡수계수에 대한 공극복셀, 회색복셀 및 명확한 고체복셀의 삼상 분할을 완료하고, 복셀 공극률을 맵핑하는 S7 단계; 공극크기 계산부(80)에서, 복셀 공극률에 대한 회색 공극 크기를 계산하는 S8 단계; 공극면적 계산부(90)에서, 복셀 공극률에 대한 복셀 투수계수를 계산하고, 복셀 공극률에 대한 회색 복셀의 만곡성을 계산하고, 복셀 공극률에 대한 회색 복셀의 공극 면적을 계산하는 S9 단계; 및 수치해석부(100)에서, 스토크스 브링크만 식에 대한 수치유동해석을 수행하는 S10 단계를 포함할 수 있다.
본 발명에 있어서, S1 단계에서는 대상 재료의 X-선 CT 이미지들을 3차원으로 중첩하고 시편 바깥과 시편의 일부 외각을 제외한 나머지 내부를 관심 체적(VOI)으로 선택하며, 관심체적의 부피(V)는 다음 수학식으로 산출될 수 있다.
[수학식 1]
Figure 112019082999326-pat00007
(여기서, dx 는 CT 이미지의 해상도 (즉 1 pixel의 길이)이다.)
본 발명에 있어서, S2 단계에서는 기공측정법으로 측정한 시편의 공극률(φ)을 곱하는 다음 수학식으로, 관심체적 내의 공극의 부피(Vp)와, 고체의 부피(Vs)가 산출될 수 있다.
[수학식 2] Vp = φV and Vs = (1-φ)V
본 발명에 있어서, S3 단계에서는 분리 직경(ds) 값은 1dx ~ 8dx 인 것이 바람직하다.
본 발명에 있어서, S4 단계에서는 누적 공극 크기 분포 곡선에서 상기 분리 직경(ds) 미만의 공극은 회색 공극이며, 상기 분리 직경(ds)을 초과하는 공극은 명확한 공극인 것이 바람직하다.
본 발명에 있어서, S5 단계에서는 X-선 흡수계수(Ip)는 다음 수학식으로 산출될 수 있다.
[수학식 3]
Figure 112019082999326-pat00008
(여기서, NI 는 관심체적에서 X-선 흡수계수가 I 인 복셀들의 개수임)
본 발명에 있어서, S6 단계에서는 회색 복셀들의 각각의 공극률인 복셀 공극률(φv)는 X-선 흡수계수와 비례한다고 가정되며, 관심체적 내의 위치(x)에서의 복셀 공극률(φv)은 다음 수학식으로 산출될 수 있다.
[수학식 4]
Figure 112019082999326-pat00009
(여기서, I(x)는 x지점에서의 X-선 흡수계수이고, Ip와 Is는 각각 명확한 공극 복셀의 가장 큰 X-선 흡수계수와 명확한 고체 복셀의 가장 작은 X-선 흡수계수임)
본 발명에 있어서, 회색 공극들의 공극 부피(Vgp)는 다음 수학식으로 산출될 수 있다.
[수학식 5]
Figure 112019082999326-pat00010
(여기서, NI 는 I의 X-선 감쇠계수를 가지는 복셀들의 수,
Figure 112019082999326-pat00011
v(I)는 I의 X-선 감쇠계수를 가지는 회색복셀의 복셀 공극률임)
본 발명에 있어서, Is 값을 (Ip+1)에서부터 1씩 더하면서, 회색 공극의 부피(Vgp)와, 관심체적 내의 공극의 부피(Vp)를 계산하여, 잠정적인 공극률(tentative porosity)과 잠정적인 회색 공극률 (tentative gray porosity)를 계산하며, 잠정적인 공극률이 실험으로 측정한 시편의 공극률과 일치하게 되는 X-선 흡수계수(Is) 값을 구할 수 있다.
본 발명에 있어서, S7 단계에서는 X-선 흡수계수가 Ip 값보다 작거나 같은 복셀을 명확한 공극 복셀로 정의하고, Ip 값을 초과하고, Is 값 미만인 복셀을 회색 복셀로 정의하고, Is 값보다 크거나 같은 복셀을 명확한 고체 복셀로 정의하여, 삼상 분할이 될 수 있다.
본 발명에 있어서, S8 단계에서는 누적 공극 크기 분포 곡선에서 분리 직경(ds) 미만의 공극은 회색 공극이며, 회색 공극의 공극 크기 분포를 작은 공극에서부터 누적시키고, 관심체적 내의 회색 공극 부피에 대해 다음 수학식으로 나타낼 수 있다.
[수학식 6]
Figure 112019082999326-pat00012
(여기서, Cgp(dgp)는 Cgp 이하의 크기를 가지는 회색 공극들의 누적 부피이며, CF(dgp) 는 공극 크기 분포 곡선에서 전체 공극 부피에 대한 dgp보다 큰 공극들의 누적 부피 비율이고, CF(ds)는 공극 크기 분포 곡선에서 전체 공극 부피에 대한 ds보다 큰 공극들의 누적 부피 비율이다. Vgp는 회색 공극의 전체 부피임)
본 발명에 있어서, 회색 복셀의 공극 크기(dgp)를 추정하기 위해서, 낮은 복셀 공극률을 가지는 회색 복셀은 작은 공극 크기를 가지고, 높은 복셀 공극을 가지는 회색 복셀은 큰 공극을 가진다고 가정될 수 있다.
본 발명에 있어서, S9 단계에서는 상기 복셀 투수계수(kv)는 다음 수학식 7로 산출되며, 상기 회색 복셀의 만곡성(τv)은 다음 수학식 8로 산출되며, 상기 회색 복셀의 공극 면적(Agp)은 다음 수학식 9로 산출될 수 있다.
[수학식 7]
Figure 112019082999326-pat00013
[수학식 8]
Figure 112019082999326-pat00014
[수학식 9]
Figure 112019082999326-pat00015
본 발명은 하드웨어와 결합되어, 본 발명에 따른 X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법을 실행시키기 위하여, 컴퓨터가 판독 가능한 기록매체에 저장된 컴퓨터 프로그램인 것이 가능하다.
본 발명에 따른 X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법은 다음과 같은 효과를 가진다.
첫째, 종래 기술은 일반적으로 이미지 내의 정보만을 사용하여 이상 (공극, 고체) 또는 삼상 (공극, 회색, 고체)을 분할하여 사용하였기 때문에, 이미지의 해상도 및 상태에 따라 그 결과의 일관성이 보장되지 못하였다. 하지만, 본 발명은 재료의 내부구조 정보인 X-선 CT 이미지와 시편의 공극 크기 정보를 결합하여 정확하고 일관적인 삼상 분할 결과를 제공하는 효과가 있다.
둘째, 삼상 분할로부터 얻어진 명확한 공극과 회색 복셀들을 이용하여 스토크스-브링크만 방정식 기반의 유동 해석을 수행하기 위해서는 회색 복셀들의 미시투수계수를 추정해야 한다. 이때 특성 길이가 필요한데, 종래 기술에서는 이 특성길이를 타당성 있게 추정하는 방법이 없었다. 하지만, 본 발명은 전체 시편의 공극 크기 분포로부터 회색 복셀들의 공극 크기 (회색 공극 크기)를 추정하는 방법을 제공하는 효과가 있다.
셋째, 추정된 회색 공극 크기를 이용하여 각 회색 복셀에서의 비 표면적(specirfic surface area)을 추정할 수 있는 효과가 있다. 각 회색 복셀에서의 비 표면적은 반응 유동 (reactive flow) 해석 시 반응 속도와 관계된 중요한 물성이다. 본 발명에서 계산된 전체 비 표면적은 실험으로 측정되는 비 표면적에 가까운 결과를 제공하여 높은 정확도를 제공하는 효과가 있다.
넷째, 도메인 생성에서 명확한 공극과 회색 공극의 비율이 유일하게 결정되는 종래 기술들과는 다르게, 본 발명에서는 분리 직경(separating diameter)이 명확한 공극과 회색 공극의 비율을 결정한다. 이로써, 동일한 CT 이미지에 대해 분리 직경을 달리하여 도메인을 구성할 수 있으며, 이에 따른 회색 공극의 밀도 및 투수 계수, 비 표면적 등이 연동되어 계산되어 정확하고 일관된 결과를 제공하는 효과가 있다.
본 발명의 효과는 이상에서 언급된 것들에 한정되지 않으며, 언급되지 아니한 다른 효과들은 아래의 기재로부터 당업자에게 명확하게 이해되어 질 수 있을 것이다.
도 1은 다공성 재료와 균열 재료들에서의 유동 해석 결과를 나타낸다.
도 2는 X-선 CT 이미징을 통한 일반적인 유동 해석 과정으로서의 격자 볼츠만 방법을 나타낸다.
도 3은 전역 문턱치 처리를 통한 이진 분할을 나타낸다.
도 4는 Berea 사암을 대상으로 한 투수 계수 실험 및 해석들의 결과를 나타낸다.
도 5는 전역적인 삼상 분할을 통한 명확한 공극, 회색, 명확한 고체 복셀들의 분할을 나타낸다.
도 6은 Berea 사암의 누적 공극 크기 분포 곡선을 나타낸다.
도 7은 분리 직경에 따른 명확한 공극의 부피 계산을 나타낸다.
도 8은 관심체적(VOI)의 CT number 히스토그램을 나타낸다.
도 9는 임의로 선택한 Is 값에 따른 공극률 및 정확한 Is 값을 추정하는 것을 나타낸다.
도 10은 삼상 분할된 관심체적(VOI) 내 복셀들의 CT number 히스토그램을 나타낸다.
도 11은 시편의 공극 크기 분포 곡선에서 추출된 회색 공극의 누적 부피 및 해당 공극 크기를 나타낸다.
도 12는 낮은 복셀 공극률을 가지는 복셀들에서부터 누적된 회색 공극 부피를 나타낸다.
도 13은 복셀 공극률에 대한 회색 공극 크기 곡선을 나타낸다.
도 14는 복셀 공극률에 대한 복셀 투수계수 곡선을 나타낸다.
도 15는 격자 볼츠만 방법에 대한 전체 해석 과정을 나타낸다.
이하, 첨부한 도면을 참조하여, 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 용이하게 실시할 수 있도록 본 발명의 실시예를 설명한다. 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 용이하게 이해할 수 있는 바와 같이, 후술하는 실시예는 본 발명의 개념과 범위를 벗어나지 않는 한도 내에서 다양한 형태로 변형될 수 있다. 가능한 한 동일하거나 유사한 부분은 도면에서 동일한 도면부호를 사용하여 나타낸다.
본 명세서에서 사용되는 전문용어는 단지 특정 실시예를 언급하기 위한 것이며, 본 발명을 한정하는 것을 의도하지는 않는다. 여기서 사용되는 단수 형태들은 문구들이 이와 명백히 반대의 의미를 나타내지 않는 한 복수 형태들도 포함한다.
본 명세서에서 사용되는 "포함하는"의 의미는 특정 특성, 영역, 정수, 단계, 동작, 요소 및/또는 성분을 구체화하며, 다른 특정 특성, 영역, 정수, 단계, 동작, 요소, 성분 및/또는 군의 존재나 부가를 제외시키는 것은 아니다.
본 명세서에서 사용되는 기술용어 및 과학용어를 포함하는 모든 용어들은 본 발명이 속하는 기술분야에서 통상의 지식을 가진 자가 일반적으로 이해하는 의미와 동일한 의미를 가진다. 사전에 정의된 용어들은 관련기술문헌과 현재 개시된 내용에 부합하는 의미를 가지는 것으로 추가 해석되고, 정의되지 않는 한 이상적이거나 매우 공식적인 의미로 해석되지 않는다.
본 발명은 컴퓨터를 포함하는 연산처리수단에 의하여 실행되는 프로그램 형태로 이루어진다. 컴퓨터는 적어도 하나의 프로세서, 송수신 장치, 및 메모리 등을 포함할 수 있다. 또한, 컴퓨터는 저장 장치, 입력 인터페이스 장치, 출력 인터페이스 장치 등을 더 포함할 수 있다. 컴퓨터에 포함된 구성 요소들은 버스(bus)에 의해 상호 연결되어 서로 통신을 수행할 수 있다.
본 발명은 컴퓨터를 포함하는 연산처리수단에 의하여 실행되는 프로그램 형태로 이루어지는 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법에 관한 것이다.
본 발명은 X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법로서, 데이터 준비부(10)에서, 대상 재료의 X-선 CT 이미지를 확보하고, 관심체적(VOI)을 추출하고, 관심체적(VOI)의 부피(V)를 계산하고, 시편 기공 측정의 시험 데이터를 준비하는 S1 단계;
제1 부피 계산부(20)에서, 관심체적 내의 공극의 부피(Vp)와, 고체의 부피(Vs)를 계산하는 S2 단계; 분리직경값 설정부(30)에서, CT 이미지 해상도를 고려한 분리 직경(ds) 값을 설정하는 S3 단계;
제2 부피 계산부(40)에서, 상기 분리 직경(ds) 값을 기준으로 관심체적의 명확한 공극의 부피(Vap)와 회색 공극의 부피(Vgp)를 산출하는 S4 단계;
제1 X-선 흡수계수 산출부(50)에서, X-선 흡수계수 히스토그램에서 낮은 X-선 흡수계수를 갖는 복셀들의 부피들을 더해가며, 명확한 공극의 부피(Vap)와 가장 근접한 X-선 흡수계수(Ip)를 구하는 S5 단계;
제2 X-선 흡수계수 산출부(60)에서, X-선 흡수계수 히스토그램에서 명확한 고체 복셀의 가장 작은 X-선 흡수계수(Is)를 산출하는 S6 단계;
맵핑부(70)에서, 관심체적(VOI)에서 X-선 흡수계수에 대한 공극복셀, 회색복셀 및 명확한 고체복셀의 삼상 분할을 완료하고, 복셀 공극률을 맵핑하는 S7 단계;
공극크기 계산부(80)에서, 복셀 공극률에 대한 회색 공극 크기를 계산하는 S8 단계;
공극면적 계산부(90)에서, 복셀 공극률에 대한 복셀 투수계수를 계산하고, 복셀 공극률에 대한 회색 복셀의 만곡성을 계산하고, 복셀 공극률에 대한 회색 복셀의 공극 면적을 계산하는 S9 단계; 및
수치해석부(100)에서, 스토크스 브링크만 식에 대한 수치유동해석을 수행하는 S10 단계를 포함하는 것을 특징으로 한다.
본 발명에 필요한 기본 정보는 다음을 포함한다.
1) 다공성 재료 및 균열 재료의 X-선 CT 이미지 및 기본정보(예, 해상도)
2) 기공 측정법(Porosimetry)으로 측정한 시편의 공극 크기 분포 곡선
3) 시편의 공극률
Figure 112019082999326-pat00016
본 발명에 따른 S1 단계는 데이터 준비부(10)에서, 대상 재료의 X-선 CT 이미지를 확보하고, 관심체적(VOI)을 추출하고, 관심체적(VOI; Volume Of Interest)의 부피(V)를 계산하고, 시편 기공 측정의 시험 데이터를 준비할 수 있다.
본 발명에 따른 S1 단계에서는 대상 재료의 X-선 CT 이미지들을 3차원으로 중첩하고 시편 바깥과 시편의 일부 외각을 제외한 나머지 내부를 관심 체적(VOI)으로 선택하며, 관심체적의 부피(V)는 다음 수학식으로 산출될 수 있다.
Figure 112019082999326-pat00017
여기서, dx 는 CT 이미지의 해상도 (즉 1 pixel의 길이)이다.
본 발명에 따른 S2 단계는 제1 부피 계산부(20)에서, 관심체적 내의 공극의 부피(Vp)와, 고체의 부피(Vs)를 계산할 수 있다.
본 발명에 따른 S2 단계에서는 기공측정법으로 측정한 시편의 공극률(φ)을 곱하는 다음 수학식으로, 관심체적 내의 공극의 부피(Vp)와, 고체의 부피(Vs)가 산출될 수 있다.
Figure 112019082999326-pat00018
공극 크기 분포 곡선을 큰 공극에서부터 누적하여 전체 공극 부피로 정규화한 누적 공극 크기 분포 곡선을 준비할 수 있다(도 6 참조). 도 6은 Berea 사암에 대한 누적 공극 크기 분포 곡선을 나타낸다.
본 발명에 따른 S3 단계는 분리직경값 설정부(30)에서, CT 이미지 해상도를 고려한 분리 직경(ds) 값을 설정할 수 있다.
본 발명에 따른 S3 단계에서는 분리 직경(ds) 값은 1dx ~ 8dx 인이 바람직하다. 여기서, 1dx는 CT 이미지의 해상도 즉 복셀의 한 변의 길이를 의미한다.
본 발명에 따른 S4 단계는 제2 부피 계산부(40)에서, 상기 분리 직경(ds) 값을 기준으로 관심체적의 명확한 공극의 부피(Vap)와 회색 공극의 부피(Vgp)를 산출할 수 있다.
본 발명에 따른 S4 단계의 경우, 누적 공극 크기 분포 곡선에서 상기 분리 직경(ds) 미만의 공극은 회색 공극이며, 상기 분리 직경(ds)을 초과하는 공극은 명확한 공극인 것이 바람직하다.
만약 Berea 사암을 X-선 CT 촬영시 이미지의 해상도가 dx = 2.43μm 이고 분리 직경 ds를 4dx(=9.72μm)로 사용한다면, 명확한 공극의 부피 Vap는 전체 공극 부피의 66.24%가 되고, 회색 공극의 부피 Vgp는 전체 공극 부피의 33.76%가 된다. 즉 Vap = 0.6624Vp 이고 Vgp = 0.3376Vp 이다(도 7 참조).
또한, 관심체적(VOI)의 X-선 흡수계수(CT number)의 히스토그램은 도 8과 같이 나타낼 수 있다.
본 발명에 따른 S5 단계는 제1 X-선 흡수계수 산출부(50)에서, X-선 흡수계수 히스토그램에서 낮은 X-선 흡수계수를 갖는 복셀들의 부피들을 더해가며, 명확한 공극의 부피(Vap)와 가장 근접한 X-선 흡수계수(Ip)를 구할 수 있다.
본 발명에 따른 S5 단계에서 X-선 흡수계수(Ip)는 다음 수학식으로 산출될 수 있다.
Figure 112019082999326-pat00019
(여기서, NI 는 관심체적에서 X-선 흡수계수가 I 인 복셀들의 개수임)
본 발명에 따른 S6 단계는 제2 X-선 흡수계수 산출부(60)에서, X-선 흡수계수 히스토그램에서 명확한 고체 복셀의 가장 작은 X-선 흡수계수(Is)를 산출할 수 있다.
본 발명에 따른 S6 단계에서는 회색 복셀들의 각각의 공극률인 복셀 공극률(φv)는 X-선 흡수계수와 비례한다고 가정되며, 관심체적 내의 위치(x)에서의 복셀 공극률(φv)은 다음 수학식으로 산출될 수 있다.
Figure 112019082999326-pat00020
(여기서, I(x)는 x지점에서의 X-선 흡수계수이고, Ip와 Is는 각각 명확한 공극 복셀의 가장 큰 X-선 흡수계수와 명확한 고체 복셀의 가장 작은 X-선 흡수계수임)
(여기서, Is는 명확한 고체 복셀의 가장 작은 X-선 흡수계수임)
또한, 회색 공극들의 공극 부피(Vgp)는 다음 수학식으로 산출될 수 있다.
Figure 112019082999326-pat00021
여기서, NI 는 I의 X-선 감쇠계수를 가지는 복셀들의 수,
Figure 112019082999326-pat00022
v(I)는 I의 X-선 감쇠계수를 가지는 회색복셀의 복셀 공극률이다.
본 발명에 따른 S6 단계에서는 Is 값을 (Ip+1)에서부터 1씩 더하면서, 회색 공극의 부피(Vgp)와, 관심체적 내의 공극의 부피(Vp)를 계산하여, 잠정적인 공극률(tentative porosity)과 잠정적인 회색 공극률 (tentative gray porosity)를 계산하며, 잠정적인 공극률이 실험으로 측정한 시편의 공극률과 일치하게 되는 X-선 흡수계수(Is) 값을 구할 수 있다(도 9 참조).
본 발명에 따른 S7 단계는 맵핑부(70)에서, 관심체적(VOI)에서 X-선 흡수계수에 대한 공극복셀, 회색복셀 및 명확한 고체복셀의 삼상 분할을 완료하고, 복셀 공극률을 맵핑할 수 있다.
본 발명에 따른 S7 단계에서는 X-선 흡수계수가 Ip 값보다 작거나 같은 복셀을 명확한 공극 복셀로 정의하고, Ip 값을 초과하고, Is 값 미만인 복셀을 회색 복셀로 정의하고, Is 값보다 크거나 같은 복셀을 명확한 고체 복셀로 정의하여, 삼상(3상) 분할이 되는 것이 가능하다.
본 발명에 따른 S8 단계는 공극크기 계산부(80)에서, 복셀 공극률에 대한 회색 공극 크기를 계산할 수 있다.
본 발명에 따른 S8 단계에서는 누적 공극 크기 분포 곡선에서 분리 직경(ds) 미만의 공극은 회색 공극이며, 회색 공극의 공극 크기 분포를 작은 공극에서부터 누적시키고, 관심체적 내의 회색 공극 부피에 대해 다음 수학식으로 나타낼 수 있다.
Figure 112019082999326-pat00023
여기서, Cgp(dgp)는 Cgp 이하의 크기를 가지는 회색 공극들의 누적 부피이며, CF(dgp) 는 공극 크기 분포 곡선에서 전체 공극 부피에 대한 dgp보다 큰 공극들의 누적 부피 비율이고, CF(ds)는 공극 크기 분포 곡선에서 전체 공극 부피에 대한 ds보다 큰 공극들의 누적 부피 비율이다. Vgp는 수학식 5에서 구한 회색 공극의 전체 부피이다.
본 발명에 따른 S8 단계에 있어서, x축을 작은 공극에서부터 누적된 VOI내의 회색 공극의 부피(Cgp) 로 나타내고, y축을 회색 공극의 크기(dgp)로 나타내면, 도 11과 같다.
본 발명에 따른 S8 단계에 있어서, 회색 복셀의 공극 크기(dgp)를 추정하기 위해서, 낮은 복셀 공극률을 가지는 회색 복셀은 작은 공극 크기를 가지고, 높은 복셀 공극을 가지는 회색 복셀은 큰 공극을 가진다고 가정할 수 있다.
따라서.
Figure 112019082999326-pat00024
이면,
Figure 112019082999326-pat00025
가 될 수 있다.
앞의 가정을 통해, 낮은 복셀 공극률을 가지는 복셀들의 공극 부피를 계산하고. 낮은 공극률을 가지는 복셀들에서부터 회색 공극 부피를 누적하면 도 12와 같다.
앞의 두 곡선인 도 11과 도 12는 누적 회색 공극 부피를 각각 x축과 y축으로 공유하므로, 서로 누적 회색 공극 부피를 기준으로 두 곡선을 하나의 곡선으로 만들 수 있다. 도 13은 복셀 공극률에 대한 회색 공극 크기 곡선을 나타낸다.
본 발명에 따른 S8 단계는 공극면적 계산부(90)에서, 복셀 공극률에 대한 복셀 투수계수를 계산하고, 복셀 공극률에 대한 회색 복셀의 만곡성을 계산하고, 복셀 공극률에 대한 회색 복셀의 공극 면적을 계산할 수 있다.
본 발명에 따른 S8 단계에 있어서, 회색 공극의 크기가 알게 되므로, 일반적인 공극률과 공극의 크기를 사용하는 공극률-투수계수 관계에 대한 경험식들을 이용하여 회색 공극들의 미시투수계수 즉 복셀 투수계수 (voxel permeability, kv)를 추정할 수 있다.
예를 들어, 휘어있는 두 면 사이의 기하구조를 가정한 공극률-투수계수 관계식으로 다음 수학식을 사용할 수 있다.
Figure 112019082999326-pat00026
여기서 τv 는 회색 공극의 만곡성 (tortuosity) 이고, 이를 구하기 위하여, 다음과 같은 공극률-만곡성 관계에 대한 경험식을 사용할 수 있다.
Figure 112019082999326-pat00027
최종적으로 도 14와 같이 복셀 공극률에 대한 복셀 투수계수 곡선을 얻을 수 있다.
반응 유동 (reative flow) 해석 시, 고체에서의 반응 속도는 비 표면적의 함수이므로 각 회색 복셀의 비 표면적 계산이 필요하다. 전술한 복셀 투수계수 계산에서 사용한 공극률-투수계수 관계식에서 가정하는 공극의 형상을 사용하거나 단순한 공극 형상들을 사용할 있다.
예를 들어, 휘어있는 두 면 사이의 기하구조를 가정하면, x지점에 있는 회색 복셀 하나의 공극의 면적은 다음 수학식과 같이 계산될 수 있다.
Figure 112019082999326-pat00028
한편, 본 발명은 컴퓨터 프로그램으로 구현될 수도 있다. 구체적으로 본 발명은 하드웨어와 결합되어, 본 발명에 따른 X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법을 실행시키기 위하여, 컴퓨터가 판독 가능한 기록매체에 저장된 컴퓨터 프로그램인 것이 바람직하다.
본 발명에 따른 장치는 컴퓨터로 읽을 수 있는 기록매체에 컴퓨터가 읽을 수 있는 코드로서 구현하는 것이 가능하다. 컴퓨터가 읽을 수 있는 기록매체는 컴퓨터 시스템에 의하여 읽힐 수 있는 데이터가 저장되는 모든 종류의 기록장치를 포함한다. 기록매체의 예로는 ROM, RAM, 광학 디스크, 자기 테이프, 플로피 디스크, 하드 디스크, 비휘발성 메모리 등을 포함한다. 또한 컴퓨터가 읽을 수 있는 기록매체는 네트워크로 연결된 컴퓨터 시스템에 분산되어 분산방식으로 컴퓨터가 읽을 수 있는 코드가 저장되고 실행될 수 있다.
본 명세서에서 설명되는 실시예와 첨부된 도면은 본 발명에 포함되는 기술적 사상의 일부를 예시적으로 설명하는 것에 불과하다. 따라서, 본 명세서에 개시된 실시예들은 본 발명의 기술적 사상을 한정하기 위한 것이 아니라 설명하기 위한 것이므로, 이러한 실시예에 의하여 본 발명의 기술 사상의 범위가 한정되는 것은 아님은 자명하다. 본 발명의 명세서 및 도면에 포함된 기술적 사상의 범위 내에서 당업자가 용이하게 유추할 수 있는 변형 예와 구체적인 실시 예는 모두 본 발명의 권리범위에 포함되는 것으로 해석되어야 할 것이다.
10 : 데이터 준비부
20 : 제1 부피 계산부
30 : 분리직경값 설정부
40 : 제2 부피 계산부
50 : 제1 X-선 흡수계수 산출부
60 : 제2 X-선 흡수계수 산출부
70 : 맵핑부
80 : 공극크기 계산부
90 : 공극면적 계산부
100 : 수치해석부

Claims (14)

  1. 컴퓨터를 포함하는 연산처리수단에 의하여 실행되는 프로그램 형태로 이루어지는 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법에 있어서,
    데이터 준비부(10)에서, 대상 재료의 X-선 CT 이미지를 확보하고, 관심체적(VOI)을 추출하고, 관심체적(VOI)의 부피(V)를 계산하고, 시편 기공 측정의 시험 데이터를 준비하는 S1 단계;
    제1 부피 계산부(20)에서, 관심체적 내의 공극의 부피(Vp)와, 고체의 부피(Vs)를 계산하는 S2 단계;
    분리직경값 설정부(30)에서, CT 이미지 해상도를 고려한 분리 직경(ds) 값을 설정하는 S3 단계;
    제2 부피 계산부(40)에서, 상기 분리 직경(ds) 값을 기준으로 관심체적의 명확한 공극의 부피(Vap)와 회색 공극의 부피(Vgp)를 산출하는 S4 단계;
    제1 X-선 흡수계수 산출부(50)에서, X-선 흡수계수 히스토그램에서 낮은 X-선 흡수계수를 갖는 복셀들의 부피들을 더해가며, 명확한 공극의 부피(Vap)와 가장 근접한 X-선 흡수계수(Ip)를 구하는 S5 단계;
    제2 X-선 흡수계수 산출부(60)에서, X-선 흡수계수 히스토그램에서 명확한 고체 복셀의 가장 작은 X-선 흡수계수(Is)를 산출하는 S6 단계;
    맵핑부(70)에서, 관심체적(VOI)에서 X-선 흡수계수에 대한 공극복셀, 회색복셀 및 명확한 고체복셀의 삼상 분할을 완료하고, 복셀 공극률을 맵핑하는 S7 단계;
    공극크기 계산부(80)에서, 복셀 공극률에 대한 회색 공극 크기를 계산하는 S8 단계;
    공극면적 계산부(90)에서, 복셀 공극률에 대한 복셀 투수계수를 계산하고, 복셀 공극률에 대한 회색 복셀의 만곡성을 계산하고, 복셀 공극률에 대한 회색 복셀의 공극 면적을 계산하는 S9 단계; 및
    수치해석부(100)에서, 스토크스 브링크만 식에 대한 수치유동해석을 수행하는 S10 단계를 포함하는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
  2. 청구항 1에 있어서, S1 단계에서는
    대상 재료의 X-선 CT 이미지들을 3차원으로 중첩하고 시편 바깥과 시편의 일부 외각을 제외한 나머지 내부를 관심 체적(VOI)으로 선택하며,
    관심체적의 부피(V)는 다음 수학식으로 산출되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
    [수학식 1]
    Figure 112019082999326-pat00029

    (여기서, dx 는 CT 이미지의 해상도 (즉 1 pixel의 길이)이다.)
  3. 청구항 1에 있어서, S2 단계에서는
    기공측정법으로 측정한 시편의 공극률(φ)을 곱하는 다음 수학식으로, 관심체적 내의 공극의 부피(Vp)와, 고체의 부피(Vs)가 산출되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
    [수학식 2] Vp = φV and Vs = (1-φ)V
  4. 청구항 1에 있어서, S3 단계에서는
    분리 직경(ds) 값은 1dx ~ 8dx 인 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
  5. 청구항 1에 있어서, S4 단계에서는
    누적 공극 크기 분포 곡선에서 상기 분리 직경(ds) 미만의 공극은 회색 공극이며, 상기 분리 직경(ds)을 초과하는 공극은 명확한 공극인 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
  6. 청구항 1에 있어서, S5 단계에서는
    상기 X-선 흡수계수(Ip)는 다음 수학식으로 산출되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
    [수학식 3]
    Figure 112019082999326-pat00030

    (여기서, NI 는 관심체적에서 X-선 흡수계수가 I 인 복셀들의 개수임)
  7. 청구항 1에 있어서, S6 단계에서는
    회색 복셀들의 각각의 공극률인 복셀 공극률(φv)는 X-선 흡수계수와 비례한다고 가정되며, 관심체적 내의 위치(x)에서의 복셀 공극률(φv)은 다음 수학식으로 산출되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
    [수학식 4]
    Figure 112019082999326-pat00031

    (여기서, I(x)는 x지점에서의 X-선 흡수계수이고, Ip와 Is는 각각 명확한 공극 복셀의 가장 큰 X-선 흡수계수와 명확한 고체 복셀의 가장 작은 X-선 흡수계수임)
  8. 청구항 7에 있어서,
    회색 공극들의 공극 부피(Vgp)는 다음 수학식으로 산출되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
    [수학식 5]
    Figure 112019082999326-pat00032

    (여기서, NI 는 I의 X-선 감쇠계수를 가지는 복셀들의 수,
    Figure 112019082999326-pat00033
    v(I)는 I의 X-선 감쇠계수를 가지는 회색복셀의 복셀 공극률임)
  9. 청구항 8에 있어서,
    Is 값을 (Ip+1)에서부터 1씩 더하면서, 회색 공극의 부피(Vgp)와, 관심체적 내의 공극의 부피(Vp)를 계산하여, 잠정적인 공극률(tentative porosity)과 잠정적인 회색 공극률 (tentative gray porosity)를 계산하며,
    잠정적인 공극률이 실험으로 측정한 시편의 공극률과 일치하게 되는 X-선 흡수계수(Is) 값을 구하는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
  10. 청구항 1에 있어서, S7 단계에서는
    X-선 흡수계수가 Ip 값보다 작거나 같은 복셀을 명확한 공극 복셀로 정의하고,
    Ip 값을 초과하고, Is 값 미만인 복셀을 회색 복셀로 정의하고,
    Is 값보다 크거나 같은 복셀을 명확한 고체 복셀로 정의하여,
    삼상 분할이 되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
  11. 청구항 1에 있어서, S8 단계에서는
    누적 공극 크기 분포 곡선에서 분리 직경(ds) 미만의 공극은 회색 공극이며,
    회색 공극의 공극 크기 분포를 작은 공극에서부터 누적시키고, 관심체적 내의 회색 공극 부피에 대해 다음 수학식으로 나타내는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
    [수학식 6]
    Figure 112019082999326-pat00034

    (여기서, Cgp(dgp)는 Cgp 이하의 크기를 가지는 회색 공극들의 누적 부피이며, CF(dgp) 는 공극 크기 분포 곡선에서 전체 공극 부피에 대한 dgp보다 큰 공극들의 누적 부피 비율이고, CF(ds)는 공극 크기 분포 곡선에서 전체 공극 부피에 대한 ds보다 큰 공극들의 누적 부피 비율이다. Vgp는 회색 공극의 전체 부피임)
  12. 청구항 11에 있어서,
    회색 복셀의 공극 크기(dgp)를 추정하기 위해서,
    낮은 복셀 공극률을 가지는 회색 복셀은 작은 공극 크기를 가지고, 높은 복셀 공극을 가지는 회색 복셀은 큰 공극을 가진다고 가정되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
  13. 청구항 1에 있어서, S9 단계에서는
    복셀 투수계수(kv)는 다음 수학식 7로 산출되며, 상기 회색 복셀의 만곡성(τv)은 다음 수학식 8로 산출되며, 상기 회색 복셀의 공극 면적(Agp)은 다음 수학식 9로 산출되는 것을 특징으로 하는
    X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법.
    [수학식 7]
    Figure 112020502855629-pat00035

    [수학식 8]
    Figure 112020502855629-pat00036

    [수학식 9]
    Figure 112020502855629-pat00037
  14. 하드웨어와 결합되어, 청구항 1에 따른 X-선 CT 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법을 실행시키기 위하여, 컴퓨터가 판독 가능한 기록매체에 저장된 컴퓨터 프로그램.
KR1020190098876A 2019-08-13 2019-08-13 X-선 ct 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법 KR102177448B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020190098876A KR102177448B1 (ko) 2019-08-13 2019-08-13 X-선 ct 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020190098876A KR102177448B1 (ko) 2019-08-13 2019-08-13 X-선 ct 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법

Publications (2)

Publication Number Publication Date
KR102177448B1 true KR102177448B1 (ko) 2020-11-11
KR102177448B9 KR102177448B9 (ko) 2023-05-11

Family

ID=73451747

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020190098876A KR102177448B1 (ko) 2019-08-13 2019-08-13 X-선 ct 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법

Country Status (1)

Country Link
KR (1) KR102177448B1 (ko)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116337708A (zh) * 2022-12-06 2023-06-27 武汉理工大学 提取多孔材料孔隙特征和流体分布的方法与系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101110787B1 (ko) * 2009-11-17 2012-02-16 한국건설기술연구원 엑스레이 씨티촬영을 통한 미세 토사의 간극비 측정방법
KR101370496B1 (ko) * 2013-09-26 2014-03-06 한국건설기술연구원 복합매질로 이루어진 시편에 대한 X-ray CT 영상의 최소 단위에 존재하는 각 순수매질의 부피비 측정방법
KR101525512B1 (ko) * 2014-06-27 2015-06-04 연세대학교 산학협력단 다공성 재료 막힘 현상 시뮬레이션 방법 및 이를 이용한 다공성 재료의 투수 성능 평가 장치 및 방법
KR20160069379A (ko) 2014-12-08 2016-06-16 공주대학교 산학협력단 암석의 박편 영상을 이용한 암석의 투수율 및 공극률 산출 장치 및 그 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101110787B1 (ko) * 2009-11-17 2012-02-16 한국건설기술연구원 엑스레이 씨티촬영을 통한 미세 토사의 간극비 측정방법
KR101370496B1 (ko) * 2013-09-26 2014-03-06 한국건설기술연구원 복합매질로 이루어진 시편에 대한 X-ray CT 영상의 최소 단위에 존재하는 각 순수매질의 부피비 측정방법
KR101525512B1 (ko) * 2014-06-27 2015-06-04 연세대학교 산학협력단 다공성 재료 막힘 현상 시뮬레이션 방법 및 이를 이용한 다공성 재료의 투수 성능 평가 장치 및 방법
KR20160069379A (ko) 2014-12-08 2016-06-16 공주대학교 산학협력단 암석의 박편 영상을 이용한 암석의 투수율 및 공극률 산출 장치 및 그 방법

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"공극 스케일에서 이방성을 가진 다공성 암석의 유체 투과율에 관한 수치적 연구", 2013 한국암반공학회 춘계 총회 및 학술발표회 논문집(pp. 29-35), 2013년 3월 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116337708A (zh) * 2022-12-06 2023-06-27 武汉理工大学 提取多孔材料孔隙特征和流体分布的方法与系统
CN116337708B (zh) * 2022-12-06 2023-11-17 武汉理工大学 提取多孔材料孔隙特征和流体分布的方法与系统

Also Published As

Publication number Publication date
KR102177448B9 (ko) 2023-05-11

Similar Documents

Publication Publication Date Title
Taud et al. Porosity estimation method by X-ray computed tomography
Garing et al. Pore-scale capillary pressure analysis using multi-scale X-ray micromotography
US6516080B1 (en) Numerical method of estimating physical properties of three-dimensional porous media
RU2639727C2 (ru) Системы и способы цифрового анализа горных пород с определением сэо многофазного потока
AU2011258594B2 (en) Method for obtaining consistent and integrated physical properties of porous media
Rabbani et al. Evolution of pore-scale morphology of oil shale during pyrolysis: a quantitative analysis
US20120163688A1 (en) System and method for multi-phase segmentation of density images representing porous media
KR101661336B1 (ko) 암석의 박편 영상을 이용한 암석의 투수율 및 공극률 산출 장치 및 그 방법
US11314909B2 (en) Methods and systems for simulating multiphase flow through porous media
Chukwudozie et al. Prediction of non-Darcy coefficients for inertial flows through the castlegate sandstone using image-based modeling
KR102177448B1 (ko) X-선 ct 이미지 및 공극 크기 분포를 이용한 다공성 재료 또는 균열 재료의 3차원 유동 평가 방법
Hyväluoma et al. Using microtomography, image analysis and flow simulations to characterize soil surface seals
Ozelim et al. Representative elementary volume determination for permeability and porosity using numerical three-dimensional experiments in microtomography data
Meira Cassaro et al. Pore‐Size Distributions of Soils Derived using a Geometrical Approach and Multiple Resolution MicroCT Images
CN108663287B (zh) 一种利用ct图像精确计算煤岩密度的方法
Munawar et al. Upscaling reservoir rock porosity by fractal dimension using three-dimensional micro-computed tomography and two-dimensional scanning electron microscope images
Sun et al. EPCI: A new tool for predicting absolute permeability from computed tomography images
Caubit et al. Evaluation of the reliability of prediction of petrophysical data through imagery and pore network modelling
Wang et al. Multiphase lattice Boltzmann modeling of cyclic water retention behavior in unsaturated sand based on X-ray Computed Tomography
Kato et al. Segmentation of multi-phase X-ray computed tomography images
Kaestner et al. Mapping the 3D water dynamics in heterogeneous sands using thermal neutrons
Scarfone et al. Analytical and numerical modelling of air trapping during wetting of unsaturated soils
CN113167713A (zh) 数字表征岩石渗透率的方法
Yang et al. Ternary segmentation and estimation of permeability for porous rocks based on 3D X-ray computed tomographic images by hidden Markov random field and Brinkman-force lattice Boltzmann model
WO2019151889A1 (en) A method for determining a three-dimensional spatial distribution of porosity in a sample of a heterogeneous porous medium

Legal Events

Date Code Title Description
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
G170 Re-publication after modification of scope of protection [patent]