Revisiting the hydromechanical formulation of a micromechanics-based phase-field model for poro-elastoplastic media

지표면 깊은 곳, 석유, 가스, 지열을 머금고 있는 다공성 암석 내부에서는 고체 암석과 그 미세한 구멍 속에 갇힌 유체 사이에서 끊임없는 침묵의 전투가 벌어지고 있습니다. 엔지니어들이 이러한 자원을 추출하려고 할 때, 그들은 종-종 암석을 강제로 벌리기 위해 유체를 주입하는데, 이를 수압 파쇄(hydraulic fracturing)라고 합니다. 수십 년 동안 과학자들은 암석을 늘어났다가 다시 돌아오는 단순한 탄성 스펀지로 취급하여 이 과정을 모델링해 왔습니다. 그러나 실제 암석은 더 복잡합니다. 특히 높은 압력에 의해 눌릴 때, 암석은 영구적으로 휘거나 부서질 수도 있고, 내부 균열을 따라 미끄러질 수도 있습니다. 고체 골격과 내부의 유체가 서로 밀고 당길 때 암석이 정확히 어떻게 부서지는지 이해하는 것은 시추를 안전하게 수행하고, 폐기물 저장소를 관리하며, 심지어 빙하가 갈라지거나 마그마가 지각을 밀고 올라오는 것과 같은 자연적인 지질학적 사건을 예측하는 데 매우 중요합니다. 과제는 불가능한 계산에 빠지지 않으면서도 암석이 부서지는 복잡한 현실을 포착할 수 있는 수학적 기술을 만드는 데 있습니다.

한 연구팀은 이 부서짐 과정을 설명하는 데 사용되는 기본 방정식들을 재검토하여, 기존 모델에서 흔히 사용되는 하나의 지름길이 상당한 오류를 초랩한다는 것을 발견했습니다. 암석 역학의 세계에서 과학자들은 '위상장(phase-field)' 접근 방식을 사용하는데, 이는 균열을 날카롭고 보이지 않는 선이 아니라 재료가 점진적으로 약화되는 흐릿한 영역으로 취급하는 방식입니다. 이를 통해 컴퓨터는 균열이 어떻게 성장하고 형태가 변하는지를 시뮬레이션할 수 있습니다. 이를 정확하게 수행하려면 모델은 암석이 당겨질 때(인장)와 눌릴 때(압축)의 두 가지 뚜렷한 거동을 모두 고려해야 합니다. 인장 시에는 미세한 균열이 열립니다. 압축 시에는 균열이 닫히고 그 내부의 거친 표면들이 서로 미끄러지며 마찰을 일으킵니다. 연구진은 산업계와 학계에서 널리 사용되는 기존 모델들이 시스템의 에너지를 설명하는 변수들을 혼동하고 있다는 사실을 발견했습니다. 구체적으로, 그들은 유체의 압력을 에너지 방정식의 직접적인 입력값으로 취급하는 것이, 구멍 속의 유체 양이 아닌 유체 압력을 사용하는 것이 수학적 '도약(jump)' 또는 불연속성을 만든다는 것을 발견했습니다. 이 도약은 모델이 인장에서 압축으로 전환될 때 암석의 강도가 급격하고 비현실적으로 변한다고 예측하게 만들며, 결과적으로 균열이 어디서 어떻게 형성될지에 대한 부정확한 예측을 초래합니다.

이를 해결하기 위해 연구진은 암석의 변형과 유체 함량 사이의 관계를 세심하게 추적하는 새로운 일관된 정식화를 개발했습니다. 그들은 암석을 유체로 채워진 작은 동전 모양의 균열들의 집합체로 상상하며 암석의 미시적 관점에서 모델을 구축했습니다. 이 미세 균열들이 다양한 응력 하에서 어떻게 열리고 닫히는지 분석함으로써, 그들은 압력이 변하더라도 암석의 강도 표면(강도가 결정되는 수학적 경계)이 매끄럽고 연속적으로 유지되도록 하는 일련의 규칙을 도출했습니다. 그들은 암석의 영구적인 변형을 가해지는 힘과 직접 연결하는 특정 유형의 유동 법칙을 사용하고, 에너지 방정식에 누락된 결합 항을 올바르게 포함함으로써, 모델이 인장에서 압축으로의 전이를 정확하게 설명할 수 있음을 증명했습니다. 이 수정은 단순한 미세 조정이 아닙니다. 이는 모델이 균열을 앞으로 밀어내는 구동력을 계산하는 방식을 근본적으로 바꿉니다.

