本文大纲：
一、当褶皱专家邂逅“奥维耶多褶皱分析小组”
二、给褶皱“画像”：从二维素描到三维CT
三、破解折叠“力学”：中性面与两种关键模式
四、量化“挤压力”：被压扁的褶皱记录了总缩短量
五、创新“运动学”：旋转剪切，一个解释相似褶皱的新钥匙
六、实战！解剖大自然的经典褶皱案例
七、微观与伴生构造：小尺度里的大乾坤
八、遗产与未来：当几何学遇见动力学
我们脚下看似亘古不变、坚硬无比的岩石，其实是一部记录了亿万次“温柔”弯曲的史诗。
山脉的脊梁、悬崖的纹理，那些波浪般起伏的岩层，并非静止的雕塑，而是动态变形的化石。
读懂这些“地球的皱纹”，就能破译板块碰撞的密码、追溯矿脉的踪迹，甚至预测能源的所在。
一位英国地质学家，远赴西班牙北部的奥维耶多，与一群数学家和不懈的野外工作者并肩二十载，他们用罗盘、显微镜、数学方程和计算机代码，共同写下了一部关于“岩石为何以及如何弯曲”的传奇。
这不仅是一场科学的远征，更是一段关于好奇心、友谊与智慧的动人故事。
一、当褶皱专家邂逅“奥维耶多褶皱分析小组”
故事始于本世纪初。
在西班牙西北部阿斯图里亚斯地区的历史名城奥维耶多，当地大学的一群地质学家组成了一个专注研究岩石褶皱（Folding）的团队——奥维耶多褶皱分析小组（Oviedo Fold Analysis Group，简称OFAG）。他们试图用更精确的数学工具，来描述和解释自然界中千变万化的岩石弯曲形态。
他们心中有一个理想的合作人选：理查德·J·莱尔（Richard J. Lisle），英国卡迪夫大学的结构地质学教授。
莱尔不仅是经典教材《现代构造地质学技术》的作者之一，更以对褶皱几何学的深刻洞察而闻名。
2004年，在西班牙国家项目的支持下，莱尔教授受邀前往奥维耶多进行了为期六个月的访问。这次访问，点燃了长期合作的火焰。
图1：友谊与合作的野外瞬间。 (a) 2007年，莱尔（中）与OFAG成员在西班牙比利牛斯山是更灵活的工具。
图2 用贝塞尔曲线“画”出褶皱。一条复杂的褶皱翼部（即一个“山坡”），可以用由四个控制点定义的贝塞尔曲线段来完美拟合。这为用数学精确描述自然界中不规则的褶皱形状提供了强大工具。
然而，方法多了反而混乱。
莱尔与OFAG团队洞察到一个关键：所有这些分类参数，本质上都与一个叫归一化面积的量有关。
简单说，就是褶皱曲线与参考线之间包围的面积比例。基于此，他们开发了一款名为“FoldProfiler”的MATLAB程序。地质学家只需在野外拍一张褶皱照片导入，程序就能自动拟合曲线，并将其归类到一个统一的图表中。
图3 褶皱形状的“统一战场”。这张图是褶皱分类的“地图”。横坐标是褶皱的宽高比（胖瘦），纵坐标就是关键的“归一化面积”。不同类型的褶皱（如尖锐的锯齿状、圆滑的正弦状、扁平的抛物线状）会落在图上不同的区域。FoldProfiler让复杂的褶皱分类变得直观、自动化。
但地球是三维的，褶皱像起伏的毯子一样在空间中延伸。
莱尔很早就指出，对于一种常见的等长褶皱（岩层像弯曲一张纸一样，长度不变），其三维曲面必须满足高斯曲率为零，这在几何上叫可展曲面和高斯曲率，为三维褶皱曲面建立了一套清晰的口头分类系统。
图4 三维褶皱的“体检报告”。根据曲面在每个点是像穹窿一样凸起（平均曲率>0，背形）还是像盆地一样凹陷（平均曲率<0，向形），以及像马鞍一样扭曲（高斯曲率<0）还是像球面一样不可展（高斯曲率>0），可以将褶皱细分为多种类型，如“同塑性背形”、“单斜背形”等。这让我们能用精确的语言描述地下结构的每一个局部形态。
为了在复杂的三维构造图上定位最重要的褶皱枢纽线（即“山脊线”或“山谷线”），他们定义了“强枢纽线”和“弱枢纽线”，分别对应曲面最大和最小弯曲程度最高的轨迹，并用FORTRAN程序进行计算和绘制。
图5 在构造等高线图上追踪“山脊线”。这张图展示了地下一个复杂褶皱面的等高线图。图中弯曲的红色线条就是计算出的“强枢纽线”，它们清晰地揭示了褶皱在三维空间中的起伏走向，而这些线并不一定经过等高线最密（最陡）的地方。这项技术对于理解非圆柱状褶皱至关重要。
三、破解折叠“力学”：中性面与两种关键模式
描述几何形态是第一步，理解其背后的力学机制才是核心。
莱尔与OFAG团队聚焦于一个经典模型：切向纵向应变，也叫中性面褶皱。你可以想象弯曲一根塑料尺：外弧被拉伸，内弧被压缩，中间有一个既不拉也不压的“中性面”。
拉姆齐等前辈定义过这个模型，但假设了面积不变。OFAG成员通过严谨的数学推演和“FoldModeler”软件的模拟发现，实际情况更复杂。他们提出了两个更符合现实的端元模式：
1. 平行切向纵向应变：岩层在弯曲时，内部原来垂直于层面的“标杆线”始终保持垂直。这会产生完美的平行褶皱（岩层厚度在弯曲中保持不变），但代价是：外弧面积必须显著增加，内弧面积必须减少。这意味着外弧可能产生张裂隙和矿脉，内弧则可能发生压溶作用形成劈理。
2. 等面积切向纵向应变：即传统的面积守恒模式，但会造成内弧在枢纽区“鼓包”。
图6 平行切向纵向应变（PTLS）的“体检数据”。 (a-c) 展示了岩层从平直到弯曲的过程，中性面呈抛物线形。(d-e) 图表显示，随着岩层倾角增大，应变椭圆轴比（R，代表变形强弱）在外弧急剧增加，面积变化率（J）在外弧大于1（膨胀），在内弧小于1（压缩）。这些定量关系是验证模型的钥匙。
他们指出，自然界中更可能是两种模式的结合或过渡：早期曲率小，可能是等面积模式；随着弯曲加剧，平行模式开始主导以容纳更大的应变。这一认识将几何观察与实际的岩石变形结构（如张裂隙、压力影、缝合线）直接联系起来，甚至能指导矿产勘探和油气储层中裂隙发育带的预测。
四、量化“挤压力”：被压扁的褶皱记录了总缩短量、
在造山带，岩石不仅被弯曲，还被整体“压扁”。
如何从单个压扁的平行褶皱中估算出区域受到的总压缩量，是个难题。传统方法只能估算“压扁”阶段的应变，忽略了之前“弯曲”阶段的缩短。
莱尔与OFAG团队发展了一套创新的两阶段法：首先，用FoldProfiler将褶皱的中线拟合为圆锥曲线，反推出弯曲前的状态和“弯曲缩短量”；然后，用一种新方法（基于Bastida等定义的参数s1, s2）或传统方法估算“压扁量”。两者结合，即得到总体的体应变。
图7 测量“褶皱的记忆”。 (a) 西班牙纳尔塞阿背形中一些小型褶皱的总应变椭圆。(b) 关键发现：这些褶皱的总缩短量与其翼间角（褶皱的张开程度）呈明显的线性负相关。这意味着，在地质图上简单地测量褶皱的紧闭程度，就能对其经历的挤压强度进行初步估算！这项技术在区域构造分析中极具实用价值。
五、创新“运动学”：旋转剪切，一个解释相似褶皱的新钥匙
为了探究更复杂的褶皱，他们必须深入应变分析的领域。除了开发用于三维应变分析的“StrainModeler”软件，团队最大的理论贡献之一是提出了一个全新的均匀变形模式：旋转剪切。
图8 像剪刀一样变形——旋转剪切。 (a) 这种变形就像一个固定剪刀的一边，转动另一边。存在一个不变的方向，和另一个长度不变但旋转的方向。(b) 与经典的简单剪切（如推卡牌）相比，旋转剪切的应变随变形进程的变化规律不同。(c) 旋转剪切可以完美地解释膝折带这种常见构造的形成。
旋转剪切成为理解相似褶皱（褶皱中各岩层形态相似，厚度在枢纽处最大、翼部最小）的关键。
传统认为相似褶皱由跨层的不均匀简单剪切形成，但这难以解释其普遍伴随的垂直缩短。Bobillo-Ares等人提出，层平行缩短加上不均匀的旋转剪切，是更合理的机制。
图9 破解“相似褶皱”之谜。(a) 西班牙一处板岩中发育的完美相似褶皱及小寄生褶皱。(b) 理论模型（层平行缩短+旋转剪切）成功拟合了自然褶皱的形态，图中的线段代表应变椭圆长轴方向。模型显示，变形伴随着显著的体积损失（压溶），这与板岩中发育的板劈理事实相符。
六、实战！解剖大自然的经典褶皱案例
依托强大的理论工具（FoldModeler等）和详实的野外工作，团队对多种特定褶皱类型发起了“总攻”。
1. 锯齿状褶皱：这种尖锐的褶皱需要多种机制顺序或同时作用。他们对英格兰康沃尔和西班牙的经典露头进行研究后，提出了一个可能的演化序列：初期层平行缩短和等面积切向应变启动弯曲；随着曲率增大，平行切向应变和弯流作用（岩层内部剪切）开始主导；后期可能叠加均匀压扁。整个过程中，层间滑动至关重要。
图10 锯齿状褶皱的演化路径。这张图通过FoldModeler模拟，展示了锯齿状褶皱如何从圆滑逐渐变尖锐。横坐标是岩层最大倾角，纵坐标是归一化振幅。曲线揭示了在ETLS-PTLS-FF机制序列下，褶皱形态的定量演化过程。
2. 平卧褶皱：这些规模巨大、近乎躺倒的褶皱是造山带的标志。团队总结认为，其形成通常需要：一个初始的不均匀性（如基底断块）作为“障碍物”；一个以简单剪切为主的变形场；以及后期重力导致的垂向压扁作用。这很好地解释了为什么许多平卧褶皱根部存在一个“根带”。
图11 大尺度平卧褶皱的“实验重现”。(a) 美国阿巴拉契亚山脉真实的平卧褶皱剖面。(b) 用多层蜡和粘土进行的物理模拟实验，因底部存在刚性基底（障碍物）而形成了极其相似的平卧褶皱形态。这证实了“障碍物”触发机制的重要性。
3. 斜层理褶皱分析：当一个先存的、斜交层面的面（如化石沙丘的斜层理）被褶皱时，它的角度变化轨迹包含了丰富的运动学信息。团队系统模拟了不同机制下“倾角-斜交角”关系图，建立了一套“诊断图谱”，成为野外判断褶皱机制的实用技巧。
图12 褶皱机制的“诊断仪”——δ-α图。图中曲线展示了在不同褶皱机制（FF弯流，PTLS平行切向应变）下，一个初始斜交层的斜交角（δ）如何随岩层倾角（α）变化。每一条曲线都有独特的形态，像指纹一样对应着特定的变形历史。
4. 不对称褶皱：通过升级FoldModeler，他们证明了初始的对称弯曲褶皱，在叠加一个非共轴的均匀应变（最大拉伸方向与褶皱轴斜交）后，会自然地演化为不对称形态，并伴随枢纽迁移。
图13 不对称褶皱的“数字拟合”。(a) 西班牙比利牛斯山一个天然不对称向斜。(b) 理论模型完美叠加在照片上。(c-d) 应变椭圆长轴方向（θ）与岩层倾角（α）的关系图，理论（红线）与实测数据（黑点）高度吻合。(e) 拉姆齐分类图上的吻合也验证了模型可靠性。
七、微观与伴生构造：小尺度里的大乾坤
褶皱从不单独出现，总伴随着一系列精彩的伴生构造。
皱纹劈理：这是在显微尺度上的褶皱，通常发育在先存的板劈理之上。团队运用旋转剪切模型，成功解释了其类似相似褶皱的几何形态，以及压溶作用导致的劈理域体积损失和微劈石域体积增加，实现了漂亮的定量拟合。
图14 显微镜下的“小山脉”——皱纹劈理。(a) 板岩中发育的对称皱纹劈理及其等倾线。(b) 通过“旋转剪切+体积补偿”模型在计算机中生成的模拟皱纹劈理，每个单元格中的线段代表应变方向。模型再现了自然界的精细结构。
“锯齿状”构造与弯曲脉体：在西班牙比利牛斯山的 turbidite 地层中，他们发现了一种奇特构造：在褶皱翼部的强硬层底部，发育一系列指向背形枢纽的楔状张性脉体，像锯齿一样。研究认为，这是在弯滑作用主导下，因孔隙流体压力升高，在强硬层底部产生局部张应力并诱发剪切破裂，进而扩展成弯曲脉体的结果。
图15 褶皱翼部的“锯齿”与“弯刀”。(a) 示意图展示强硬层底部的“锯齿状”阶梯及相关的弯曲脉体（底部组和顶部组）。(b) 野外实例，背形枢纽在左侧。这些精美构造是特定应力条件下的忠实记录者。
先存褶皱中的劈理：在西班牙布雷拉的海岸，他们研究了后期劈理如何叠加在已经紧闭的褶皱之上。他们提出了“转换机制”：后期变形中，强硬层像活塞一样相对滑动，挤压夹在其中的软弱层，迫使软弱层发生体积变化（压溶）并形成新的劈理，同时强硬层本身发生张裂。
图16 倒转褶皱中，后期劈理在正常翼和倒转翼与层理的不同交切关系。
图17 弱能干岩石（灰色所示）通过平移机制发生褶皱变形（强能干岩石以点状图案表示）。（a）褶皱初始形态。（b）平移机制作用下形成的变形形态。
八、遗产与未来：当几何学遇见动力学
理查德·莱尔与奥维耶多团队的这段长达二十年的合作，在结构地质学史上留下了浓墨重彩的一笔。他们以严谨的几何学为基石，结合系统的运动学模拟，为纷繁复杂的褶皱现象建立了一套可量化、可预测的解释框架。
他们的工作，关于切向纵向应变两种模式的区分，关于褶皱总缩短量的估算公式，关于旋转剪切的新概念，关于多种褶皱的演化序列，已被广泛写入教科书和综述论文，成为该领域的标准知识。
展望未来，这份遗产是通往“动力学”模拟的桥梁。他们精确描述的位移场和应变场，可以为岩石赋予特定的流变学性质，进而计算应力场，并最终预测裂缝的发育。这将使褶皱研究从“解释过去”迈向“预测属性”，在水文地质、工程地质和油气储层评价中发挥更大的实际价值。
科学的故事，始于对一块弯曲岩石的好奇，成于跨国界、跨学科的智慧碰撞，最终融入我们对这个星球更深刻的理解之中。这，或许就是这段奥维耶多往事，留给我们最珍贵的启示。
声明 | 文中图片及文章版权属于原作者及出版社，本公众号只做基于文章数据和方法的科学解读与二次创作之用，如有侵权，联系即删。
来源 | Bastida F, Aller J, Bobillo-Ares N C, et al. Richard Lisle in Oviedo (Spain): Contributions to the analysis of folding–a review[J]. Journal of Structural Geology, 2025: 105553.
编辑 | 从白垩纪路过
审核 | 微的尘
出品 | 概论地球科学
