岩石弹性模量测试得到的数据应如何进行有效分析和处理？

发布时间：2026-05-17 09:45:00

最近更新：2026-05-17 09:45:00

发布来源：微析技术研究院

本文包含AI生成内容，仅作阅读参考。如需专业数据知识指导，请联系微析在线工程师。

相关服务热线： 156-0036-6678 微析检测业务区域覆盖全国，专注为高分子材料、金属、半导体、汽车、医疗器械等行业提供大型仪器测试、性能测试、成分检测等服务。

岩石弹性模量是反映岩石抵抗弹性变形能力的关键力学参数，直接影响地下工程支护设计、油气井稳定性评价等工程实践的可靠性。然而，由于岩石本身的非均质性、测试仪器的系统误差及操作过程的人为影响，测试数据往往存在离散性与异常值。如何对这些数据进行有效分析与处理，提取真实、可靠的弹性模量特征值，成为岩石力学试验中的核心问题。本文结合试验数据处理的实际流程，从预处理、正态性检验、计算方法、重复性验证等多个维度，系统阐述岩石弹性模量测试数据的分析处理逻辑。

数据预处理：剔除异常值的核心逻辑

岩石弹性模量测试数据的异常值主要来源于三方面：一是仪器误差，如传感器漂移导致应力或应变读数偏差；二是操作失误，如试样安装不垂直使加载力偏心，或应变片粘贴不牢导致应变数据跳变；三是岩石自身的非均质性，如试样中存在大型裂隙、包裹体，导致局部变形异常。这些异常值会严重干扰后续分析，因此预处理的第一步是识别并剔除它们。

常用的异常值检验方法有三种：格拉布斯检验（Grubbs’ Test）适用于小样本（n<30），计算每个数据的G值（|x_i - 均值|/标准差），并与对应显著性水平（如α=0.05）的临界值对比，若G值超过临界值则判定为异常；肖维纳特检验（Chauvenet’s Criterion）通过计算允许的误差范围（W=1+2/√n），若数据与均值的偏差超过W倍标准差则剔除；箱线图法（Box Plot）是可视化方法，以上四分位数（Q3）+1.5×四分位距（IQR）和下四分位数（Q1）-1.5×IQR为阈值，超出范围的点即为异常值。

例如，某花岗岩试样的弹性模量测试结果为38、40、39、45、41 GPa（n=5），计算均值为40.6 GPa，标准差为2.7 GPa。若采用格拉布斯检验（α=0.05，临界值G=1.672），数据45的G值为(45-40.6)/2.7≈1.63，未超过临界值，暂不剔除；若结果为48 GPa，G值为(48-40.6)/2.7≈2.74，远大于临界值，需剔除。

数据正态性检验：后续分析的基础前提

岩石弹性模量测试数据的统计分析（如显著性检验、方差分析）多基于“数据服从正态分布”的假设。若数据不满足正态性，直接使用参数检验会导致结果偏差。因此，正态性检验是后续分析的必要前提。

常用的正态性检验方法包括：Shapiro-Wilk检验适用于小样本（n≤50），通过计算W统计量（范围0-1），若W值接近1且p值大于0.05，说明数据符合正态分布；Kolmogorov-Smirnov检验（K-S检验）适用于大样本（n>50），比较样本累积分布与理论正态分布的差异，p值大于0.05则接受正态假设；Q-Q图（Quantile-Quantile Plot）是可视化工具，将样本分位数与理论正态分位数绘制在坐标系中，若数据点大致沿45度线分布，说明正态性良好。

例如，某砂岩的10个试样弹性模量结果为28、30、29、31、32、30、27、33、30、29 GPa，Shapiro-Wilk检验的W值为0.96，p值为0.85（大于0.05），说明数据符合正态分布；若某泥岩的测试结果为10、12、15、8、20 GPa，Q-Q图显示点偏离45度线，需进行对数转换（转换后结果为2.30、2.48、2.71、2.08、3.00），再检验时p值为0.72，满足正态性要求。

弹性模量计算：从原始曲线到特征值提取

岩石的应力-应变曲线通常分为四个阶段：压密阶段（初始非线性，裂隙闭合）、弹性阶段（线性，应力与应变成正比）、塑性阶段（非线性，微裂隙扩展）、破坏阶段（应力骤降）。弹性模量的计算核心是识别弹性阶段，并提取该阶段的应力-应变斜率。

切线法是最常用的方法，在弹性阶段选取两个点（如应力为0.3σ_c和0.7σ_c，σ_c为岩石单轴抗压强度），计算斜率Δσ/Δε作为弹性模量。该方法适用于弹性阶段明显的硬岩（如花岗岩）。割线法针对压密阶段较长的岩石（如砂岩），取原点到某一应力点（通常为0.5σ_c）的连线斜率，能综合反映压密与弹性变形的共同作用。卸载-再加载法适用于软岩（如泥岩），加载至某一应力（如0.4σ_c）后卸载至0，再加载时的曲线斜率即为弹性模量，可消除初始压密的影响。

例如，某花岗岩的σ_c为150 MPa，弹性阶段在45-105 MPa之间，取应力50 MPa（应变0.0002）和100 MPa（应变0.0004），切线法计算的弹性模量为(100-50)/(0.0004-0.0002)=250000 MPa=250 GPa。某砂岩的σ_c为80 MPa，压密阶段至40 MPa，取原点到40 MPa（应变0.0005）的割线斜率，弹性模量为40/0.0005=80000 MPa=80 GPa。

