《Journal of Structural Geology》：Experimental fluid intrusion in visco-elasto-plastic media with implications for structural interpretations of weak geological systems
编辑推荐：
断层阵列、储层盖层、岩浆房等地质系统的变形同时包含塑性分量与粘性分量，所形成的地质结构在形态上可能高度相似，例如粘性指进与塑性破裂。由于粘性应力与弹性应力的相互作用，两种机制的转换难以通过力学手段精确量化。实验方法能够探测这一转换过程，从而观测在难以数值求解的
断层阵列、储层盖层、岩浆房等地质系统的变形同时包含塑性分量与粘性分量，所形成的地质结构在形态上可能高度相似，例如粘性指进与塑性破裂。由于粘性应力与弹性应力的相互作用，两种机制的转换难以通过力学手段精确量化。实验方法能够探测这一转换过程，从而观测在难以数值求解的力学机制中渐进变形的特征。研究人员采用Hele-Shaw细胞（一种由两块有机玻璃板构成的准二维装置）与一种聚合物凝胶开展实验。该凝胶具有屈服应力与剪切稀化行为，可通过改变交联密度及pH值调控其流变特性。注入第二种流体后，在粘性指进与塑性破裂之间的过渡区域形成了树枝状结构。实验旨在观察在同时具有粘性与弹性组分的明确体系中，从中心注入点分叉的多条指状体与裂缝的运动学特征。单个指状体或裂缝主要以塑性失稳的形式成核，这种失稳由粘性力与界面张力的共同作用局部化，并与弹性响应相抗衡。基于无量纲化的力平衡关系建立的相图表明，在弱弹性系统中，粘塑性转换发生在狭窄的应力差范围内。该实验结果可辅助解译部分中尺度（露头尺度）地质结构。例如，加拿大不列颠哥伦比亚省维多利亚市一处著名露头的岩浆混合结构中，中尺度岩浆单元界面处同时发育了粘性指进与塑性破裂。在该实验背景下，这两种通常同期发生的变形占据狭窄的应力区间，处于中等刚度与侵入应力的场域中。其他地质现象，如熔岩流、断层带以及大陆架甲烷与二氧化碳逃逸管道，则位于该力学空间的不同区域，各自对应弹性—塑性—粘性转换的不同地质表现形式。
研究背景方面，地质结构通常近似为塑性破坏或粘性流动及其组合的结果。冰川与熔岩可作为粘塑性材料流动，形成裂隙与透镜状包体；大陆架甲烷或二氧化碳沿管状通道逃逸；中尺度（露头尺度）的脉体与矿脉因多种过程作用而呈现不同的扩展速率。然而，许多系统（如多孔沉积物中的注入脉、断层带、岩浆体）的结构虽看似由流体间粘性指进形成，却存在于具整体弹性的变形介质中。这类变形与结构的力学解译困难，原因在于弹性、塑性与粘性应力在主导变形机制的转换点趋于相近的量级，导致形态相似的中尺度结构可能源自完全不同的失稳类型，例如粘性指进与塑性破裂的数值模拟极具挑战性。化学工程领域已广泛开展相关研究，其中Carbopol（卡波姆）凝胶作为可调的“智能”材料，其强度可由交联密度控制，并已应用于地球动力学模拟、盐构造变形及断层慢滑事件研究。本研究采用Carbopol 934，聚焦于指状与裂缝网络从注入系统内部失稳向外围介质演化的运动学过程，利用Hele-Shaw细胞容纳凝胶层，并以甘油水溶液作为注入流体，通过观测结构演化将实验结果与地质现象关联，并探讨弹性、塑性与粘性应力收敛于狭窄力学场的意义。
关键技术方法方面，研究人员使用Hele-Shaw细胞开展流体注入实验，通过调节Carbopol 934凝胶浓度（1%与3%重量比）及pH值（3至6）控制其流变性质，采用TA Discovery HR-2流变仪测定其符合Herschel-Bulkley模型的参数，包括屈服应力τ0、一致性指数k与流动行为指数n，以及剪切模量G。注入流体为95%甘油与水的混合物，粘度恒定，以恒定流量注入凝胶层，实验全程背光成像并记录压力变化，利用粒子图像测速（PIV）追踪凝胶内示踪气泡位移以分析应变场，并通过重复注入测试体系的滞后效应。
研究结果部分，在“Rheology and non-dimensional characterization”小节中，研究人员利用Herschel-Bulkley关系式描述凝胶从弹性/塑性固体到剪切稀化流体的转变，并通过毛细数Ca、韦森伯格数Wi与宾汉数Bn三个无量纲数对粘性、弹性与塑性应力进行表征，引入有效稳定数Ca?衡量系统抵抗变化的能力，并结合Saffman-Taylor不稳定性解释指进的形成机制。“Methods”小节显示，实验在透明Hele-Shaw细胞中完成，凝胶制备通过调节pH控制交联密度，流变测试验证其在低应变下表现为弹性储能模量G′主导，高应变下转为损耗模量G″主导，注入过程中压力随时间平方根增长，符合弹性与粘性力的相互作用规律。“Results”小节表明，软凝胶中形成宽大指状体，部分未贯穿细胞而局限于体积内；随凝胶刚度增加，主侵入通道两侧出现更多近直角分支，且分支在主通道通过后继续生长；注入速率显著影响模式：极低流速下呈“池状”累积—压缩—破裂循环，中等流速形成主通道与正交分支的指进模式，高流速则出现指状体与滑移面混合的复杂网络；指宽测量结果显示多数情况符合屈服应力主导的稳定机制而非牛顿流体表面张力机制；指端形态在软、硬凝胶中均接近抛物线形，但硬凝胶侧分支在弹性回弹作用下会收缩；PIV应变分析揭示主指侵入以推挤为主，侧指生长伴随垂直方向的撕裂；重复注入实验显示初始损伤网络不会完全愈合，二次注入会开辟新路径，体现类似断裂的不可逆性。
讨论部分，“Rheological phase space”小节提出由Ca、Wi、Bn构成的三维相空间，指出显著的侧分支与网状结构集中在Ca与Wi共同主导且Bn>1的区域，即塑性屈服应力小于粘性应力时，随着稳定应力超过粘性应力，破坏局限于弱区，反之则形成广泛的局部破坏模式；弹性引起的滞后效应与损伤不可逆性支持部分侵入具备真实破裂特征。“Application to outcrop magma-mingling structures”小节将实验相空间应用于加拿大不列颠哥伦比亚省维多利亚市Holland Point露头，该露头保存了侏罗纪Bonanza弧的岩浆混合结构，包括岩墙、不规则包体、波动界面与岩脉，尺度与实验匹配，推断其形成于Ca、Bn接近1且Wi中等的应力区间，反映了岩浆系统中粘性、弹性与塑性应力的短暂平衡。“A framework for rheology of weak geological systems”小节进一步构建了以体系刚度（rigidity）为横轴、流体侵入应力为纵轴的框架，将熔岩流、气体渗漏、流体化沉积置于低刚度端，断层、裂缝、盖层破裂置于高刚度端，岩浆混合结构位于中部，为地质结构的应力定量解译提供了力学基础。
结论部分，Hele-Shaw细胞内的Carbopol凝胶注入实验为流体侵入成因的地质结构与系统提供了流变学框架，具体表现为：高粘度/低表面能与有限体强度导致粘性指进；低粘度/高表面能与低体强度形成脉网；成比例的粘度与表面能配合高体强度可能导致流体汇聚；高粘度/低表面能配合高体强度产生复杂破裂结构。该框架允许引入其他稳定或失稳应力，通过已知应力推算未知应力。以Holland Point露头为例，岩浆混合结构兼具裂缝与粘性边界特征，形成于流体、内聚与表面应力相近且时空略有变化的阈值区间。这一对比可推广至更广泛的地质系统，依据其刚度（中尺度的流体或固体行为）与流体相（气体、熔体等）侵入应力的关系进行解译。
