基于微观岩石力学试验和NWA13618陨石的小行星岩石力学参数分析

未来人类的地外天体勘探、太空资源开发、地外基地建设等活动都离不开岩土工程技术的发展。目前,人类获得太空岩石样本仅有太空飞行器取样和收集陨石两种途径。由于陨石存量稀少、价格昂贵、尺寸小并且形状任意,因此难以加工成传统力学试验与模拟（mechanical testing&simulation,简称MTS）材料试验机等宏观岩石力学试验所需要的高质量标准岩样。基于微观岩石力学试验和统计分布模型,发展了适用于小尺寸陨石的力学参数测量新技术。首先,应用矿物自动定量分析系统（TESCAN integrated mineral analyzer,简称TIMA）获得NWA13618陨石成岩矿物组成、含量和分布。然后,利用网格纳米压痕技术进行大量压痕试验获得多点弹性模量;之后,利用高斯混合模型,统计求解NWA13618陨石中4种主要矿物的力学参数。橄榄石、辉石、铁镍金属、长石的弹性模量分别为116.73、101.77、87.24、70.74 GPa;最后,基于岩石力学试验结果,采用Mori-Tanaka均质化方法获得NWA13618陨石宏观cm尺度弹性模量为90.48 GPa。提出的微观岩石力学试...     

温度和压强对花岗岩中矿物微观物理力学特性的影响

岩石工程性质取决于岩石中不同矿物的微观特征,本次研究将从微观角度计算北山花岗岩中主要矿物（石英、黑云母和钠长石）的物理性质和力学性质。矿物的微观物理性质使用体积和密度来表征,矿物的微观力学性质使用杨氏模量和泊松比来表征;矿物的初始晶体结构、初始晶胞参数、单晶胞模型和超晶胞模型使用公开数据库获得;矿物晶胞的稳定构象由几何优化方法来实现;使用分子力学模拟,得到了不同矿物稳定构象的微观物理性质参数和微观力学性质参数;使用分子动力学模拟,分析了温度和压强对这些微观物理力学性质参数的影响。结果表明,花岗岩的3种主要矿物中,晶胞体积大小顺序是钠长石、黑云母、石英,晶胞密度大小顺序是黑云母、石英、钠长石;石英的杨氏模量为161.70~168.78GPa、泊松比为0.25,黑云母的杨氏模量为164.85~579.93GPa、泊松比为0.16~0.31,钠长石的杨氏模量为110.72~112.49GPa、泊松比为0.27;压强为0.0001GPa、温度从300K变化到500K时,石英晶胞体积增大0.17%、密度减小0.15%、杨氏模量和泊松比没有明显变化,黑云母晶胞体积增大0.24%、密度减小0.26%、杨氏模量增大140.55%、泊松比减小26.92%,钠长石晶胞体积减小3.76%、密度增大3.91%、杨氏模量增大319.71%、泊松比减小7.41%;温度为298K、压强从0.010GPa变化到0.020GPa时,石英晶胞体积、密度、杨氏模量和泊松比没有明显变化,黑云母晶胞体积和密度没有明显变化、杨氏模量增大53.79%、泊松比减小23.81%,钠长石晶胞体积和密度没有明显变化、杨氏模量减小36.24%、泊松比增大29.41%。研究结果对分析岩石宏观物理力学特性的微观机制具有一定的参考价值。     

水岩循环对变质砂岩微观结构损伤与宏观力学劣化影响机理研究

为研究水岩循环过程中变质砂岩矿物成分及微观结构变化对宏观物理力学参数弱化的影响,通过水岩循环试验,对变质砂岩微观矿物成分、结构在水岩循环过程中的演化机制及其与宏观物理力学参数相互关系开展了研究。研究结果表明,岩石各项物理力学指标随着水岩循环次数的增加而不断衰减,其衰减速率随循环次数增加而减缓。水岩循环后岩石抗压强度试验表明,应力-应变曲线中初始裂隙压密阶段随着循环次数的增加而出现并且变得明显,岩石试样破坏特征有随着循环次数的增加由剪切破坏向拉张破坏转变的趋势。结合XRD和SEM试验手段对矿物成分及微观孔隙平面面积的变化进行了统计,分析认为绿泥石、云母、斜长石等矿物的水化溶蚀作用以及微观孔隙的水化扩展,是造成岩石损伤的主要原因。因此,提出利用矿物成分劣化度D_m和微观结构劣化度D_c的概念来表达岩石水岩循环微观损伤量化指标,并建立了水岩循环作用下受矿物成分及微观结构影响下的岩石物理力学参数弱化的回归模型。该模型的提出能够为水岩循环下相关物理力学参数的选取提供参考。     

