大兴安岭北端洛古河东花岗岩的地球化学、SHRIMP锆石U-Pb年龄和岩石成因

大兴安岭北端漠河县洛古河东岩体主要岩石类型为二长花岗斑岩、正长花岗斑岩和石英二长斑岩,内部可见闪长质微粒包体,属高钾钙碱性I型花岗岩。花岗岩的元素地球化学和锆石SHRIMP铀-铅年代学研究结果表明,洛古河东岩体形成于早白垩世,其花岗斑岩体的锆石SHRIMP铀-铅年龄为129.8±2.2Ma。花岗岩的SiO₂含量介于68.03%~74.32%之间,Al₂O₃含量介于13.06%~14.55%之间,Na₂O/ K₂O介于0.45~0.86之间,铝饱和指数为0.94~1.11,Mg^#指数介于18~42之间且多小于30。稀土元素总量为160.00×10^-6 ~ 235.15×10^-6,δEu介于0.31~0.52,(La/Yb)_N介于8.99~17.87,为轻稀土富集型。岩体Sr含量低,介于118×10^-6 ~ 268×10^-6之间,而Y含量高,介于16.9×10^-6~ 26.1×10^-6之间,Sr/Y比值低,介于5.62~13.81之间,属低锶高钇型岩石。在原始地幔标准化的微量元素蛛网图中,Rb、Th、U、K、Zr、Hf和轻稀土元素（如La、Ce、Nd和Sm等）富集,Ba、Sr、P和Ti等元素强烈亏损,Nb和Ta具有中等-弱亏损。主量、稀土和微量元素特征表明,岩石具后碰撞花岗岩类的地球化学特征,属后碰撞花岗岩。岩体ε_Nd (t)值介于-3.45~-2.64,平均-3.01;亏损地幔Nd模式年龄介于969~1131Ma之间,平均1018Ma;锶初始比值 (I_Sr)介于0.702486~0.707269之间,平均0.705434;钾长石²⁰⁶Pb/²⁰⁴Pb、²⁰⁷Pb/²⁰⁴Pb和²⁰⁸Pb/²⁰⁴Pb比值变化范围分别为18.5939~18.6721、15.6019~15.6058和38.4058~38.5249,平均值分别为18.6426、15.6035和38.4613;岩体中的钾长石氧同位素组成很低,δ¹⁸O (‰)值介于-8.1 ~ 4.1之间,多为负值,表明洛古河东岩体为低¹⁸O花岗岩。Nd、Sr、Pb和O同位素组成显示洛古河东岩体形成于含有较多幔源成分的源区物质的部分熔融作用,推测源区主要为Rodinia超大陆会聚过程中（中元古代—新元古代之交）形成的初生地壳。由于古亚洲洋和蒙古-鄂霍茨克洋分别于古生代末期和二叠纪—中侏罗世闭合,因此大兴安岭北端早白垩世花岗岩应该形成于中朝-蒙古大陆与西伯利亚大陆碰撞造山过程的后碰撞阶段。     

辽北西丰地区晚侏罗世营厂岩体LA-ICP-MS锆石U-Pb年龄及岩石地球化学特征

营厂岩体分布于辽北西丰县和隆镇境内，研究结果表明，该岩体主要为似斑状二长花岗岩，锆石²⁰⁶Pb/²³⁸U年龄加权平均值为156.1 Ma±1.5 Ma，形成时代为晚侏罗世。营厂岩体w（Na₂O）/w（K₂O）平均值为1.07，属高钾钙碱性系列；里特曼指数σ平均值为2.15，属钙碱系列；铝饱和指数A/CNK平均值为1.16，A/CN平均值为1.99，属过铝型；Rb、Sr、Ba、Th、Hf含量偏高，Cr、Ni、Co、V、Rb、Zr、Nb等元素含量低，为"S"型花岗岩特征。稀土元素总量平均值为w（ΣREE）＝82.80×10^-6，w（LREE）/w（HREE）=15.76，w（La）_N/w（Yb）_N=24.32，表明其轻稀土相对富集，重稀土相对亏损，稀土配分曲线右倾。δEu=0.80，铕弱负异常，说明成因与板块运动有关，w（Sm）/w（Nd）=0.20，w（Ba）/w（Sr）=1.60，w（K）/w（Rb）=249.15，表明具有壳源性质。营厂岩体具有埃达克岩的特点，其属于受太平洋板块向欧亚板块俯冲作用的影响，在挤压环境下加厚下地壳部分熔融，从而形成活动陆缘型花岗岩。     

