多元储量计算积分法的多元类算法.docxVIP

多元储量计算积分法的多元类算法 1 新型数学方法的算法探索 目前，各种存储算法都是算法的类似方法，因此有些明显的不足（例如，计算过程中的误差较大或过程复杂）很难克服。主要原因是储量计算块段通常具有扭曲特性，现代数学没有相应的精确算法。为此,笔者经多年实践与探索,有幸发现了简柱体新型快捷算法、降维变换、箱变坐标系、变值函数等新型数学方法,从而探求了一种新型储量计算方法或称CS法。 CS法的算法分类可归纳如下,见图1。 其中,一元法来自简柱体新型算法,多元法来自箱变坐标系等新型数学方法。一元非线性类算法只需将一元线性类算法适当扩展即可,为便于交流本文拟对多元类算法的部分要点作一简要介绍或提示。由于箱变坐标系等新型数学方法要点也一并作一简要说明。 2 箱变量坐标系部分要点 由于箱变坐标系是以正箱体为基础和原型,因此,要了解箱变坐标系应首先了解正箱体。 (1) 正实际模型类型 所谓正箱体是指有四个侧面同垂(或正交)于某一平面的六面空间体(如图2)。这类形体为常见的几何形体,正箱体的四个侧面可两两称为横侧面和纵侧面,其余两面叫顶面和底面。其中,四个侧面可为平面或直线型曲面(如柱面等),顶底面可为平面、曲面或扭面。这里的扭面可由两条不在同一平面上的导线和一条母线所决定,一般来说,当导线和母线均为直线时则为线性扭面,当导线和母线不全为直线时既可为非线性扭面也可为线性扭面(如图3)。 正箱体可有八种基本类型(如图3),图中所示均为标准正箱体[即有三个共点侧面(含顶面或底面)相互垂直,否则可通过降维变换而成],其中,当正箱体在基面上的投影为矩形时可叫方箱体,当正箱体为等高方箱体时便为常规箱体(即长方体)。正箱体的数学特征可由其长距、宽距和高距变化特征完全确定,并将其高距变化视为沿长距和宽距方向具有某种规律的连续变化。其中,长距和宽距的变化曲线也可全为外凸或外凹,高距的变化曲线也可全为上凸或上凹。 上述的正箱体类型对应于常规储量计算中的各类单元层。由于单元层的划分通常由平行断面(如常规的断面法)或平行工程(或分段平行,如常规的投影法)所控制,加之相邻工程间的矿体界线(含品位曲线)常按直线或近于直线连接,故其形态及其正厚度(如水平厚度或铅垂厚度等)和线品位的数学图象常为简柱体(用于一元法)或线性直箱体(用于多元法)。在CS法的各类基本算法中,正厚度和线品位是基础,据此首先建立各侧面的相应函数关系,进而求得单元层的各类基本算式。为了突出正厚度和线品位的变化特征,CS法的“线性类”特指二者在四个侧面上均为线性变化,否则称为“非线性类”。 (2) 基型变值及变值系数 箱变坐标系仍由坐标轴、坐标原点和坐标单位所构成,但其坐标单位却处处可变(这是其区别于常规坐标系的本质所在)。其坐标网是以标准正箱体为基准或原型(基箱体),按照等比对应或变值对应进行网状分割而成(如图4中的虚线),分割后的平面图形常为直边或曲边的梯形网,空间图形常为相应的正箱网,这里不象常规坐标系那样只能是方形网,其中,与基箱体对应的坐标网可叫基准网。在箱体坐标系中,基箱体的三个互垂侧面可叫基准侧面或基侧面,常与三个坐标平面相对应,另外的三个侧面可叫变值侧面或变侧面,常用于确定变值系数。当基箱体的三个基侧面与三个坐标平面相互重合时,该箱变坐标系称为最简箱变坐标系,此时三个坐标轴上的长距、宽距和高距分别称为基长、基宽和基高。 在箱变坐标系中,不同位置的同名坐标单位是可变的(叫变值单位),其对应坐标值叫变值坐标,相应坐标系属变值坐标系;而常规坐标系中的同名坐标单位处处相等,其对应坐标值可叫等值坐标,相应坐标系可叫等值坐标系。变值坐标的变值单位与某一基准单位(可叫基值单位,常为同名坐标轴上的坐标单位,其对应坐标值叫基值坐标)之比叫坐标系数或变值系数(相当于坐标自身的变化率),即:变值系数=变值单位/基值单位=变值坐标/基值坐标,或:变值坐标=基值坐标×变值系数。此式可叫变值公式,是进行坐标变换的算法依据。变值系数可为首项为任意常数的连续型初等函数(CS法常取首项为1的整幂函数),其中,首项(KA)对应于方箱体或等值部分,余项(R)对应于变值部分。当变值系数为1时,则变值坐标=基值坐标=等值坐标,当变值系数均为1时,则变值坐标系变为等值坐标系(此时的基箱体为直角六面体)。这里,变值系数的作用是将各类正箱体的坐标单位数变为与之相应的方箱体对应相同,或是将密度均匀的正箱体变为密度不均匀的方箱体,进而按照方箱体进行各种数学运算。 箱变坐标系通常包含着等值坐标系和三种坐标,相应地,求变值坐标的函数叫变值函数,求基值坐标的函数叫基值函数,求等值坐标的函数可叫等值函数(此即常规的数学函数)。三种函数均可由基值坐标或变值坐标来表示,即函数的自变量均可为基值坐标或变值坐标,只不过等值函数的基值坐标与其变值坐标相同而