黏土矿物组构对水化作用影响评价

泥页岩水化过程中黏土矿物起着主要的作用,但目前还缺少相关的定量描述研究。为此,利用SEM、X-粉晶衍射技术,对采自不同地区的泥岩和板岩在干燥状态、不同浸泡时间后的微观结构和含量进行分析及研究,分析黏土矿物的组构随浸泡时间的变化过程,并定量地分析其对水化过程的影响规律,研究泥页岩的软化过程。研究认为,（1）泥页岩、板岩中影响水化过程的内在因素是蒙脱石、高岭石、绿泥石和伊利石等黏土矿物的含量;（2）黏土矿物定向排列时水化作用效果较无序排列时的明显;（3）泥页岩的水化作用是一个渐进的过程,即在水化过程中流体介质首先使表面的黏土矿物软化,产生掉块;随着时间延长,流体沿裂缝、层理面向岩石内部渗透,黏土土矿物发生膨胀;无围岩限制作用时坍塌掉块持续发生;（4）泥浆中加入一定的无机盐有助于抑制水化过程的发生;压力增加将抑制水化过程的产生;温度升高则能促进泥页岩的水化过程。     

敦煌莫高窟围岩吸水特性及其影响因素分析

为了研究敦煌莫高窟典型围岩吸水特性及其关键影响因素,本文利用自主研发的岩石水理作用智能测试系统,对从敦煌莫高窟采集的3种不同岩性的岩石样品进行了无水压条件下的液态水吸附实验。依据吸水实验获得的3种围岩样品吸水随时间动态变化曲线,分析比较了不同种类岩石样品的吸水特性,建立了其吸水过程函数。同时,通过X射线衍射实验、电子显微镜扫描实验和压汞实验,分别获得了岩石样品的矿物成分及含量,微观孔隙结构特征和孔隙率等重要物理性质参数,综合分析讨论了黏土矿物及微观孔隙结构这两大重要因素对莫高窟3种典型岩石样品吸水特性的影响。研究结果表明,莫高窟3种围岩样品吸水动态变化的总体趋势有相似性,均表现为减速吸水过程,并可以用负指数函数进行拟合。3种岩样吸水能力差异显著,按大小排序依次为:黑砂岩>泥岩>粉砂岩。对于此次研究的3种不同岩性的莫高窟围岩样品而言,黏土矿物（主要是伊利石）的含量不是制约其吸水性的主要因素,决定其吸水量大小、吸水速率快慢的主要因素是岩石样品的微观孔隙结构特征尤其是孔隙率的大小。     

多尺度内聚颗粒模型破碎过程研究

岩石破碎过程涉及到多个变量（应力、应变和孔隙率等）变化和裂纹的产生、拓展和聚集,是研究岩石破碎机制的重要途径。为了克服现有研究中岩石颗粒模型未考虑岩石内部特征问题,针对岩石内部颗粒非均匀分布、聚集的特点,开展了岩石轴压破碎试验和岩石岩相分析试验,并在此基础上,构建符合真实岩石内部特征的多尺度内聚颗粒模型。根据离散元的颗粒黏结模型（BPM）理论,求解了多尺度内聚颗粒模型不同粒级颗粒间黏结键的力学关系,发现与二级颗粒形成的黏结键断裂判据为 ≥ 2 GPa,三级颗粒之间形成的黏结键断裂判据为 ≥ 6 GPa,并基于该判据建立了用于模拟颗粒模型破碎的演化模型。通过模拟轴压破碎试验,破碎演化模型可以从细观角度得到颗粒模型各颗粒间黏结键承受力的实时变化和岩石破碎过程中黏结键从外而内的断裂顺序;从宏观角度得到岩石内部裂纹呈V形从上表面两端延伸并相交于岩石中部。通过与岩石轴压破碎试验结果对比发现,模拟试验得到的岩石裂纹特征与岩石轴压破碎试验结果相似,验证模型的可靠性,实现了从细观和宏观角度分析岩石破碎过程。     

