发明名称 矿井突水水源综合快速判识设备及方法
摘要 本发明公开了一种矿井突水水源综合快速判识设备及方法,包括水位传感器、分光光路模块、水温传感器、模数转换器、处理器和显示器,水位传感器、分光光路模块和水温传感器的信号输出端通过模数转换器连接处理器的信号输入端,处理器的信号输出端连接显示器。本发明将传统的水化学分析判别法、水位分析判别法和水温分析判别法利用上述步骤中记载的方法复合为一个矿井突水水源综合判别模型,使得上述三个传统的方法能够相互弥补各自存在的缺陷,增加矿井突水水源综合判别结果的可信度,降低矿井突水水源综合判别过程的出错率,从而得到更精确的矿井突水水源综合判别结果。
申请公布号 CN103592420B 申请公布日期 2015.06.03
申请号 CN201310571052.X 申请日期 2013.11.15
申请人 武汉长盛煤安科技有限公司;中国矿业大学(北京) 发明人 武强;刘春生;胡加斌;刘鎏
分类号 G01N33/18(2006.01)I 主分类号 G01N33/18(2006.01)I
代理机构 武汉开元知识产权代理有限公司 42104 代理人 潘杰;李满
主权项 一种利用矿井突水水源综合快速判识设备进行矿井突水水源综合快速判识的方法,所述矿井突水水源综合快速判识设备,它包括水位传感器(1)、分光光路模块(3)、水温传感器(4)、模数转换器(8)、中央处理器(9)和显示器(2),其中,所述水位传感器(1)、分光光路模块(3)和水温传感器(4)的信号输出端通过模数转换器(8)连接中央处理器(9)的信号输入端,中央处理器(9)的信号输出端连接显示器(2);它还包括PH值电极(5),所述PH值电极(5)的信号输出端通过模数转换器(8)连接中央处理器(9)的信号输入端;它还包括导电率电极(6),所述导电率电极(6)的信号输出端通过模数转换器(8)连接中央处理器(9)的信号输入端;它还包括氧化还原电位电极(7),所述氧化还原电位电极(7)的信号输出端通过模数转换器(8)连接中央处理器(9)的信号输入端;其特征在于,它包括如下步骤:步骤1:在中央处理器(9)内存储P个参考水源样本,其中,P=m*Z,m代表参考水源样本中含水层种类的个数,Z代表参考水源样本中第i种含水层包含的参考样品个数,i=1,2,…,m;步骤2:在中央处理器(9)内根据P个参考水源样本建立含水层种类的评价集合U,有U={a<sub>1</sub>、a<sub>2</sub>、a<sub>3</sub>...,a<sub>i</sub>,...a<sub>m</sub>},m代表参考水源样本中含水层种类的个数,a<sub>i</sub>代表第i种含水层的判别元素因素集合,其中,所述判别元素因素集合包括含水层中的K<sup>+</sup>离子浓度加上Na<sup>+</sup>离子浓度、Ca<sup>2+</sup>离子浓度、Mg<sup>2+</sup>离子浓度、Cl<sup>‑</sup>离子浓度、SO<sub>4</sub><sup>2‑</sup>离子浓度、HCO<sub>3</sub><sup>‑</sup>离子浓度、含水层的水位和含水层的水温,因此a<sub>i</sub>={a<sub>1i</sub>、a<sub>2i</sub>、a<sub>3i</sub>...,a<sub>ki</sub>,...a<sub>8i</sub>}其中a<sub>8i</sub>代表第i种含水层中第8个判别元素因素的统计平均值;步骤3:对某个待测矿井突水水源进行取样,通过分光光路模块(3)测量所述待测矿井突水水源样本中K<sup>+</sup>离子、Na<sup>+</sup>离子、Ca<sup>2+</sup>离子、Mg<sup>2+</sup>离子、Cl<sup>‑</sup>离子、SO<sub>4</sub><sup>2‑</sup>离子和HCO<sub>3</sub><sup>‑</sup>离子的浓度,通过水位传感器(1)测量所述待测矿井突水水源样本中的水位,通过水温传感器(4)测量所述待测矿井突水水源样本中的水温,并将上述所有的浓度、水温和水位数据发送给中央处理器(9);其中,对于所述待测矿井突水水源样本,K<sup>+</sup>离子的浓度加Na<sup>+</sup>离子的浓度记为y<sub>1</sub>,Ca<sup>2+</sup>离子的浓度记为y<sub>2</sub>,Mg<sup>2+</sup>离子的浓度记为y<sub>3</sub>,Cl<sup>‑</sup>离子的浓度记为y<sub>4</sub>,SO<sub>4</sub><sup>2‑</sup>离子的浓度记为y<sub>5</sub>,HCO<sub>3</sub><sup>‑</sup>离子的浓度记为y<sub>6</sub>,水位记为y<sub>7</sub>,水温记为y<sub>8</sub>,则得到待测矿井突水水源的判别元素因素实测值集合y={y<sub>1</sub>,y<sub>2</sub>,y<sub>3</sub>,y<sub>4</sub>,y<sub>5</sub>,y<sub>6</sub>,y<sub>7</sub>,y<sub>8</sub>};步骤4:由于上述各个待测矿井突水水源的判别元素因素,对待测矿井突水水源属于那种参考水源判断时的重要性不同,首先在中央处理器(9)中建立待测矿井突水水源相对于参考水源的判别因子权重模糊矩阵E={e<sub>k1</sub>,e<sub>k2</sub>,…e<sub>km</sub>},k=1,2,…,8;其中,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>e</mi><mi>ki</mi></msub><mo>=</mo><mfrac><msub><mi>y</mi><mi>k</mi></msub><msub><mi>a</mi><mi>ki</mi></msub></mfrac><mo>/</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>8</mn></munderover><mfrac><msub><mi>y</mi><mi>k</mi></msub><msub><mi>a</mi><mi>ki</mi></msub></mfrac></mtd></mtr><mtr><mtd><msub><mi>e</mi><mi>ki</mi></msub><mo>&GreaterEqual;</mo><mn>0</mn></mtd></mtr><mtr><mtd><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mn>8</mn></munderover><msub><mi>e</mi><mi>ki</mi></msub><mo>=</mo><mn>1</mn></mtd></mtr><mtr><mtd><msub><mi>a</mi><mi>ki</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>Z</mi></munderover><msub><mi>x</mi><mi>kip</mi></msub><mo>/</mo><mi>Z</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000624109680000021.GIF" wi="838" he="1223" /></maths>上述公式(1)中y<sub>k</sub>为所述待测矿井突水水源样本的第k个因素的实测值k=1,2,…,8,a<sub>ki</sub>为参考水源样本中第i种含水层的第k个因素的统计平均值,其中Z代表参考水源样本中第i种含水层包含的参考样品个数,k=1,2,…,8;i=1,2,…,m,x<sub>kip</sub>为参考水源样本的第i种含水层中包含的第p个参考水源样本的第k个评判因素的数值;其次得出待测矿井突水水源相对于参考水源的判别因子权重模糊矩阵E:<img file="FDA0000624109680000031.GIF" wi="885" he="402" />然后在中央处理器(9)内建立待测矿井突水水源样本相对于参考水源样本的隶属度模糊关系矩阵R={r<sub>ki</sub>},r<sub>ki</sub>表示待测矿井突水水源样本第k个评判因素相对于参考水源样本中第i种含水层的隶属度,由下列隶属函数计算:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>R</mi><mrow><mo>(</mo><msub><mi>r</mi><mi>ki</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mn>1</mn><msub><mi>b</mi><mi>ki</mi></msub></mfrac><msup><mrow><mo>(</mo><msub><mi>y</mi><mi>k</mi></msub><mo>-</mo><msub><mi>a</mi><mi>ki</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup></mrow>]]></math><img file="FDA0000624109680000032.GIF" wi="615" he="197" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>b</mi><mi>ki</mi></msub><mo>=</mo><msup><mrow><mo>[</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mi>Z</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>kip</mi></msub><mo>-</mo><msub><mi>a</mi><mi>ki</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>Z</mi><mo>-</mo><mn>1</mn></mrow></mfrac><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000624109680000033.GIF" wi="1056" he="528" /></maths>上述公式(3)中b<sub>ki</sub>为参考水源样本中第k个评判因素对第i种含水层的标准差,y<sub>k</sub>为所述待测矿井突水水源样本的第k个因素的实测值k=1,2,…,8,a<sub>ki</sub>为参考水源样本中第i种含水层的第k个因素的统计平均值,x<sub>kip</sub>为参考水源样本的第i种含水层中包含的第p个参考水源样本的第k个评判因素的数值,Z代表参考水源样本中第i种含水层包含的参考样品,于是对每一个待测矿井突水水源样本,均能得到一个模糊关系矩阵:<img file="FDA0000624109680000041.GIF" wi="941" he="452" />其中,T为转置符号;最后在中央处理器(9)内对待测矿井突水水源相对于参考水源的判别因子权重模糊矩阵E的m行8列数据的同一列数据相加求平均值,得到一个1行8列的待测矿井突水水源相对于参考水源的判别因子权重均值矩阵E1,得到待测矿井突水水源样本隶属于各种含水层的概率集合B,B=E1*R<sup>T</sup>,得到B={b<sub>1</sub>,b<sub>2</sub>,…b<sub>m</sub>}<sup>T</sup>,b<sub>i</sub>为该样品隶属于各个含水层的程度,i=1,2,…,m,在待测矿井突水水源样本隶属于各种含水层的概率集合B中概率最大值对应的含水层种类即为待测矿井突水水源对应的含水层种类。
地址 430312 湖北省武汉市盘龙城经济开发区巨龙大道第一企业社区38号