塔里木北缘Rodinia超大陆聚合事件——来自新元古代早期花岗岩成因的证据

对塔里木板块拼入Rodinia超大陆聚合时限进行限定,报道塔里木板块北缘库鲁克塔格隆起奥依库如克岩体花岗闪长岩的年代学、岩石学、地球化学和锆石Hf同位素特征。花岗闪长岩的LA-ICP-MS锆石U-Pb年龄为(855±29)Ma,为新元古代早期岩浆作用产物。全岩化学组成上,该岩石为准铝质-弱过铝质、钙碱性系列"I"型花岗岩,轻重稀土元素分馏显著,基本无Eu异常,富集Rb,Ba,K,Sr等大离子亲石元素,亏损P,Ti,Nb,Ta等高场强元素。锆石Hf同位素分析结果显示:锆石εHf(t)值为-15.12～-11.52;二阶段模式年龄TDM2=2.46～2.66 Ga。研究表明,奥依库如克岩体形成于板块俯冲消减环境,指示Rodinia超大陆的汇聚。     

满洲里南部白音高老组流纹岩锆石U-Pb定年及岩石成因

满洲里南部白音高老组火山岩主要由流纹岩组成,含少量珍珠岩和流纹质凝灰岩。LA-ICPMS锆石U-Pb定年结果显示,流纹岩形成于141~139Ma的早白垩世早期。岩石地球化学研究表明,火山岩具有高硅富碱、贫钙镁和高FeO^T/MgO比值的特征; 稀土丰度总量较高(∑REE介于103×10^-6~488×10^-6),轻重稀土分馏明显［(La/ Yb)_N=4.12~30.94)］,Eu负异常显著(δEu=0.12~0.46); 微量元素以富集Rb、Th、U、K,强烈亏损Ba、Sr、P、Ti,中等亏损Nb、Ta和高Ga/Al值为特征,与A-型花岗岩特征相似。锆石¹⁷⁶Hf/¹⁷⁷Hf比值介于0.282785~0.282970之间,ε_Hf(t)值均为正值,介于3.78~9.98之间。流纹岩岩浆来源于斜长石稳定区玄武质下地壳物质的部分熔融,形成于非造山板内伸展构造环境。     

内蒙古边家大院锡多金属矿床成矿岩体岩石地球化学特征及成矿潜力评价

内蒙古边家大院锡多金属矿床是大兴安岭南段多金属成矿带的代表性矿床之一,其西矿区主要发育斑岩型Sn-Cu-Mo矿体,矿体发育于石英斑岩体内。文章通过对石英斑岩开展全岩地球化学、锆石Hf同位素以及锆石微量元素地球化学分析研究,确定了该含矿岩体岩浆性质、来源及演化历史,探讨了成岩成矿关系,并进一步评估了该岩体成锡、铜矿潜力。研究表明,边家大院石英斑岩为准铝质-弱过铝质,高钾钙碱性花岗岩。稀土元素具有轻稀土元素富集、重稀土元素亏损,Eu负异常明显的特征。微量元素具有富集大离子亲石元素,亏损高场强元素的特征。结合其较低的Zr/Hf和Nb/Ta比值以及较高的Rb/Sr比值判断其经历了高分异结晶演化。根据锆石微量元素地球化学特征,确定边家大院石英斑岩源于还原性（ΔFMQ-0.15）、高温贫水（> 750oC）岩浆。边家大院石英斑岩ε_Hf（t）为-0.86~5.99,T_DM2=809~1240 Ma,为中元古代-新元古代年轻地壳部分熔融形成。根据锆石微量元素Ce⁴⁺/Ce³⁺,Ce_N/Ce_N*,Ce/Nd和Eu_N/Eu_N*指针性比值及岩浆特性判断,该石英斑岩有利于锡金属富集成矿,而成铜钼矿潜力小。     

江西冷水坑矿床含矿花岗斑岩的Sr-Nd及锆石Hf-O同位素研究