单矿物纳米孔隙特征及其在页岩储集层表征中的意义

为研究单矿物纳米孔隙特征及其在页岩储集层表征中的意义,针对11种单矿物开展了N_2吸附实验,结果表明,高岭石、黄铁矿、钙质及长英质矿物以大孔(50 nm)为主,孔隙形态为楔形;其余黏土矿物介孔(2~50 nm)发育,以平行板状孔隙为主;除铁绿泥石和伊利石主要由微孔(2nm)贡献比表面积之外,单矿物多由介孔贡献;黄铁矿、钙质及长英质矿物中以钠长石孔隙最为发育,黏土矿物中则铁绿泥石高岭石伊利石伊蒙混层钙蒙脱石。同时,基于页岩样品SEM结果分析,当页岩孔隙结构特征与单矿物相近时,通过矿物组成加权构建了页岩吸附等温线及孔径分布。其构建效果受有机孔(或TOC)含量的影响。     

热处理后花岗岩纳米压痕试验研究

随着温度的升高,花岗岩宏观力学参数往往经历略升高、小幅降低、大幅降低几个阶段。利用纳米压痕测试仪器,对热处理后的花岗岩中主要矿物成分进行纳米压痕试验,揭示花岗岩宏观力学参数随温度演化的微观机理,（1）温度在 300 ℃以内微观结构没有显著变化,基本上没有微裂纹产生,在300～500 ℃之间时岩体内部开始产生微裂隙,温度超过500 ℃时岩体内部产生大量微裂隙,并逐渐扩展增大;（2）石英在300 ℃以内弹性模量和硬度略有增加,超过300 ℃弹性模量和硬度开始下降。超过500 ℃晶体结构发生 相到 相转变,弹性模量和硬度急剧下降,长石在400～500 ℃之间时,弹性模量和硬度开始下降,超过800 ℃弹性模量和硬度急剧下降,云母在800 ℃以内弹性模量和硬度有所增加;超过800℃时弹性模量和硬度开始下降;（3）花岗岩宏观力学性能与主要矿物成分的力学性能和微观结构均相关,室温～300 ℃范围内主要受到前者控制,300～500 ℃范围内时受到两者联合控制,超过500 ℃后主要受到后者控制。研究结果可为研究花岗岩和其他岩石的温度效应提供实验数据和理论支持。     

台湾海峡表层沉积物中黏土矿物特征及物质来源

利用X射线衍射仪,对比分析台湾海峡及周边表层沉积物的黏土矿物特征,讨论长江、福建和台湾河流对海峡黏土矿物的贡献程度。福建河流（闽江）的黏土矿物以高岭石为主（>40％）,其次是伊利石（30％）和绿泥石,不含蒙脱石;长江和浙闽泥质区的黏土矿物特征是伊利石含量最高（61％～66％）,其次是绿泥石（16％～17％）和高岭石（10％～13％）,少量蒙脱石（6％～9％）。台湾河流（浊水溪）黏土矿物组合为伊利石（>70％）和绿泥石,不含高岭石和蒙脱石为特征。台湾海峡高岭石含量由西向东明显减少,说明高岭石主要来自福建的河流。〖JP2〗台湾海峡50 m等深线以东区域的黏土矿物具有高的伊利石和绿泥石含量,低的蒙脱石和高岭石含量,且伊利石的化学指数和结晶度值都明显低于长江和福建的河流,反映台湾河流的物质来源。在台湾暖流作用下,台湾河流为台湾海峡中北部提供了大量的伊利石和绿泥石。福建近海的黏土矿物主要是由福建河流提供的,但在平潭岛周边海域的黏土矿物可能有长江物质的混合。〖JP〗     

储层岩石微观孔隙结构对岩石力学特性及裂缝扩展影响研究

