[논문 리뷰] Fluid Flow Complexity in Fracture Networks: Analysis with Graph Theory and LBM

H. O. Ghaffari, M. H. B. Nasseri외 1명|arXiv (Cornell University)|2011. 07. 25.

Lattice Boltzmann Simulation Studies참고 문헌 21인용 수 15

한 줄 요약

이 연구는 라티스 보르츠만 방법(LBM)과 그래프 이론을 사용하여 합성 균열 네트워크 내 유체 흐름의 복잡성을 조사한다. 유동성 전이가 허브 성장률에 의해 지배되며, 변환된 네트워크에서 평균 링크 길이와 유동성 간에 역상관관계가 있음을 밝혀내었고, 모듈성과 모티프 스케일링을 통해 공간적으로 비균일한 균열 시스템에서의 유체 흐름 제도를 식별하였다.

ABSTRACT

Through this research, embedded synthetic fracture networks in rock masses are studied. To analysis the fluid flow complexity in fracture networks with respect to the variation of connectivity patterns, two different approaches are employed, namely, the Lattice Boltzmann method and graph theory. The Lattice Boltzmann method is used to show the sensitivity of the permeability and fluid velocity distribution to synthetic fracture networks' connectivity patterns. Furthermore, the fracture networks are mapped into the graphs, and the characteristics of these graphs are compared to the main spatial fracture networks. Among different characteristics of networks, we distinguish the modularity of networks and sub-graphs distributions. We map the flow regimes into the proper regions of the network's modularity space. Also, for each type of fluid regime, corresponding motifs shapes are scaled. Implemented power law distributions of fracture length in spatial fracture networks yielded the same node's degree distribution in transformed networks. Two general spatial networks are considered: random networks and networks with "hubness" properties mimicking a spatial damage zone (both with power law distribution of fracture length). In the first case, the fractures are embedded in uniformly distributed fracture sets; the second case covers spatial fracture zones. We prove numerically that the abnormal change (transition) in permeability is controlled by the hub growth rate. Also, comparing LBM results with the characteristic mean length of transformed networks' links shows a reverse relationship between the aforementioned parameters. In addition, the abnormalities in advection through nodes are presented.

연구 동기 및 목표

• 합성 균열 네트워크 내 연결 패턴이 유체 흐름 복잡성에 미치는 영향를 이해하기 위해.
• 그래프 이론을 통해 네트워크 구조(예: 유동성 및 속도 분포)와 같은 유체 운반 특성 간의 관계를 정량화하기 위해.
• 특히 모듈성과 허브 형성과 같은 비정상적인 유체 흐름 행동을 제어하는 핵심적인 구조적 특징을 규명하기 위해.
• 예측 모델링을 위해 공간적 균열 네트워크 구조와 그에 상응하는 그래프 이론적 표현 간의 매핑을 수립하기 위해.
• 균열 길이 분포와 공간적 구성(무작위 대비 허브 중심/손상 영역)이 흐름 제도 전이에 미치는 영향을 조사하기 위해.

제안 방법

• 유동성 및 속도 필드를 계산하기 위해 합성 균열 네트워크 내에서 라티스 보르츠만 방법(LBM)을 사용해 유체 흐름을 시뮬레이션한다.
• 균열을 노드로, 교차점을 간선으로 간주하여 공간적 균열 네트워크를 추상적 그래프로 변환한다.
• 노드의 차수 분포, 모듈성, 부분그래프 분포와 같은 그래프 특성을 분석하여 구조적 복잡성을 캡처한다.
• 공간 네트워크의 균열 길이에 힘의 법칙 분포를 적용하고, 이를 변환된 그래프의 등가 차수 분포로 매핑한다.
• 모듈성 공간 내에서 흐름 제도 영역을 식별하고, 해당 모티프 형태를 스케일링하여 분류한다.
• 특히 그래프 표현에서의 평균 링크 길이와 LBM 시뮬레이션에서 유도된 유동성 간의 역관계를 드러내기 위해 LBM 결과와 주요 구조적 지표를 비교한다.

실험 결과

연구 질문

• RQ1 합성 균열 네트워크의 연결 패턴이 유체 유동성과 속도 분포에 어떻게 영향을 미치는가?
• RQ2 그래프 표현에서의 구조적 특징(예: 모듈성 및 허브 형성)이 실제 유체 흐름 행동을 어느 정도 반영하는가?
• RQ3 그래프 표현에서의 평균 링크 길이와 LBM 시뮬레이션에서 유도된 유동성 간의 관계는 무엇인가?
• RQ4 다른 공간적 구성(무작위 네트워크 대비 허브 중심 손상 영역)이 흐름 제도 전이에 어떻게 영향을 미치는가?
• RQ5 공간 네트워크에서의 힘의 법칙 분포를 가진 균열 길이가 해당 그래프 네트워크에서 등가의 차수 분포를 생성할 수 있는가?

주요 결과

• 유동성의 이례적인 변화는 네트워크 내 허브 성장률에 의해 수치적으로 제어되며, 중심 노드의 결정적 역할을 시사한다.
• 그래프 표현에서의 평균 링크 길이와 LBM 시뮬레이션에서 유도된 유동성 간에 역관계가 존재한다.
• 흐름 제도는 네트워크의 모듈성 공간 내에서 구분 가능한 영역으로 매핑될 수 있으며, 이는 흐름 행동의 분류를 가능하게 한다.
• 각 흐름 제도와 관련된 모티프 형태는 스케일 가능하며, 그래프 내 특정 구조적 구성과 대응한다.
• 공간 네트워크에서 힘의 법칙 분포를 가진 균열 길이가 변환된 그래프에서 등가의 노드 차수 분포를 생성하며, 이는 구조적 매핑의 타당성을 검증한다.
• 허브 중심(손상 영역 유사) 네트워크는 무작위 네트워크와 비교해 훨씬 다른 유동 특성을 보이며, 특히 흐름 비균형성과 유동성 이질성 측면에서 뚜렷한 차이를 보인다.