[发明专利]页岩微米孔隙成像方法及装置

权利要求书

1.页岩微米孔隙成像方法，其特征在于，所述方法包括：

采用同步辐射平行X射线束对页岩进行CT成像，得到所述页岩的成像数据；

根据所述页岩的成像数据建立稀疏正则化模型；

采用梯度下降算法对所述稀疏正则化模型进行求解，得到重构的所述页岩内部的微观图像；

其中，采用梯度下降算法对所述稀疏正则化模型进行求解，得到重构的所述页岩内部的微观图像，包括：

通过以下公式确定所述稀疏正则化模型的梯度；

g(m)≈L^T(Lm-d)+αγ(m)

$\mathrm{\γ}\left(m\right)={(\frac{m_1}{\sqrt{{\left(m_1\right)}^T{m}_1+\mathrm{\ϵ}}},\frac{m_2}{\sqrt{{\left(m_2\right)}^T{m}_2+\mathrm{\ϵ}}},\mathrm{...},\frac{m_n}{\sqrt{{\left(m_n\right)}^T{m}_n+\mathrm{\ϵ}}})}^T$

其中，g(m)表示所述稀疏正则化模型的梯度，L代表离散拉顿变换，m表示重构的所述页岩内部的微观图像，d表示所述页岩的成像数据，L^T表示离散拉顿变换的转置，α表示正则化因子，ε表示大于零的常量，n表示m的维数，m_n表示m的第n个分量，(m_n)^T表示m_n的转置；

通过以下公式确定所述稀疏正则化模型的海森矩阵；

H(m)≈L^TL+αχ₃(m)

其中，H(m)表示所述稀疏正则化模型的海森矩阵，L代表离散拉顿变换，L^T表示离散拉顿变换的转置，α表示正则化因子，m表示重构的所述页岩内部的微观图像，ε表示大于零的常量，n表示m的维数，m_n表示m的第n个分量，(m_n)^T表示m_n的转置，i表示介于1和n之间的常整数；

根据梯度下降算法中的非单调梯度算法通过以下公式确定迭代步长；

${\mathrm{\μ}}_k^{Ray1}={\mathrm{\β}}_1{\mathrm{\μ}}_k^1+{\mathrm{\β}}_2{\mathrm{\μ}}_k^2$

${\mathrm{\μ}}_k^1=\frac{(g_{k-1},g_{k-1})}{(g_{k-1},H_k{g}_{k-1})},{\mathrm{\μ}}_k^2=\frac{(g_{k-1},H_k{g}_{k-1})}{(g_{k-1},H_k^T{g}_{k-1})}$

其中，表示所述迭代步长，k表示迭代次数，β₁和β₂均表示正实数，H表示所述海森矩阵，g表示所述梯度；

按照所述迭代步长和所述梯度采用梯度下降算法对所述稀疏正则化模型进行迭代求解，得到重构的所述页岩内部的微观图像。