储层岩石是一种具有复杂孔隙结构的天然非均质材料,微观结构影响岩石的宏观力学特性和破裂特性,认识储层岩石微观孔隙结构及其对岩石力学特性的影响对油气开采和储层改造具有重要意义。为此,利用CT扫描试验和数字图像处理技术,定量表征储层岩石微观孔隙结构;据此建立微观孔隙结构的三维模型,根据单轴压缩试验获取模型的细观参数,并在PFC3D中进行单轴压缩模拟,分析微观孔隙结构对岩石力学特性和裂缝扩展的影响。结果表明,储层岩石微观孔隙结构概率分布满足对数正态分布,几何形状复杂,孔隙分布具有分形特征;微观孔隙的存在造成单轴抗压强度的降低,分形维数越大,单轴抗压强度越小,两者呈近似线性关系;孔隙结构对裂缝起裂位置、扩展及贯通方向具有决定性作用,裂缝起裂位置多出现在孔隙尖端处;在平行于加载方向的平面内,裂缝多沿孔隙尖端近30°方向和平行于加载方向扩展;在垂直于加载方向的平面内,裂缝沿裂隙轴线方向延伸或使孔洞直径增大;加载过程中裂缝易造成孔隙之间的贯通。     

基于纳米划痕试验的砂岩结构面宏-微观摩擦系数关系研究

基本摩擦系数是影响岩体结构面抗剪强度的重要参数。常规室内试验方法所确定的基本摩擦系数往往受矿物成分、温度等因素影响,为系统揭示结构面的基本摩擦性质,分别研究了砂岩宏观摩擦系数与微观摩擦系数,并建立了两者的关系。首先,通过X射线衍射、纳米压痕试验确定砂岩矿物组分与力学参数。其次,采用倾斜试验、直剪试验开展摩擦系数的宏观尺度研究,对规格为10 cm×10 cm×5 cm的平直结构面试样开展直剪试验,分别施加1、2、3、8、12 MPa的恒定法向应力,研究表明,摩擦系数随法向应力增加呈对数降低趋势,随剪切速率增加呈对数增长趋势。再次,采用纳米划痕试验开展摩擦系数的微观尺度研究,研究表明,低荷载条件下,长石矿物摩擦系数随荷载增加表现出先降低后保持不变的趋势,石英矿物摩擦系数随荷载增加表现出先降低后增加的趋势;此外,随着剪切速率的增加,两种矿物的摩擦系数在低荷载条件下呈增长的趋势,在高荷载条件下基本保持稳定。最后,基于摩擦系数的宏-微观摩擦系数试验结果,采用速度-状态摩擦（RSF）定律建立红砂岩基本摩擦系数与矿物摩擦系数的线性回归方程,并通过直剪试验验证了该经验关系的可靠性,其误差范围为0.17...     

储层孔隙系统“水—岩”反应结晶动力学研究进展

“水—岩”反应作为储层成岩作用研究的重要组成部分,其研究结果对解释储层非均质性成因机制、综合评价储层品质等方面都具有重要意义。近年来随着基础理论、分析测试手段、物理实验方法和模拟技术的进步,该领域研究取得了诸多进展。成岩物理模拟实验的进步实现了对流体—围岩/矿物体系“水—岩”作用过程的宏观尺度观测和研究,也明确了各物理化学参数对体系内矿物溶解—结晶过程的影响,但分辨率限制了其在微观（纳米）尺度解释许多现象成因机制和约束条件方面的应用。晶体生长理论的建立和发展为解释储层孔隙系统中“水—岩”作用过程的结晶动力学原理奠定了基础,特别是近年来伴随纳米科技而发展起来的各种在线和非在线测量技术大大提高了储层孔隙系统结构和内部晶体生长情况的观测精度,将“水—岩”作用研究分辨率提升至纳米级,这为从微观尺度了解孔隙系统中流体—矿物体系的溶解—结晶（沉淀）相平衡过程及其控制因素提供了方案。储层孔隙系统中流体的结晶动力学行为与孔隙介质和流体性质密切相关,是成核自由能、矿物表面电化学特征、传质速率等多因素综合作用的结果。对不同成岩环境和孔隙系统中孔隙流体结晶动力学行为的系统研究,有助于了解孔隙系统空间结构、底衬表面化学特性、表面能效应等对矿物晶体成核生长与溶解过程的影响,为进一步深入解释储层非均质性的结晶动力学原理奠定了理论基础。     

长期注水对致密砂岩油藏岩石力学性质影响机制研究

