致密砂岩水压瞬态致裂液化格子波尔兹曼及有限元数值研究.pdfVIP

第33卷第4期 地 震 V01．33，No．4 2013年10月 EARTHQUAKE ()ct．，2013 致密砂岩水压瞬态致裂液化格子波尔 兹曼及有限元数值研究。 朱伯靖1’2…，刘旭耀1’2，程惠红1‘2，柳畅1’2”，刘善琪1’2， 任天翔1’2，李永兵1’2，石耀霖1’2 (1．中国科学院计算地球动力学重点实验室，北京100049；2．中国国科学院大学地球科学 100049；3．RockMechanics rtmentofEarthSciences， 学院，北京 I。aboratory，Depa UniVeofDurham，DurhamDHl3I。E，UK；4．I。aboratoireDe rsity Geologie， Ecole NormaleSup∈rieure，Paris，France) 摘要：水压瞬态致裂液化在地应力测量、地震破坏评估和机理研究、油气(天然气、页岩气)及 地热资源开发等地学领域都具有重要的理论和应用价值。自20世纪60年代至今在理论和实 际应用方面取得了一些奠基性和开创性成果，但鉴于问题复杂性，瞬态致裂液化机理至今尚不 清楚。本研究应用格子波尔兹曼及有限元多孔介质流固耦合物理模型，对地震波载荷作用下致 密砂岩水压瞬态致裂液化过程进行数值模拟研究。首先，以鄂尔多斯盆地某油田延长组致密砂 岩为例，利用X射线CT断层成像技术，应用基于量子力学第一性原理格子波尔兹曼方法，建 立致密砂岩数字岩芯模型。进而，推导格子波尔兹曼及有限元多孑L介质流固耦合数值模型公 式，建立致密砂岩水压致裂液化物理模型。最后，应用地震波载荷边界和初始条件，模拟瞬态 致裂液化流固耦合过程，讨论了地震波载荷幅值、频率及作用时间对致密砂岩孔隙结构(孔隙 度大小及连通性)、致密砂岩破裂最大主应力之间关系，得到了地震波作用下致密砂岩致裂液 化准则。 关键词：致密砂岩水压致裂液化；地震波；X射线CT断层成像分辨率；格子波尔兹曼及有限元 模型；并行CPu—GPU计算 中图分类号：P315 文献标识码：A 文章编号：1000—3274(2013)04一0071—13 引言 水压致裂液化作为地应力测量，石油、天然气、页岩气、地热等资源开发的基本技术 手段，具有重要的应用价值；同时作为地震破坏评估及地震机理研究重要观测途径之一， 具有重要的理论意义。自20世纪60年代至今，在水压致裂液化理论和实际应用领域都取 *收稿日期：2013一01—11；修改回日期：2013一0530 基金项目： 院长基金A(Y1510lKYOO) 作者简介： 朱伯靖(1974)，男，河北永清人，副教授．博士，主要从事断裂力学及固体地球物理学等研究。 万方数据 地 震 33卷 得了一些奠基性和开创性的成果[1~4]，但鉴于水压致裂问题本身复杂性(破坏为动态瞬态 过程，且伴随岩石本身物性变化)、相关学科发展局限性(塑性力学体系虽已建立，但仍处 于发展完善阶段；塑性断裂力学，特别是动态断裂力学理论体系尚处于建立阶段；水压致 裂液化破坏过程力学模型应变、应力、时问、材料属性问定量关系处于经验公式水平)，客 观上大大限制了对该问题进一步的研究，致使研究水平停留在静态／准静态假设基础上的 岩石弹性破坏最大主应力大小与方向确定上；对于随着水压增加，破坏临界点是否存在非 浸润过程、塑性破坏过程是如何发生的、水压增大到临界点以波形式如何影响、破坏过程 伴随岩石物性变化等问题都没有进行研究，随着高性能计算技术，特别是高性能计算硬