江西冷水坑铅锌矿床是武夷成矿带北段的典型矿床,矿区内含矿花岗斑岩锆石U-Pb年龄为（162.0±2.0）Ma,岩浆分异程度较高,具有高SiO₂（69.46%~75.52%）、高K₂O（K₂O/Na₂O>3.22）、强过铝质（ASI=1.20~1.96）的特征。岩体发育较明显的Eu负异常,富集强不相容元素Rb、Th、U,而亏损高场强元素Ba、Sr、Nb、Ta、P、Ti。全岩同位素测试结果显示,冷水坑含矿花岗斑岩具有极高的⁸⁷Sr/⁸⁶Sr初始比值（0.74005~0.76518）与低的ε_Nd值（-11.48~-10.78）,锆石Hf-O同位素测定值变化范围较小,¹⁷⁶Hf/¹⁷⁷Hf介于0.282317~0.282460,具有负的ε_Hf（t）值（-12.39~-7.62）和相对较高的δ¹⁸O值（7.23‰~8.81‰）,指示岩浆主要来自于成熟地壳,可能起源于北武夷地区变质基底,源岩为基底变质岩,这可能与古太平洋板块俯冲形成的挤压背景下的地壳熔融有关。     

藏南冈底斯岩基东段晚白垩世早期朗县杂岩的岩石成因和深部地球动力学过程

藏南冈底斯岩基晚白垩世早期岩浆岩保存了有关新特提斯洋的俯冲演化过程和大陆地壳生长的重要信息。本文对朗县杂岩中出露的晚白垩世早期中酸性岩（岩脉）开展了全岩元素地球化学、锆石U-Pb地质年代学、同位素（Sr、Nd和Hf）组成的研究。闪长岩和花岗闪长岩年龄为92.4~86.9Ma,花岗岩（脉）年龄为91.9~88.6Ma,均为晚白垩世早期岩浆作用的产物。闪长岩和花岗闪长岩具有高钾钙碱性偏铝质特征,具有较高的锆石Hf（ε_Hf（t）=+8.3~+13.2,平均值+10.9）和全岩Nd（ε_Nd（t）=+3.2~+2.7）同位素组成,为受俯冲板片流体和大洋沉积物熔体共同交代的地幔楔部分熔融形成,闪长岩岩浆源区形成深度较浅且受板片流体的交代程度较高,花岗闪长岩岩浆源区更深,受沉积物熔体交代程度较高。花岗岩（脉）在主量元素、稀土元素和微量元素组成上显示明显差异,可划分为两类,第一类花岗岩属于低钾钙碱性系列,Na₂O/K₂O>3.0,稀土总量较高,具有明显的Eu负异常,Sr/Y低（<7.2）。在主量元素组成上,该类花岗岩与闪长岩和花岗闪长岩形成较好的线性演化关系,表明它们可能是上述中酸性岩浆演化的产物。第二类花岗岩属于高钾钙碱性系列,Na₂O/K₂O较低（均<1.0）,铝饱和指数较高（A/CNK=1.01~1.02）,发育角闪石,稀土总量较低,具有微弱负或无Eu异常（Eu/Eu^*=0.88~1.12）,Sr/Y比值（33.8~55.4）较高,锆石Hf（ε_Hf（t）=+4.1~+10.8）和全岩Nd（ε_Nd（t）=+0.8）同位素组成都较低,为新生下地壳部分熔融形成。结合冈底斯岩基晚白垩世早期（100~87Ma）岩浆岩已有的研究结果,认为新特提斯洋板片自早白垩世以高角度俯冲,在晚白垩世早期俯冲板片发生回撤,导致软流圈物质上涌增强,诱发熔流体交代过的地幔楔较高程度的部分熔融形成镁铁质岩浆。这些镁铁质岩浆上升侵位到下地壳,发生不同程度的分离结晶作用并诱发新生下地壳部分熔融,形成晚白垩世岩石地球化学性质各异的岩浆岩。     

塔里木盆地北缘二叠纪巴立克立克组凝灰岩年代学、地球化学及构造意义

