静力触探数据智能预处理（3）

前言

将钻探、物探、静力触探三种数据放在一起对比分析，三种方法间存在尚不明确的物性联系。利用某规范进行土的分类，得出了全是砂土的错误结论，非专业编写，仅供参考。

1、对比结果

三种数据对比结果如下

2、土的分类

某规范内容

按照上述规程进行编程计算，代码片如下

u0 = 0.01*(deepth+18.5)*1.05*9.79363; % 单位为0.1mPa
Qt = qc + (1 - 0.75)*ut*10; % 孔压修正锥尖阻力,单位mPa
P = 2*10^3*deepth*9.79363*10^-6; % 土的总自重压力，单位mPa
Bq = (10*ut - 10*u0)./(qc - P);  %孔隙压力参数比
% Bq(Bq<0) = 0.4; % 任意设的一个阈值
Bqpre1 = 0.2321*qc + 0.06;
Bqpre2 = 0.0714*(qc - 1);

name = cell(length(Qt),1);
for i = 1:length(Qt)
    if qc(i) <= 0.8
        name(i) = cellstr('软土');
    elseif Bq(i) > Bqpre1
        name(i) = cellstr('黏土');
    elseif (Bq(i) <= Bqpre1 ) & (Bq(i) >= Bqpre2)
        name(i) = cellstr('粉质黏土');
    elseif Bq(i) < Bqpre2 & Bq(i) > 0.02
        name(i) = cellstr('粉土');
    else
        name(i) = cellstr('砂土');
    end
end

计算结果显示全部是砂土，令人匪夷所思。