The grain packing texture of underground sandstone tends to be a major controlling factor in porosity loss during compaction when the effective stress and buried depth increase. However, existing models used to predict porosity loss during sandstone compaction mostly disregard the grain packing texture of sandstone. Based on the micromechanical parameters of discrete element method (DEM) numerical simulation reported in literature, we designed one-component, binary, and ternary packing textures with different grain size distributions and subsequently performed compaction under triaxial servo simulation. We monitored the porosity loss, grain displacement, and force acting on the grain contact point during compaction. Based on the packing texture in sandstone, we propose a new method that considers varying grain sizes, grain size contents, and packing texture types to determine porosity loss during compaction without grain crushing and plastic deformation. The applicability of the method under theoretical conditions was evaluated with 5, 31, and 53 types of one-component, binary, and ternary packing textures. The correlation coefficient between the predicted and simulated values of the void ratio change (△e) was 0.999 and 0.985 for the binary and ternary packing textures, respectively. The reliability of the proposed method was verified using nine groups of physical experimental data derived from literature. The results showed a correlation coefficient of 0.88 between the measured and predicted values of △e and porosity loss (△ϕ). Errors in physical experimental data were derived mainly from grain shape and the crushing of coarse grains. Although the predicted physical experimental data △e and △ϕ were inferior to DEM simulation data, our model was deemed reliable since it showed high correlation between predicted and measured values of △e and △ϕ. Furthermore, the proposed method characterized the influence of the micromechanical process of grain rotation and grain packing texture on reservoir quality during compaction, thereby establishing its importance in predicting reservoir quality of underground sandstone.

当有效应力和埋深增加时，地下砂岩的颗粒堆积结构往往是压实过程中孔隙度损失的主要控制因素。然而，现有的用于预测砂岩压实过程中孔隙度损失的模型大多忽略了砂岩的颗粒堆积结构。基于文献报道的离散元法（DEM）数值模拟的微观力学参数，我们设计了具有不同粒度分布的单组分、二元和三元堆积织构，并随后在三轴伺服模拟下进行压实。我们监测了孔隙率损失、颗粒位移和压实过程中作用在颗粒接触点上的力。基于砂岩中的堆积结构，我们提出了一种考虑不同粒度、粒度含量、和填充纹理类型，以确定在没有颗粒压碎和塑性变形的情况下压实过程中的孔隙率损失。该方法在理论条件下的适用性通过 5、31 和 53 种单组分、二元和三元填充结构进行评估。对于二元和三元填充结构，孔隙比变化 (△e) 的预测值和模拟值之间的相关系数分别为 0.999 和 0.985。使用文献中的九组物理实验数据验证了所提出方法的可靠性。结果表明，△e 和孔隙度损失（△φ）的实测值和预测值之间的相关系数为 0.88。物理实验数据的误差主要来自颗粒形状和粗颗粒的破碎。虽然预测的物理实验数据△e 和△φ 不如DEM 模拟数据，但我们的模型被认为是可靠的，因为它显示△e 和△φ 的预测值和测量值之间的高度相关性。此外，所提出的方法表征了压实过程中颗粒旋转和颗粒堆积结构的微观机械过程对储层质量的影响，从而确立了其在预测地下砂岩储层质量方面的重要性。