基于MATLAB编程的碎屑岩粒度计算对比:以江西省会昌盆地上白垩统周田组为例  

Comparison of Particle Size Calculation of Clastic Rocks Based on MATLAB Programming:A Case Study of the Upper Cretaceous Zhoutian Formation in the Huichang Basin,Jiangxi Province,China

作　　者：陈念楠 李满根[1,2] 宋志杰 刘东兴 范鹏飞 吴思楷 魏广富 刘颖[1] CHEN Nian-nan;LI Man-gen;SONG Zhi-jie;LIU Dong-xing;FAN Peng-fei;WU Si-kai;WEI Guang-fu;LIU Ying(School of Earth Sciences,East China University of Technology,Nanchang 330013,China;State Key Laboratory of Nuclear Resources and Environment,East China of Technology,Nanchang 330013,China)

机构地区：[1]东华理工大学地球科学学院,南昌330013 [2]东华理工大学,核资源与环境国家重点实验室,南昌330013

出　　处：《科学技术与工程》2024年第11期4409-4415,共7页Science Technology and Engineering

基　　金：中国铀业有限公司-东华理工大学核资源与环境国家重点实验室联合创新基金(2022NRE-LH-12);中国核工业地质局(地LCEQ01)。

摘　　要：计算机强大的图像分析和统计功能在碎屑岩粒度测量过程中能够有效解决计算结果精度不足以及工作负荷大等问题。通过MATLAB编程语言对会昌盆地上白垩统周田组样品进行灰度化、二值化、自适应中值滤波、彩色标注等图像处理方式,准确测量出其粒度参数。结果显示:①MATLAB直接计算法与图解法所求得平均粒径分别为4.008φ和3.103φ(φ为粒级);标准偏差分别为0.8φ和0.803φ;偏度分别为-0.195φ和-0.041φ;峰度分别为0.911φ和0.916φ。两种计算方式中标准偏差和峰度误差较小,可以互相替代,而标准偏差和偏度误差程度较大,不可互相替代;②MATLAB直接计算方法与镜下岩矿鉴定均显示周田组整体磨圆度较差,碎屑颗粒呈棱角状-次棱角状,通过沉积环境判别公式反映周田组形成于河流相沉积环境中;③通过MATLAB编程语言,能够精准识别碎屑岩中颗粒轮廓及统计各项参数,对沉积物粒度研究提供了极大的便利,避免了人为主观因素干扰,提高了工作效率。The powerful image analysis and statistical capabilities of computers can effectively solve problems such as insufficient precision in clastic rock particle size measurements and the heavy workload during the process.The MATLAB programming language was employed to process the samples of the Zhoutian Formation in the Upper Cretaceous of the Huichang Basin through techniques such as grayscale transformation,binarization,adaptive median filtering,and color annotation,accurately measuring their grain size parameters.The results indicate the following.The average grain size obtained using MATLAB's direct calculation method and the graphical method are 4.008φand 3.103φ,respectively.The standard deviations are 0.8φand 0.803φ,the skewness values are-0.195φand-0.041φ,and the kurtosis values are 0.911φand 0.916φ.The standard deviation and kurtosis errors are small,allowing for mutual substitution between the two calculation methods,while the standard deviation and skewness errors are relatively large,making them unsuitable for mutual substitution.Both MATLAB's direct calculation method and microscopic rock mineral identification indicate a relatively poor roundness of the Zhoutian Formation as a whole,with angular to sub-angular clastic particles.The sedimentary environment discrimination formula reflects the formation of the Zhoutian Formation in a fluvial depositional environment,consistent with the geological background and microscopic observation results.Through the MATLAB programming language,it is possible to accurately identify particle outlines and calculate various parameters in clastic rocks,providing great convenience for sediment grain size research,eliminating subjective interference,and improving work efficiency.

关 键 词：粒度分析 图解法 会昌盆地 周田组 MATLAB 

分 类 号：P575.9[天文地球—矿物学] P588.21[天文地球—地质学]