数据重复性验证：确保结果可靠性的关键步骤

岩石的非均质性导致同一批次试样的弹性模量结果存在离散性，重复性验证的目的是判断这种离散性是否在合理范围内。常用指标是相对标准偏差（RSD），计算公式为RSD=(标准差/均值)×100%，反映数据的相对离散程度。

根据岩石力学试验规范，若RSD≤10%，说明数据重复性良好，可取均值作为该岩样的弹性模量；若RSD>10%，需查找原因：一是试样制备问题，如试样尺寸不一致（高径比偏差超过5%）或表面不平整；二是测试条件问题，如加载速率波动（规范要求0.5-1 MPa/s）或环境湿度变化；三是岩石本身的不均质性，如试样包含不同层理或裂隙。

例如，某石灰岩的5个试样弹性模量结果为45、47、46、44、48 GPa，均值为46 GPa，标准差为1.6 GPa，RSD=(1.6/46)×100%≈3.5%，重复性良好；若结果为45、50、40、55、40 GPa，均值为46 GPa，标准差为6.7 GPa，RSD≈14.6%，需检查试样：发现其中两个试样包含垂直裂隙，导致弹性模量偏低，剔除后重新计算，结果为45、47、46、48 GPa，RSD降至2.8%。

影响因素关联分析：识别数据波动的来源

岩石弹性模量的波动受多种因素影响，需通过关联分析识别主要影响因素，为结果修正提供依据。常见影响因素包括：矿物成分（石英、长石等硬矿物含量越高，弹性模量越大）、结构特征（层理、裂隙发育程度越高，弹性模量越小）、含水率（水的软化作用使弹性模量降低）、测试条件（加载速率越快，弹性模量越高）。

相关性分析是常用方法，通过计算皮尔逊相关系数（Pearson Correlation Coefficient）判断变量间的线性关系：系数绝对值>0.7为强相关，0.3-0.7为中等相关，<0.3为弱相关。例如，某砂岩的石英含量与弹性模量的相关系数为0.82（强正相关），含水率与弹性模量的相关系数为-0.75（强负相关），加载速率与弹性模量的相关系数为0.5（中等正相关）。

回归分析可量化影响程度，建立弹性模量与影响因素的数学模型。例如，针对某泥岩，建立多元线性回归模型：E=0.6×石英含量 - 1.5×含水率 + 2×加载速率 + 8（E单位为GPa，石英含量单位为%，含水率单位为%，加载速率单位为MPa/s）。若某试样的石英含量为50%，含水率为3%，加载速率为1 MPa/s，预测弹性模量为0.6×50 -1.5×3 +2×1 +8=30-4.5+2+8=35.5 GPa，与实际测试结果36 GPa一致。

数据标准化：不同测试条件下的结果可比

不同实验室的测试条件（如试样尺寸、加载设备、环境湿度）差异会导致弹性模量结果不可比，数据标准化的目的是将非标准条件下的结果转换为标准条件下的结果，确保数据的通用性。

试样尺寸标准化是最常见的需求。根据《岩石力学试验方法标准》（GB/T 23561-2009），标准试样为直径50mm、高度100mm（高径比2:1）。若使用非标准尺寸试样（如直径40mm、高度80mm），需引入尺寸效应系数n（通常取0.1-0.2），转换公式为E标准=E测试×(D标准/D测试)^n，其中D为试样直径。例如，D测试=40mm，D标准=50mm，n=0.1，E测试=30 GPa，则E标准=30×(50/40)^0.1≈30×1.023=30.7 GPa。

加载速率标准化需基于回归模型。若某实验室的加载速率为2 MPa/s（非标准），而标准速率为1 MPa/s，根据之前的回归模型E=0.6×石英含量 -1.5×含水率 +2×加载速率 +8，加载速率从2 MPa/s降至1 MPa/s，弹性模量减少2×(2-1)=2 GPa，因此标准结果=测试结果-2 GPa。例如，测试结果为32 GPa，标准结果为30 GPa。

可视化呈现：直观传递数据特征

可视化是数据分析的重要环节，能将抽象的数值转化为直观的图形，快速传递数据特征。针对岩石弹性模量数据，常用的可视化方法包括：

1. 应力-应变曲线对比图：将同一批次试样的应力-应变曲线绘制在同一坐标系中，对比弹性阶段的斜率差异，直观反映试样的均质性。例如，5个花岗岩试样的曲线均较陡且重合度高，说明均质性好；而5个泥岩试样的曲线斜率差异大，说明均质性差。

2. 箱线图：展示数据的中位数、四分位数、异常值，反映数据的离散程度。例如，某砂岩的箱线图显示中位数为30 GPa，上四分位数32 GPa，下四分位数28 GPa，无异常值，说明数据集中；某石灰岩的箱线图显示上四分位数45 GPa，下四分位数35 GPa，异常值55 GPa，说明数据离散。

3. 散点图：展示弹性模量与影响因素的关系，反映变量间的相关性。例如，弹性模量与含水率的散点图显示点从左上到右下，说明含水率越高，弹性模量越小；弹性模量与石英含量的散点图显示点从左下到右上，说明石英含量越高，弹性模量越大。

4. 柱状图：对比不同岩性或不同条件下的弹性模量平均值，反映群体差异。例如，花岗岩（40 GPa）、砂岩（30 GPa）、泥岩（15 GPa）的柱状图，直观展示硬岩与软岩的弹性模量差异；不同含水率（0%、2%、5%）下的柱状图，展示含水率对弹性模量的影响。