塔里木盆地北缘广泛分布的二叠纪火山岩,对于理解天山造山带晚古生代构造演化具有重要研究意义。本文以哈拉峻地区巴立克立克组凝灰岩为研究对象,开展了地球化学、年代学和Lu-Hf同位素研究,结果显示凝灰岩富硅(SiO₂=74.59％~76.18％)、富碱(Na₂O+K₂O=7.48％~7.79％)、低铝(Al₂O₃=13.80%~14.29%),属于钾玄质过铝质流纹岩类;具轻稀土元素富集、重稀土元素相对亏损的特征,具Eu弱正异常,相对富集Rb、Th等大离子亲石元素(LILE),亏损P、Nb、Ta、Ti等高场强元素(HFSE),Nb/Ta值平均为16.60,Sr/Y值平均为2.42,Mg^#值平均为34.71;LA-ICP-MS锆石U-Pb年代学分析结果显示巴立克立克组凝灰岩的结晶年龄为280.47±2.58 Ma,εHf(t)值为-0.73~2.12,二阶段模式年龄(t_DM2)为1 318~1 145 Ma。结合塔里木盆地北缘早二叠世的构造演化、岩浆活动及沉积分布等信息,指出巴立克立克组沉积于南天山洋壳向塔里木板块俯冲的板内裂谷盆地中,具非造山板内岩浆岩特征,可能为地幔柱作用背景下,岩石圈地幔底侵导致地壳物质部分熔融的产物,继承了壳源和幔源的双重源区特征。     

苏鲁超高压岩石部分熔融时间的准确限定：来自含黑云母花岗岩中锆石U-Pb定年、REE和Lu-Hf同位素的证据

在北苏鲁超高压变质带中,广泛分布强变形的新三叠纪含黑云母花岗岩和伟晶岩脉。锆石中矿物包体激光拉曼测试、阴极发光图像分析、不同性质锆石微区LA-(MC)-ICP-MS和SHRIMP U-Pb定年、REE及Lu-Hf同位素测试等综合研究结果表明,北苏鲁威海地区超高压正片麻岩在构造折返的 (高压)麻粒岩相升温减压阶段,发生部分熔融 (深熔)作用形成了花岗质岩浆,并在临近角闪岩相退变质作用之前结晶结束形成了含黑云母的花岗岩。该类含黑云母花岗岩中的锆石成因复杂,可划分为三类锆石微区。第一类为强发光效应 (白色)继承性岩浆锆石 (微区I),具有典型的岩浆结晶环带,矿物包体为Qtz+Kfs+Ap,记录的²⁰⁶Pb/²³⁸U年龄为790~782Ma;第二类为新生锆石微区 (微区II),发光强度相对较弱 (灰色-灰白色),也具有较明显的岩浆结晶环带,矿物包体为Qtz+Kfs+Ab+Ap,记录的²⁰⁶Pb/²³⁸U年龄为222~217Ma,加权平均年龄为219±2Ma,表明苏鲁超高压地体的部分熔融 (深熔)作用发生在新三叠纪,这组年龄比苏鲁地体超高压变质时代 (235~225Ma)明显偏新,指示部分熔融 (深熔)作用的时代要晚于苏鲁地体的超高压变质时代;第三类锆石微区 (微区III)围绕第二类锆石微区分布,发光强度最弱 (黑色),也具有典型的岩浆结晶环带,矿物包体十分少见,为Qtz+Ap,记录的²⁰⁶Pb/²³⁸U年龄集中于216~209Ma之间,加权平均年龄为214±2Ma,应代表部分熔融 (深熔)而成的岩浆结晶结束的年龄,这组年龄比苏鲁地体构造折返晚期角闪岩相退变质时代(210~200Ma)偏老,表明新生岩浆结晶结束的时间要早于角闪岩相退变质时代。继承性岩浆结晶锆石 (微区I) ¹⁷⁶Hf/¹⁷⁷Hf(t)=0.281975~0.281984,¹⁷⁶Hf/¹⁷⁷Hf=0.00196~0.00221,ε_Hf(t)=-11.6~-11.8,相应的t_DM2=2160~2170Ma,与研究区周围新元古代 (795~730Ma)正片麻岩继承性岩浆结晶锆石的Lu-Hf同位素特征十分相似,表明新元古代正片麻岩是新三叠纪部分熔融 (深熔)成因的含黑云母花岗岩的母岩。新三叠纪新生岩浆结晶锆石的核部 (微区II)和边部 (微区III)具有类似的¹⁷⁶Hf/¹⁷⁷Hf(t)、¹⁷⁶Hf/¹⁷⁷Hf比值和ε_Hf(t)值,¹⁷⁶Hf/¹⁷⁷Hf(t)=0.282110~0.282199,¹⁷⁶Hf/¹⁷⁷Hf=0.00041~0.00183,ε_Hf(t)=-15.8~-19.1,t_DM2=1980~2130Ma,表明新三叠纪由超高压正片麻岩部分熔融而成的岩浆自结晶开始到结束是在一个相对封闭体系条件下完成的。晚期角闪岩相退变质作用对其Lu-Hf同位素体系也未造成破坏。     

