致密砂岩水压瞬态致裂液化格子波尔兹曼及有限元数值研究
摘要
水压瞬态致裂液化在地应力测量、地震破坏评估和机理研究、油气(天然气、页岩气)及地热资源开发等地学领域都具有重要的理论和应用价值。自20世纪60年代至今在理论和实际应用方面取得了一些奠基性和开创性成果,但鉴于问题复杂性,瞬态致裂液化机理至今尚不清楚。本研究应用格子波尔兹曼及有限元多孔介质流固耦合物理模型,对地震波载荷作用下致密砂岩水压瞬态致裂液化过程进行数值模拟研究。首先,以鄂尔多斯盆地某油田延长组致密砂岩为例,利用X射线CT断层成像技术,应用基于量子力学第一性原理格子波尔兹曼方法,建立致密砂岩数字岩芯模型。进而,推导格子波尔兹曼及有限元多孔介质流固耦合数值模型公式,建立致密砂岩水压致裂液化物理模型。最后,应用地震波载荷边界和初始条件,模拟瞬态致裂液化流固耦合过程,讨论了地震波载荷幅值、频率及作用时间对致密砂岩孔隙结构(孔隙度大小及连通性)、致密砂岩破裂最大主应力之间关系,得到了地震波作用下致密砂岩致裂液化准则。
Hydrofracturing-liquefaction mechanism is the basic theory for understanding the in-situ stress measurement,the petroleum-gas(nature gas and shale gas)-geothermic development,and earthquake damage evaluation and mechanism study.Due to the complexity and challenge of this issue,mechanism of dynamic hydrofracturing-liquefaction is still not clear,even if a lot of fundamental and landmark achievements had been obtained in this field since the 1960′s.In this work,transient hydrofracturing-liquefaction of tight sandstone under earthquake wave was studied by using hybrid lattice Boltzmann and finite element fluid-solid coupled modeling on parallel GPU-CPU environment.First,the tight sandstone sample from the Ordos Basin Triassic formation was selected and different tomography resolution data was obtained by X-ray CT digital technology.Then,based on the D3Q27 lattice Boltzmann and finite element method,the fluid-solid coupled porous medium physical modeling was established,and the relatively hydrofracturing-liquefaction of tight sandstone was obtained for the first time.Last,the hydrofracturing-liquefaction process was simulated,and the relationship between the amplitude,frequency and time of the earthquake wave and the tight sandstone pore-network evolution,maximum tight sandstone fracturing-liquefaction stress criteria was explored.
Hydrofracturing-liquefaction mechanism is the basic theory for understanding the in-situ stress measurement,the petroleum-gas(nature gas and shale gas)-geothermic development,and earthquake damage evaluation and mechanism study.Due to the complexity and challenge of this issue,mechanism of dynamic hydrofracturing-liquefaction is still not clear,even if a lot of fundamental and landmark achievements had been obtained in this field since the 1960′s.In this work,transient hydrofracturing-liquefaction of tight sandstone under earthquake wave was studied by using hybrid lattice Boltzmann and finite element fluid-solid coupled modeling on parallel GPU-CPU environment.First,the tight sandstone sample from the Ordos Basin Triassic formation was selected and different tomography resolution data was obtained by X-ray CT digital technology.Then,based on the D3Q27 lattice Boltzmann and finite element method,the fluid-solid coupled porous medium physical modeling was established,and the relatively hydrofracturing-liquefaction of tight sandstone was obtained for the first time.Last,the hydrofracturing-liquefaction process was simulated,and the relationship between the amplitude,frequency and time of the earthquake wave and the tight sandstone pore-network evolution,maximum tight sandstone fracturing-liquefaction stress criteria was explored.
引文
[1]Grgic D,Amitrano D.Creep of a porous rock and associated acoustic emission under different hy-drous conditions[J].Journal of Geophysical Research,2009,114:B10201.
[2]Matsuki K.Theoretical examination of the method for measuring three-dimensional in-situ stresseswith anelastic strain recovery of rock core[J].Journal of Mining and Materials Processing Instituteof Japan,1992,108(1):41-45.
[3]Kilb D,Gomberg J,Bodin P.Aftershock triggering by complete Coulomb stress changes[J].Jour-nal of Geophysical Research,2002,107(B4):ESE 2-1-ESE 2-14.
[4]Saar M O,Manga M.Seismicity induced by seasonal groundwater recharge at Mt.Hood,Oregon[J].Earth and Planetary Science Letters,2003,214(3):605-618.
[5]Aboudi J.Micromechanical analysis of fully coupled electro-magneto-thermo-elastic multiphase com-posites[J].Smart materials and structures,2001,10(5):867-877.
[6]Pérez-Aparicio J L,Sosa H.A continuum three-dimensional,fully coupled,dynamic,non-linear fi-nite element formulation for magnetostrictive materials[J].Smart materials and structures,2004,13(3):493-502.
[7]Zhang T,Shi B,Guo Z,et al.General bounce-back scheme for concentration boundary condition inthe lattice-Boltzmann method[J].Physical Review E,2012,85(1):016701.
[8]Fang D N,Jiang B,Hwang K C.A model for predicting effective properties of piezocomposites withnon-piezoelectric inclusions[J].Journal of elasticity and the physical science of solids,2001,62(2):95-118.
[9]Fang D,Berkovits A.Fatigue design model based on damage mechanisms revealed by acoustic emis-sion measurements[J].Journal of engineering materials and technology,1995,117(2):200-208.
[10]Ma J,Wu K,Jiang Z,et al.SHIFT:An implementation for lattice Boltzmann simulation in low-po-rosity porous media[J].Physical Review E,2010,81(5):056702.
[11]Jiang Z,Wu K,Couples G D,et al.The Impact of Pore Size and Pore Connectivity on Single‐Phase Fluid Flow in Porous Media[J].Advanced Engineering Materials,2011,13(3):208-215.
[12]Wang C Y.Liquefaction beyond the near field[J].Seismological Research Letters,2007,78(5):512-517.
[13]Popescu R.Finite element assessment of the effects of seismic loading rate on soil liquefaction[J].Canadian geotechnical journal,2002,39(2):331-344.
[14]Ghosh B,Madabhushi S P G.A numerical investigation into effects of single and multiple frequencyearthquake motions[J].Soil Dynamics and Earthquake Engineering,2003,23(8):691-704.
[15]Ishihara K.Soil behaviour in earthquake geotechnics[M].Clarendon Press;Oxford UniversityPress,1996.350.
[2]Matsuki K.Theoretical examination of the method for measuring three-dimensional in-situ stresseswith anelastic strain recovery of rock core[J].Journal of Mining and Materials Processing Instituteof Japan,1992,108(1):41-45.
[3]Kilb D,Gomberg J,Bodin P.Aftershock triggering by complete Coulomb stress changes[J].Jour-nal of Geophysical Research,2002,107(B4):ESE 2-1-ESE 2-14.
[4]Saar M O,Manga M.Seismicity induced by seasonal groundwater recharge at Mt.Hood,Oregon[J].Earth and Planetary Science Letters,2003,214(3):605-618.
[5]Aboudi J.Micromechanical analysis of fully coupled electro-magneto-thermo-elastic multiphase com-posites[J].Smart materials and structures,2001,10(5):867-877.
[6]Pérez-Aparicio J L,Sosa H.A continuum three-dimensional,fully coupled,dynamic,non-linear fi-nite element formulation for magnetostrictive materials[J].Smart materials and structures,2004,13(3):493-502.
[7]Zhang T,Shi B,Guo Z,et al.General bounce-back scheme for concentration boundary condition inthe lattice-Boltzmann method[J].Physical Review E,2012,85(1):016701.
[8]Fang D N,Jiang B,Hwang K C.A model for predicting effective properties of piezocomposites withnon-piezoelectric inclusions[J].Journal of elasticity and the physical science of solids,2001,62(2):95-118.
[9]Fang D,Berkovits A.Fatigue design model based on damage mechanisms revealed by acoustic emis-sion measurements[J].Journal of engineering materials and technology,1995,117(2):200-208.
[10]Ma J,Wu K,Jiang Z,et al.SHIFT:An implementation for lattice Boltzmann simulation in low-po-rosity porous media[J].Physical Review E,2010,81(5):056702.
[11]Jiang Z,Wu K,Couples G D,et al.The Impact of Pore Size and Pore Connectivity on Single‐Phase Fluid Flow in Porous Media[J].Advanced Engineering Materials,2011,13(3):208-215.
[12]Wang C Y.Liquefaction beyond the near field[J].Seismological Research Letters,2007,78(5):512-517.
[13]Popescu R.Finite element assessment of the effects of seismic loading rate on soil liquefaction[J].Canadian geotechnical journal,2002,39(2):331-344.
[14]Ghosh B,Madabhushi S P G.A numerical investigation into effects of single and multiple frequencyearthquake motions[J].Soil Dynamics and Earthquake Engineering,2003,23(8):691-704.
[15]Ishihara K.Soil behaviour in earthquake geotechnics[M].Clarendon Press;Oxford UniversityPress,1996.350.
版权所有：© 2023 中国地质图书馆 中国地质调查局地学文献中心