연구팀은 새로운 모델이 잘 작동하는지 확인하기 위해 매우 다른 두 가지 시나리오를 대상으로 테스트를 진행했습니다. 첫째, 유체를 주입하여 암석에 균열을 만드는 전형적인 수압 파쇄 벤치마크를 시뮬레이션했습니다. 그들은 자신들의 결과를 이상적인 조건에 대해 도출된 정밀한 수학적 해답인 알려진 해석적 해(analytical solution)와 비교했습니다. 새로운 모델은 기존의 혼합된 정식화들보다 이 해에 훨씬 더 가깝게 일치했으며, 균열이 성장함에 따라 필요한 압력과 균열의 폭을 정확하게 예측했습니다. 기존 모델들은 필요한 유체 압력을 과대평가하는 경ks향이 있었는데, 이러한 결함은 비효율적이거나 위험한 시추 작업으로 이어질 수 있습니다. 둘째, 암석이 양쪽에서 눌리는 상황을 시뮬레이션했는데, 이는 암석이 엄청난 압력을 받는 지하 환경을 모사합니다. 여기서 암석은 단순히 툭 부서지는 것이 아니라, 재료가 스스로를 지나쳐 미끄러지는 전단 띠(shearing band)를 형성하며 소성 변형을 일으켰습니다. 새로운 모델은 기계적 압착에 의한 전단 균열과 유체 주입에 의한 인장 균열을 모두 포착하며 이 복잡한 거동을 성공적으로 재현했습니다.

시뮬레이션 결과는 암석의 소성 변형이 균열 성장의 중요한 장벽 역할을 한다는 것을 보여줍니다. 암석이 영구적으로 변형될 수 있게 되면 더 많은 에너지를 흡수하게 되며, 이는 균열을 열린 상태로 유지하고 성장시키기 위해 내부의 유체 압력이 더 높아야 함을 의미합니다. 이 추가적인 저항은 균열의 전체적인 프로파일을 변화시켜, 순수 탄성 모델이 제시하는 것보다 더 좁고 더 많은 힘을 필요로 하는 균열을 만듭니다. 또한, 압축 테스트에서 모델은 균열이 형성됨에 따라 암석 내의 공극 압력이 어떻게 변하는지를 보여주었습니다. 암석이 변형되어 유체가 이동할 공간을 만들면서, 국부적으로 압력이 떨어지고 더 많은 유체가 손상 구역으로 유입됩니다. 이러한 상호작용은 유체 흐름과 암석의 기계적 파괴가 긴밀하게 결합되는 피드백 루프를 생성하며, 새로운 모델은 이러한 역동성을 높은 충실도로 포착합니다.

수학적 기초를 바로잡음으로써, 연구진은 다공성 매체를 다루는 엔지니어와 과학자들에게 더 신뢰할 수 있는 도구를 제공했습니다. 새로운 정식화는 서로 다른 유형의 암석 파괴 사이의 전이를 매끄럽고 물리적으로 현실적이게 만들어, 기존의 시도들을 괴롭혔던 인위적인 도약을 제거합니다. 이러한 정확성은 석유 및 가스 생산 최적화부터 탄소 포집을 위한 안전한 지하 저장소 설계에 이르기까지 다양한 응용 분야에 필수적입니다. 이 연구는 고체와 유체가 상호작용하는 복잡한 시스템에서도, 근저에 깔린 물리학에 대한 세심하고 일관된 접근 방식이 재료가 부서지는 진정한 본질을 밝혀낼 수 있음을 입증합니다. 이제 이 모델은 이론적 역학과 발밑의 복잡하고 소성적인 현실 사이의 간극을 메우며, 더 높은 확신을 가지고 지구 지각의 거동을 예측할 준비가 되었습니다.