注水开发是致密砂岩油藏增产的重要手段,长时间注水会改变地层岩石物理力学性质,影响油井产量。为探究长时间注水对地层岩石物理力学性质影响的微观机制,选取同一储层已注水15a及尚未注水岩芯进行对比试验分析。通过试验得到了注水前后地层岩石在弹性力学参数、矿物成分、微观孔隙结构等方面的差异。长时间注水后,岩石矿物成分及内部结构发生变化,主要表现为黏土矿物及方解石含量降低,岩石颗粒间填充物及胶结物大量减少,中小孔隙发育为大孔隙,孔隙度增大,进而导致致密砂岩力学性质弱化,变形能力增大。     

基于吸附等温线的表面分形研究及其地球科学应用

矿物岩石的表面微形貌和孔隙结构是影响其地球化学行为的关键因素，从纳米尺度上表征这一特征对地球化学动力学研究和材料研发有着重要的意义。重点介绍了基于吸附等温线的分形研究方法，以表征纳米尺度上矿物或岩石表面的不规则性和微孔隙结构。从该方法的物理化学原理出发，对比分析了其适用范围和样品限制。在综合当前煤岩学、土壤学、材料学等领域的应用研究成果的基础上，提出了该分形研究方法在地球科学研究中的应用前景和发展趋势。     

使用微观试验确定软土中晶质矿物的类型与含量

土中矿物成分的类型与含量是软土物理力学性质研究的基础,可以使用微观试验进行研究。本次研究以上海软土为例,使用X射线荧光（XRF）分析得到土中不同元素的含量,使用X射线衍射（XRD）分析和XRF分析得到了晶质矿物的类型,根据元素含量和矿物类型得到不同矿物的含量,对不同土层中矿物含量的分布特征进行了分析。结果表明,上海软土主要矿物中的元素为O,Si,Al,Fe,C,K,Mg,Ca和Na;所有土层中,元素含量最多的是O,其次是Si,主要矿物成分为石英、钠长石、方解石、白云石、伊利石、鲕绿泥石和黑云母,石英的含量最高、含量为35.8%~50.8%,其它矿物的含量均不超过21%;黏性土中,元素Al要略多于砂性土或粉性土,伊利石和鲕绿泥石含量明显高于砂性土或粉性土,钠长石和方解石含量要明显小于砂性土或粉性土。研究结果对从微观上研究不同土的物理力学性质有一定的参考价值。     

基于岩石物理模型的碳酸盐岩储层微观孔隙特征分析方法

碳酸盐岩储层的微观孔隙结构普遍具有类型复杂和纵向变化快的特征。由于地球物理资料尺度的限制,现有的解释手段不能有效刻画这些关键特征。笔者采用临界孔隙度模型作为理论依据,建立了微观孔隙结构与宏观弹性参数响应之间的定量关系;以岩石物理实验数据为约束,提出了一种基于概率统计学的方法估算储层关键参数并识别岩石孔隙类型。在YS1井的应用表明：新方法有效挖掘了声波测井数据中隐含的微观信息,参数预测结果达到较高的精度（平均相对误差小于10%）,可以为优质储层段的划分提供可靠的参考依据。     

复杂储层多尺度数字岩石评价

复杂储层岩石矿物组成非均质性强,孔喉结构细小.储集空间有效性评价、岩石结构精细评价及流体赋存状态与运移规律评价是决定复杂储层油气勘探成效的关键.针对复杂储层的储集空间（孔喉、裂缝）、岩石结构（矿物、有机质）、流体特征3方面,建立了复杂储层多尺度数字岩石评价技术及工作流程.储集空间表征方面:二维大面积分析技术可建立跨越6~7个数量级的多尺度选取及非均质性评价;多尺度CT及FIB-SEM联用可精确刻画孔喉和裂缝的三维空间分布;电化学和显影剂技术可以有效地帮助分析微观孔隙连通性.固体组分分析方面:XRF及Qemscan联用可定量评价矿物组成与分布;三维FIB-SEM技术可以实现有机质形态和分布的定量分析.流体特性方面:荷电效应可用于微量残留有机流体的识别与表征;通过合成孔径、润湿性、表面微结构均可调控的纳米材料,开展地层条件下页岩油赋存及流动物理模拟研究,确定了单一因素对页岩油赋存及可动孔径下限的影响;利用分子模拟研究油气在无机、有机质纳米孔隙中的聚集机理与扩散潜力.复杂储层多尺度数字岩石评价技术体系和一系列具体应用可以有效地填补常规储层分析手段的不足,为页岩油气、致密砂岩油气储层以及深部油气储层等复杂储层有效性评价和含油气性定量评价提供技术支撑.      