粤西新兴岩体的形成时代与成因研究：对古特提斯洋东支关闭时间的约束

华南板块南缘二叠纪与三叠纪之交的构造属性仍存在较大争议。对新兴岩体详细的野外调查和研究发现,新兴岩体的侵位时间为晚三叠世（240~224 Ma）而非侏罗纪,其主要岩性为细中粒-细粒斑状黑云母二长花岗岩。地球化学分析显示,新兴花岗岩具有高钾钙碱性、过铝质-强过铝质花岗岩特征,具轻稀土元素富集、重稀土元素相对亏损的右倾稀土配分模式,富集大离子亲石元素（Rb、U）而亏损高场强元素（Nb、Ta、Zr、Hf、Ti）。Sr-Nd、Lu-Hf同位素分析显示,新兴花岗岩ε_Nd（t）值介于-11.5~-10.5,ε_Hf（t）值介于-2.9~-10.3,具有壳源源区特征。本次研究表明,粤西地区印支期构造-岩浆活动可能开始于~250 Ma,华南板块南缘海西-印支期岩浆作用自晚二叠世（大容山岩体）一直延续到晚三叠世（新兴岩体）,且晚三叠世仍存在强烈的岩浆活动。新的证据支持古特提斯洋东段分支的关闭时间在~250 Ma,而印支板块和华南板块的陆陆碰撞拼贴一直延续到240~224 Ma。     

Geochemistry of Thermal Waters of Continental Margin of Far East of Russia

Studied waters belong to warm(T=30-50℃),alkaline(pH=8.9-9.3),low mineralized(TDS235 mg/1)Na-HCO_3 or Na-SO_4-HCO_3 thermal waters with high content of SiO_2(up to 81 mg/l)and F(up to 3.9 mg/1),occur on modern volcano-tectonic rejuvenated areas of Eastern Sikhote-Alin orogenic belt.Low~3He concentration as well as N_2/O_2 and N_2/Ar ratios exclude influence of deep mantle fluid.New rare earth element data constrain our understandmg of water-rock interaction occurring in the water source region.Meteoric origin of waters is proved by stable isotope values varying from-71‰to-136.1‰and from-10.8‰to-18.8‰forδ~2U andδ~(18)O respectively.REE patterns reflect high pH,resultfing from water-rock interaction and oxidative conditions.Calculations of deep aquifer temperature using Na-K and quartz geothermometers show 116.8-131.1°C and 82.2-125.8℃respectively.Presence of deep faults both with abnormal thermal gradient(~45-50 K/km)define unique geochemical shape of thermal waters of Sikhote-Alin,area,where no present volcanic activity is registered.     

黄土路基沉降曲线特征分析

在分析郑西客运专线试验段路基沉降实测资料的基础上,提出沉降变形常规预测方法受观测时间及单次观测结果的影响比较大,不利于用来预测最终沉降。通过分析现场沉降变形与时间关系曲线,提出路基沉降的发展主要受固结的影响,其发展是一个从沉降增加到趋缓的过程,可以明显分为固结变形及流变变形,并用非饱和土等效固结理论加以解释,提出了一种利用s- 曲线特征预测最终沉降量的方法,经检验对比发现其预测结果受观测时间及单次观测结果的影响比较小,优于其它常规预测方法     

沉管隧道基槽爆破施工对既有堤岸稳定性影响的数值仿真分析

结合爆破数值模拟中爆破荷载的特点,分析了各种荷载施加方法的特点及爆破震动场模拟的要求,基于圣维南荷载等效原理,提出了爆破荷载的等效施加方法,即将作用在炮孔壁上的爆破荷载等效后,施加在同排炮孔中心线或者面上来模拟爆破荷载的作用。结果表明,等效施加方法在近区存在一定差异,但在中远区吻合得较好。利用该等效荷载施加方法,采用三维实体模型计算分析了某沉管隧道基槽爆破施工对既有堤岸稳定性影响,计算了周围不同位置、不同深度位置处的振动响应,并在此基础上,分析研究了减小爆破振动对既有结构物影响的减振措施,可为类似工程提供借鉴。     

Structure of sediments of kaolinite

