花岗岩裂隙储层各向异性高阶有限差分数值模拟  

Anisotropic Numerical Simulation of Fracture Reservoir in Granite with High-order Finite Difference

作　　者：张华[1] 龚育龄[1] 

机构地区：[1]东华理工大学放射性地质与勘探技术国防重点学科实验室,江西抚州344000

出　　处：《东华理工大学学报（自然科学版）》2012年第4期416-421,共6页Journal of East China University of Technology(Natural Science)

基　　金：东华理工大学放射性地质与勘探技术国防重点学科实验室开放基金(2011RGET03);国防科学技术工业局"甘肃北山重点花岗岩体展布及内部结构地球物理探测研究项目";江西省教育厅青年科学基金项目(GJJ12391)

摘　　要：在花岗岩型高放废物地质处置库中,花岗岩中裂隙带的发育程度直接影响到预选区的稳定性和安全性,必须及时有效地探测出裂隙带的发育位置,而地震勘探方法是有效预测和确定裂隙带分布的地球物理勘探方法之一。利用交错网格高阶有限差分算法,采用完全匹配层吸收边界条件(PML),对不同厚度各向同性裂隙带和VTI型裂隙带进行了弹性波场数值模拟,获得裂隙带层状介质的波场快照和地震单炮记录,结果显示,弹性波在花岗岩各向同性和各向异性裂隙带中的传播特征几乎相同,然而随着裂隙带厚度的变化,波场的动力学特征会随之发生改变。In this study a high order finite difference algorithm is developed. The elastic wave numerical simula- tion of isotropic and VTI （transversely isotropic with a vertical symmetry axis） fracture zone is processed with the absorbing boundary condition for perfectly matched layers （PML） under the different thickness. The wave field snapshot and single shot record of fracture zone layered medium is obtained. The results show that elastic wave propagation present no obvious differences between isotropic and Anisotropy fracture zone in granite. The dynamic characteristics of wave field could alter with the changes of fracture zone thickness.

关 键 词：高放废物处置库 有限差分 数值模拟 花岗岩 

分 类 号：P315[天文地球—地震学]