脆性颗粒材料的动态多尺度模型研究

脆性颗粒材料的多尺度模型一般包含微观尺度的基本粒子、细观尺度的颗粒和宏观尺度的颗粒堆积体3个尺度。基于离散元方法（DEM）构建多尺度模型,并将该模型应用于动态加载。首先,对多尺度模型所涉及的两种接触模型和两种黏结模型的参数进行分析,详细讨论微细观模型参数与宏观材料常数之间的联系。然后,选用Hertz-Mindlin接触模型[1]和平行键黏结模型,建造石英砂的动态多尺度模型。通过选择合适的强度和局部阻尼参数发现,模型宏细观尺度上的动态压缩响应与对石英砂的相关试验结果吻合很好。利用多尺度模型和选定的参数,探讨与动态加载密切相关的局部阻尼机制对多尺度模型各个尺度上力学响应的影响。结果表明,阻尼越大则颗粒材料对波的衰减能力越强,但过高的阻尼会使团簇强度和模型的宏观压缩曲线都表现出异常的加载速度效应（后者实际是阻尼引起的微惯性效应）。另外,高阻尼会过度衰减颗粒破碎过程产生的应力波,从而阻碍颗粒破碎。最后,应用改进的动态多尺度模型,对脆性颗粒材料的动态破碎特性进行研究,发现该模型不但能给出与试验相吻合的颗粒级配曲线,还能揭示出颗粒破碎过程中微裂纹分布的空间不均匀性,即颗粒破碎过程中波的产生机制和衰减机制相互作用导致的微裂纹聚团分布的现象。     

冻融红砂岩损伤演化多尺度分析

以红砂岩为研究对象,进行冻融循环、CT扫描及力学特性试验,采用图像处理技术结合遗传算法寻优模型实现了0、5、10、20、40 次冻融循环后 CT 扫描图像的去噪、增强、分割及三维重构处理,通过对同一对象跨尺度的损伤识别与对比研究,建立了基于细观损伤的弹性模量劣化预测公式,并从材料细观结构的物理本质诠释了冻融红砂岩宏观力学行为。结果表明：基于图像最大熵值的遗传算法能够快速精确地选取阈值进行图像分割,实现对岩石细观结构中基质和缺陷的识别;随着冻融次数增加,岩石孔隙率上升、孔隙分维下降,细观尺度上呈现出孔隙扩展、数量增多,但结构复杂程度下降的演化行为;传统方法以有效承载面积、弹性模量为度量基准定义的宏、细观损伤变量未能全面考虑损伤物理机制和材料内部结构信息,宏细观损伤演化曲线差异较大;基于2种物理机制定义细观损伤变量和考虑岩石天然损伤定义宏观损伤变量,实现了损伤的宏-细观结合。最后通过冻融循环过程中细观结构演化与宏观力学响应之间的关系,提出了弹性模量劣化预测公式,并分析冻融砂岩孔隙大小及孔隙结构形态变化在损伤过程中占据的不同主导作用,根据细观结构的物理机制解释宏观砂岩冻融破坏的力学机制。     

饱和细粒土固结过程中矿物组分与有机质的三维表征

饱和细粒土颗粒细小, 比表面积大, 易与土体中的孔隙水、有机质相互结合形成团聚体, 这种结合形式对土体的物理力学性质有着重要影响。本文从三维微观结构角度探讨饱和细粒土的多物质结构特征对其固结演化性质的作用。土样分成4种组分:孔隙, 有机质, 多矿物组分和伊利石为主的矿物组分。通过同步辐射显微CT技术联合数据约束模型, 获取多物质三维结构表征。伊利石为主的矿物组分, 与孔隙呈共存结构, 即使在1600kPa高压作用下仍不易排出孔隙水。有机质与伊利石为主的矿物组分易结合形成有机黏土颗粒, 有机质不仅吸附在组分表面, 而且嵌入组分结构内部; 在土体固结过程中, 有机质存在富集特征, 40~400μm直径的有机质总体积逐渐增大, 削弱土颗粒间的连结强度。本研究采用定性与定量分析相结合, 获取了固结过程中矿物组分与有机质的结构表征, 为特殊土体的微宏观性质分析提供一种可行性研究方法。