In an assembly of clay particles placed in a fluid, each particle is typically subjected to: (1) double-layer repulsive forces; (2) van der Waals attractive forces; and (3) contact mechanical forces. The study presented here concerns an approximate, quantitative analysis of clay suspensions, with considerations to the first two - the physico-chemical forces. Using recent theories to calculate the physico-chemical forces between two clay particles in an approximate model of an assembly, the equilibrium void ratio of a clay suspension is computed. The mechanical forces are ignored in the analysis. The results serve to verify the validity of physico-chemical theories employed and help interpret experimental data more fundamentally in terms of the system variables.     

不同形式诱导缝的等效强度承载性能的数值分析

将考虑尺寸效应规律的混凝土等效断裂模型应用到不同形式诱导缝等效强度公式中,得到了修正后的不同形式诱导缝的等效强度公式。同时,得到了诱导缝等效强度的变化曲线,并比较了单双向间隔诱导缝的优劣,得到双向间隔诱导缝优于单向间隔诱导缝的结论,为高拱坝中诱导缝的设计提出了建议。     

土的可恢复剪胀的一种解释

从土的各向异怀角度对土的可恢复剪胀现象进行了解释。基于各向异性情况下的土体弹性本构关系理论分析，认为土的可恢复剪胀现象可部分归因于土的各向异性引起的弹性剪胀。借助有关土体弹性参数实验结果，研究了应力诱导各迥异性对土体弹性剪胀的影响，结果表明：随土体应力诱导各向异性的增大，土体的弹性剪胀也增大。从土体弹性剪胀角度研究了土的卸荷体缩条件，认为土体卸荷体缩取决于加载应力路径的应力增量比，给出了土体出现卸荷体缩的区域。     

岩石坚硬程度对核电厂地基-基础地震响应特征的影响分析

考虑核电厂地基-基础的动力相互作用,应用显式动力有限差分法分析了地震作用下极软岩、较软岩、坚硬岩上核电厂建筑结构基础的地震响应特征,比较了岩石坚硬程度对基础加速度反应谱的影响。研究表明：随着岩石坚硬程度的提高,核电厂建筑物结构基础的地震响应有增加的趋势;在周期轴上,基础处的加速度反应谱曲线会随着岩石坚硬程度的提高逐渐向短周期（高频段）方向移动。在高频段,建造于较坚硬岩石上基础结构的加速度反应谱值偏大;在中等频段,建造于较软岩石上基础结构的加速度反应谱值偏大;在低频段,岩石坚硬程度对加速度反应谱的影响不显著。     

Strains of scattering of near-field of a point source

Three dimensional scattering of near-field is studied for dilatation and rotation in the time domain. The perturbation method is applied to solve the equation of motion for the first order scattering from a weak inhomogeneity in an otherwise homogeneous medium. The inhomogeneity is assumed close enough to the point source so that the near-field intermediate wave is dominating over the far-field sphericalP andS pulses. The integral expressions are derived to relate dilatation and rotation of scattering to the radial fluctuations of velocities and density in the inhomogeneity. These integrals are solved to calculate the strains of scattering from (a part of) an inhomogeneous spherical shell of arbitrary curvature. Variable curvature may allow the shape of inhomogeneity volume element to change uniformly from spherical to rectangular. Rotation of scattering from a spherical shell is independent ofP wave velocity inhomogeneity. Dilatation of scattering does not involveS wave velocity inhomogeneity but its gradient. The back scattering results are obtained as a special case. Strains are computed numerically, for hypothetical models to study the effects of various parameters viz., velocity inhomogeneity, distance of source from inhomogeneity and from receiver, and thickness of inhomogeneity. The curvature of the spherical shell is varied to study the effects of the shape of inhomogeneous volume element on scattering.     

土石坝心墙水力劈裂分析

简要分析了水力劈裂的发生条件、力学机理。基于水压楔劈机理,利用有限元方法建立了一种水力劈裂的发生判定方法。该方法假定心墙预先存在局部渗透弱面（裂缝）,通过将裂缝位置的单元材料改为裂缝软材料,考虑库水进入裂缝后对裂缝周围土体的作用,建立水力劈裂分析的平面有限元模型,确定裂缝端部垂直于裂缝面的正应力,进而依据该正应力判断水力劈裂发生的可能性,该方法同时可模拟水力劈裂的发展过程。