基于ArcGIS与MATLAB的资源量计算方法--以松辽盆地油页岩为例

摘要

基于钻孔资料，采用体积法对松辽盆地青一段、嫩一段、嫩二段油页岩进行了资源量计算，依据分布面积最大的嫩二段建立基准坐标体系，确定了计算范围，采用地理信息系统 ArcGIS对3个层段油页岩面积进行了子块划分与子块面积计算，青一段划分为147个子块，嫩一段划分为197个子块，嫩二段划分为246个子块，子块单位实际面积为537 km2。采用趋势面法，在嫩二段坐标体系下以钻孔地理坐标为自变量，其对应的厚度趋势值为因变量，采用算法软件 MATLAB求得3层油页岩段厚度趋势值方程，并对油页岩厚度空间展布进行拟合，获得了较为理想的拟合度。在确定子块面积与厚度的基础上，对每个子块加权求和分别得出青一段资源量为12724亿 t、嫩一段为8873.1亿 t、嫩二段为9872.4亿 t。%Based on bore hole data,the volumetric method to calculate the resources of oil shale in the 1 st member of Qingshankou formation and the 1st-2nd member of Nenjiang formation in Songliao Basin were adopted.A reasonable area dividing scheme was established taking the 2nd member of Nenjiang formation as criterion when the area distribution of the three sets of oil shale was determined.Sub-block dividing and sub-block area calculation were proceeded on the three researching formations using geographic information system ArcGIS and the actual area of each sub-block was 537 km2 .1st member of Qingshankou formation,1st member of Nenjiang formation and 2nd member of Nenjiang formation were divided into 147 sub-blocks,197 sub-blocks and 246 sub-blocks respectively.According to the theory of trend surface method,oil shale thickness tendency value equation was obtained using algorithm software MATLAB where the spatial distribution coordinates were considered as independent variable and the relative thickness tendency values as dependent variable.Further more,this paper did fitting analysis on oil shale thickness spatial distribution,and the fitting results were satisfied which ensured the accuracy of calculation.The oil shale resources of the three formations were obtained after the determination of the area and thickness of those sub-blocks.The results are 1 272 400 million tons,887 310 million tons and 987 240 million tons for 1st member of Qingshankou formation,1st member of Nenjiang formation and 2nd member of Nenjiang formation respectively.