CN106869990B - 基于lvq-cpso-bp算法的煤体瓦斯渗透率预测方法 - Google Patents
基于lvq-cpso-bp算法的煤体瓦斯渗透率预测方法 Download PDFInfo
- Publication number
- CN106869990B CN106869990B CN201710121485.3A CN201710121485A CN106869990B CN 106869990 B CN106869990 B CN 106869990B CN 201710121485 A CN201710121485 A CN 201710121485A CN 106869990 B CN106869990 B CN 106869990B
- Authority
- CN
- China
- Prior art keywords
- lvq
- prediction
- cpso
- algorithm
- layer
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 230000035699 permeability Effects 0.000 title claims abstract description 83
- 238000004422 calculation algorithm Methods 0.000 title claims abstract description 75
- 239000003034 coal gas Substances 0.000 title claims abstract description 69
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000013528 artificial neural network Methods 0.000 claims abstract description 49
- 239000003245 coal Substances 0.000 claims abstract description 40
- 238000012549 training Methods 0.000 claims abstract description 27
- 238000005457 optimization Methods 0.000 claims abstract description 22
- 239000002245 particle Substances 0.000 claims abstract description 18
- 239000013598 vector Substances 0.000 claims abstract description 16
- 238000007906 compression Methods 0.000 claims description 22
- 230000006835 compression Effects 0.000 claims description 19
- 210000002569 neuron Anatomy 0.000 claims description 16
- 230000006870 function Effects 0.000 claims description 9
- 241001269238 Data Species 0.000 claims description 6
- 238000011156 evaluation Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000000739 chaotic effect Effects 0.000 claims description 3
- 239000006185 dispersion Substances 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 241000894007 species Species 0.000 claims description 3
- 230000008595 infiltration Effects 0.000 claims description 2
- 238000001764 infiltration Methods 0.000 claims description 2
- 238000009738 saturating Methods 0.000 claims 1
- 210000003127 knee Anatomy 0.000 abstract description 5
- 238000013139 quantization Methods 0.000 abstract description 4
- 238000010835 comparative analysis Methods 0.000 abstract description 2
- 230000008859 change Effects 0.000 description 7
- 238000012360 testing method Methods 0.000 description 7
- 238000005065 mining Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 244000287680 Garcinia dulcis Species 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000011478 gradient descent method Methods 0.000 description 1
- 230000002687 intercalation Effects 0.000 description 1
- 238000009830 intercalation Methods 0.000 description 1
- 238000004452 microanalysis Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 210000004205 output neuron Anatomy 0.000 description 1
- 230000002028 premature Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F7/00—Methods or devices for drawing- off gases with or without subsequent use of the gas for any purpose
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Physics & Mathematics (AREA)
- Mining & Mineral Resources (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Linguistics (AREA)
- Computer Hardware Design (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Geology (AREA)
- Biomedical Technology (AREA)
- Artificial Intelligence (AREA)
- Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Geometry (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Geochemistry & Mineralogy (AREA)
- Software Systems (AREA)
- General Life Sciences & Earth Sciences (AREA)
Abstract
基于LVQ‑CPSO‑BP算法的煤体瓦斯渗透率预测方法,提出基于学习向量量化神经网络(LVQ)分类、混沌粒子群算法(CPSO)优化、BP神经网络预测的LVQ‑CPSO‑BP煤体瓦斯渗透率预测方法。确定临界值将煤层埋深划分为两层;基于有效应力与瓦斯渗透率之间存在拐点关系,确定拐点值将有效应力划分为两段;采用LVQ将4个微观样本参数依据拐点特征进行分类识别,采用BP神经网络进行学习训练并输出预测结果,并用CPSO对BP神经网络的权值和阈值进行优化;最后,基于样本案例对本发明构建的LVQ‑CPSO‑BP算法进行了预测结果验证,并与BP算法、GA‑BP算法及PSO‑BP算法预测的结果进行了对比分析。
Description
技术领域
本发明涉及煤体瓦斯渗透率预测领域,特别涉及基于LVQ-CPSO-BP算法的煤体瓦斯渗透率预测方法。
背景技术
煤矿瓦斯灾害是煤矿开采过程中的主要灾害之一,受地质条件的影响,使得矿井瓦斯在煤层中多以失稳状态渗流运动。由于瓦斯渗透率与地应力、温度和瓦斯压力等所处的地球时间空间特性息息相关,致使瓦斯渗透率呈现出时变性、非线性、模糊性等特点,如何精确预测煤体瓦斯渗透率的动态变化,对预防矿山开采中煤矿瓦斯灾害的发生有重要意义。
国内外研究人员针对煤体瓦斯渗透率的变化进行了研究,吕闰生等研究了不同结构下煤体瓦斯渗透率的差异性,揭示了煤体结构差异性对瓦斯渗透率的影响;魏建平等研究了水分和受载条件下含瓦斯煤的渗透性影响因素;尹光志,蒋长宝等研究了瓦斯压力、地应力与渗透率的关系;王登科等研究了瓦斯的渗流特性,得出了渗透率新的计算方法;许江等研究了温度、低孔隙压力和不同瓦斯压力的环境下对煤体瓦斯渗透率特性的影响;Perera等研究了煤体瓦斯气体在非零横向应变条件下的渗透率;赵阳升等研究了三维应力下吸附作用对煤岩体气体渗流规律。由于开采过程中煤体瓦斯渗透率的不确定性以及模糊性特征,国内外学者建立的理论模型和经验公式是不适用于煤体瓦斯渗透率预测的。因此,以智能算法为主的预测方法开始得到了应用,其中BP神经网络(back propagation neuralnetwork)是用于煤体瓦斯渗透率预测的主要方法之一,一定程度上解决了矿山安全预测问题。但是以BP神经网络为主的单个预测模型算法只能针对特定问题进行解答,因此以组合算法为主架构的混合预测模型研究得到迅速发展,即考虑单一预测模型的优缺点及适用性,通过构建智能优化算法结合神经网络分析煤体瓦斯渗透率在不同条件下的变化特点,建立一种优化算法与神经网络结合的双组合算法,提高了预测精度,然而,获得更高精度的煤体瓦斯渗透率预测在算法上需要更多的组合和改进。
综上所述,本发明提出采用学习向量量化神经网络算法(Learning VectorQuantization,即LVQ)分类、混沌与粒子群组合生成的混沌粒子群算法(Chaos ParticleSwarm Optimization,CPSO)优化、BP神经网络预测的多组合算法对煤体瓦斯渗透率进行预测,通过煤层埋深宏观分层、有效应力微观分段、权值和阈值并行优化、样本数据学习训练、检验样本预测验证,并与BP单一算法、遗传和BP双组合算法(Genetic AlgorithmGA-BP)、粒子群和BP双组合算法(Particle Swarm Optimization,PSO-BP)预测的结果进行对比分析,充分证明了LVQ-CPSO-BP多组合算法在煤体瓦斯渗透率预测精度提高方面可行、有效。
发明内容
为解决上述现有技术存在的问题,本发明的目的在于提供基于LVQ-CPSO-BP算法的煤体瓦斯渗透率预测方法。总体预测值与实测值吻合度高,尤其当有效应力减小时,预测精度更高。
为达到上述目的,本发明的技术方案为:
基于LVQ-CPSO-BP算法的煤体瓦斯渗透率预测方法,由如下步骤构成:
步骤一、煤层埋深划分及选择:对煤层埋深进行研究:煤层埋深分为浅埋层和非浅埋层,而浅埋层一般在150m以上,处于瓦斯风化带,本发明不做讨论;而非浅埋层是指150m以下,最深达到1500m左右,大部分煤层处于200m-700m之间,因此本发明忽略浅埋层,将非浅埋层自定义分为两层,设临界值为700m,以浅为中埋层,以深为深埋层;
步骤二,主要影响因素选择:选择有效应力、温度、瓦斯压力、抗压强度作为主要影响因素;
有效应力与煤体瓦斯渗透率之间的关系出现拐点,有效应力小于拐点时,煤体瓦斯渗透率变化较快,大于拐点时,煤体瓦斯渗透率变化平缓,将2MPa设为有效应力拐点;
步骤三、在煤层埋深非浅埋层宏观划分基础上,采用LVQ算法对有效应力样本数据按拐点进行微观分类;
步骤四、建立基于CPSO-BP算法的煤体瓦斯渗透率预测子模型;
步骤五、通过LVQ-CPSO-BP预测***的训练,获得煤体瓦斯渗透率预测结果。
进一步的,为了评价所建立预测模型的可信度及精确度,对误差和精度进行量化,采用式(1)-式(4)预测误差和精度评价指标进行煤体瓦斯渗透率预测结果定量评估:
1)预测误差的可靠性用平均相对误差:(mean relative error,MRE):
2)预测值偏离真值误差的大小用平均绝对误差(mean absolute error,MAE):
3)样本数据离散程度用均方根误差(root mean square error,RMSE):
4)预测模型精确度(forecast model accuracy,FMA):
其中:Xi为真实值,X′i为与预测值,M为样本数量。
进一步的,所述步骤三中,采用LVQ以有效应力拐点为界将其分为煤体瓦斯渗透率变化平缓型和陡峭型两种类型,以样本数据中有效应力为(1-5)MPa的平均值温度平均值瓦斯压力平均值抗压强度平均值作为LVQ的输入,经过LVQ分类识别训练后,获得两类分类样本数据。
进一步的,所述步骤三中,煤体瓦斯渗透率预测学习样本向量用X表示,样本集为:
[X]=[x1,x2,x3,x4]T=[σ,T,P,σm]T (5)
其中σ为有效应力(MPa)、T为温度(℃)、P为瓦斯压力(MPa)、σm为抗压强度(MPa);
输入层输入集为:
竞争层的样本用Y表示为:
Y=(y1,y2,…ym)T,yi∈{0,1},i=1,2,…,m (7)
输入层和竞争层之间的权值矩阵用W1:
竞争层到输出层之间的权值矩阵用W2:
网络的期望值输出用d表示为:
d=(d1,d2,…dl)T (10)
输出层的输出用O表示,输出集为:
[O]=[o1,o2]T=[1,0]T (11)
其中,o1为陡峭型类别,o2为平缓型类别。
进一步的,所述步骤四中,在LVQ将样本数据分为两大类后,采用CPSO-BP预测,CPSO算法采用Logistic映射混沌***,迭代公式为:
zi+1=μzi(1-zi),i=0,1,…,μ∈(2,4] (12)
式中:μ为控制变量;
当μ=4,0≤z0≤1时,Logistic完全处于混沌状态,利用μ=4时的混沌运动特性改进标准粒子群算法,即
zi+1=4zi(1-zi),i=0,1,… (13)
粒子群初始化阶段利用公式(13)产生出优良的初始群体作为初始粒子群,以减轻随机初始化粒子群对算法优化性能的影响。
进一步的,所述步骤四中,将BP神经网络结构的层数设为3层,设R、S1、S2分别为输入层、隐含层、输出层神经元个数,其中的输入层神经元为LVQ分类后的样本数据,隐含层为S函数,最佳隐含节点数为3个,输出层为煤体瓦斯渗透率。BP神经网络的权值和阈值与粒子群个***置向量P的维数N之间的映射关系如公式(14)所示:
N=S1×R+S1×S2+S1+S2 (14)
进一步的,所述步骤三中、经过影响因素筛选,以有效应力,温度,瓦斯压力,抗压强度与渗透率的对应数据为样本进行学习;
采用Matlab2015神经网络工具箱结合样本数据编写了BP、GA-BP、PSO-BP、CPSO-BP四种组合方式渗透率预测模型代码,设定BP神经网络迭代次数1500次,学习率0.01,最小均方根误差10-5;CPSO优化神经网络的粒子取W=0.7298,c1=c2=1.4962,粒子的上界值为5,下界值为-5。
进一步的,所述步骤五中、将样本数据输入LVQ进行训练,首先进行网络初始化:设定较小随机数作为初始权值W1和W2;确定初始学***均向量。通过期望值最小寻找获胜神经元。以胜出神经元小于本发明指定的有效应力平均值(约2.5MPa)的划分为正确分类,否则调整分类时的权值,更新学习效率,直到t=tm,获得分类结果。
进一步的,所述步骤五中,以有效应力、温度、瓦斯压力、抗压强度的平均值作为输入样本,规定以有效应力平均值为拐点的临界值左边和右边的输出类型依次为1,0,得到分类识别结果;
将样本数据输入LVQ进行分类后的两组样本数据输入CPSO-BP进行训练,首先初始化BP神经网络3层结构,输入层神经元为4:有效应力、温度、瓦斯压力、抗压强度;隐含层神经元为5,输出层神经元为1:煤体瓦斯渗透率;通过混沌初始化种群大小、维度D及粒子S速度和位置及均方误差粒子群适应度函数自寻优,将寻找出的最优解赋给BP网络的权值及阈值,继续BP网络的学习训练预测。
进一步的,本发明采用四种预测模型进行了40组样本数据的学习训练,预测模型一为BP神经网络预测,预测模型二为GA-BP神经网络预测,预测模型三为PSO-BP神经网络预测,预测模型四为CPSO-BP;预测模型中的输入数据均为LVQ分类后的数据,其中类型1训练样本为26组;类型0训练样本为14组。
进一步的,采用LVQ模型对有效应力拐点划分具有较高的精度,分类识别正确率最高可达96.15%。
进一步的,依据建立的误差评估方法对煤体瓦斯渗透率的四种预测结果与试验真实值之差进行误差计算,LVQ-CPSO-BP算法的煤体瓦斯渗透率预测精度为95.26%,优于其他三种模型,所建立的LVQ-CPSO-BP预测模型有效、适应性更强、预测精度更高。
相对于现有技术,本发明的有益效果为:
本发明基于LVQ-CPSO-BP算法的煤体瓦斯渗透率预测方法,从宏观与微观结合角度,筛选出影响预测精度的5个主要因素—1个宏观(煤层埋深)和4个微观(有效应力、温度、瓦斯压力、抗压强度),提出一种基于学习向量量化神经网络(LVQ)分类、混沌粒子群算法(CPSO)优化、BP神经网络预测的LVQ-CPSO-BP煤体瓦斯渗透率预测方法。首先,基于我国煤层埋深值的统计数据,从宏观上确定临界值将煤层埋深划分为两层;基于有效应力与瓦斯渗透率之间存在拐点关系,从微观上确定拐点值将有效应力划分为两段;接着,采用LVQ将4个微观样本参数依据拐点特征进行分类识别,采用BP神经网络进行学习训练并输出预测结果,并用CPSO对BP神经网络的权值和阈值进行优化;最后,基于样本案例对本发明构建的LVQ-CPSO-BP算法进行了预测结果验证,并与BP算法、GA-BP算法及PSO-BP算法预测的结果进行了对比分析。结果表明:LVQ分类正确识别率较高,CPSO-BP算法预测精度较好,且优于其他三种算法。LVQ-CPSO-BP算法总体预测值与实测值吻合度高,尤其当有效应力减小时,预测精度更高。
附图说明
图1为渗透率随有效应力变化图,其中,(a)拐点为1.5MPa,(b)拐点为3MPa,(c)拐点为2MPa。
图2为煤体瓦斯渗透率预测总体架构框图。
图3 LVQ网络结构图。
图4 CPSO-BP算法流程框图。
图5有效应力拐点分类预测结果图,其中:(a)类型1预测结果(b)类型0预测结果。
具体实施方式
下面结合附图和具体实施方式对本发明技术方案做进一步详细描述:
实验例:
如图1所示,
基于LVQ-CPSO-BP算法的煤体瓦斯渗透率预测方法,由如下步骤构成:
主要影响因素筛选及分层临界值与应力拐点值确定
煤体瓦斯渗透率与地质时间空间特性息息相关,随着煤矿开采过程的变化,呈现出时变性、非线性、模糊性等特点,表现出动态变化特征。其影响因素宏观上主要体现在地质构造、地质煤层结构、煤层开采深度、煤质特征等方面,微观上主要体现在有效应力、温度、瓦斯压力、抗压强度、孔隙结构等方面,忽略以上因素与渗透率之间的关系,则不能真实反映内在联系与规律,全面考虑以上因素又会造成参数繁多、交叉干扰、模型复杂、难以预测。因此对影响因素进行划分与筛选及分类,将是提高煤体瓦斯渗透率预测精度的有效手段。
在宏观方面,从已发表的研究文献可知:煤层埋深变化将产生不同的地应力,对瓦斯渗透率影响较大,而地质构造、煤质特征相对影响较弱,因此本发明只对煤层埋深进行研究:煤层埋深分为浅埋层和非浅埋层,而浅埋层一般在150m以浅,处于瓦斯风化带,本发明不做讨论。而非浅埋层是指150m以深,最深达到1500m左右,大部分煤层处于200m-700m之间,因此本发明忽略浅埋层,将非浅埋层自定义分为两层,设临界值为700m,以浅为中埋层,以深为深埋层。
在微观方面,从已发表的研究文献可知:有效应力、温度、瓦斯压力、抗压强度、孔隙结构影响较大,但是文献指出,孔隙随应力变化而变化,因此选择有效应力、温度、瓦斯压力、抗压强度作为主要影响因素。其中的有效应力根据文献的试验数据具有如图1所示的特征:有效应力与煤体瓦斯渗透率之间的关系出现拐点,有效应力小于拐点时,煤体瓦斯渗透率变化较快,大于拐点时,煤体瓦斯渗透率变化平缓。由于试验条件不同,出现拐点的位置也不同,图1(a)中拐点为1.5MPa,图1(b)中拐点为3MPa,为了体现通用性,本发明提出将两组试验数据中的有效应力拐点的平均值定义为拐点,图1(c)中为2MPa的拐点充分体现了本发明定义拐点值的合理性。
LVQ-CPSO-BP预测***架构
LVQ具有对大样本量进行分类识别的功能,CPSO具有对不确定性及模糊性量的权值和阈值进行智能优化的功能,BP神经网络具有对随机及非线性量进行自学习预测功能。
因此,在煤层埋深非浅埋层宏观划分基础上,采用LVQ算法对有效应力样本数据按拐点进行微观分类,建立基于CPSO-BP算法的煤体瓦斯渗透率预测子模型。通过LVQ-CPSO-BP预测***的训练,获得煤体瓦斯渗透率预测结果。
预测精度评估方法
为了评价所建立预测模型的可信度及精确度,对误差和精度进行量化,采用式(1)-式(4)预测误差和精度评价指标进行煤体瓦斯渗透率预测结果定量评估。
1)预测误差的可靠性用平均相对误差:(mean relative error,MRE):
2)预测值偏离真值误差的大小用平均绝对误差(mean absolute error,MAE):
3)样本数据离散程度用均方根误差(root mean square error,RMSE):
4)预测模型精确度(forecast model accuracy,FMA):
其中:Xi为真实值,X′i为与预测值,M为样本数量。
LVQ-CPSO-BP算法
LVQ算法
LVQ属于多层前馈型神经网络,用于输入样本分类识别。网络由输入层、竞争层和输出层三层组成,当样本数据作为LVQ的输入时,竞争层的神经元通过胜者为王竞争规则产生获胜神经元,输出为1,失败者输出为0。与获胜神经元所在组相连接的输出神经元输出为1,其他输出为0,以达到输出样本分类识别的目的。LVQ网络结构如图3所示。
CPSO-BP算法
BP算法是一种常见的有导师监督型算法,通过真实值与预测值之间的误差,反向传播给输入层、隐含层和输出层之间的权值和阈值,由于BP算法本质上为梯度下降法,运算后期易陷入局部最优、收敛速度慢,而PSO算法是一种迭代的优化工具,具有全局寻优的特点,它在优化神经网络时,无需导数等梯度信息,可很好地解决神经网络的结构和权值优化问题,但是PS0算法存在寻优过程中种群多样性丧失问题,容易导致算法出现早期收敛的现象,因此将混沌思想与PS0算法相结合,应用混沌算法对PSO算法的粒子进行搜索,使得两种算法的优点相融合,形成CPSO算法。采用CPSO算法优化BP神经网络的权值和阈值,既达到了全局寻优的目的,又提高了收敛速度。CPSO-BP算法流程如图4所示。
LVQ-CPSO-BP预测模型及验证
LVQ-CPSO-BP预测模型
LVQ分类识别模型
煤体瓦斯渗透率随着有效应力的变化呈现出两种变化规律,拐点之前变化陡峭,拐点之后变化平缓,显而易见,值大又陡峭的煤体瓦斯渗透率安全隐患更大。因此本发明采用LVQ以有效应力拐点为界将其分为煤体瓦斯渗透率变化平缓型和陡峭型两种类型,以样本数据中有效应力为(1-5)MPa的平均值温度平均值瓦斯压力平均值抗压强度平均值作为LVQ的输入,经过LVQ分类识别训练后,获得两类分类样本数据。
煤体瓦斯渗透率预测学习样本向量用X表示,样本集为:
[X]=[x1,x2,x3,x4]T=[σ,T,P,σm]T (5)
其中σ为有效应力(MPa)、T为温度(℃)、P为瓦斯压力(MPa)、σm为抗压强度(MPa)
输入层输入集为:
竞争层的样本用Y表示为:
Y=(y1,y2,…ym)T,yi∈{0,1},i=1,2,…,m (7)
输入层和竞争层之间的权值矩阵用W1:
竞争层到输出层之间的权值矩阵用W2:
网络的期望值输出用d表示为:
d=(d1,d2,…dl)T (10)
输出层的输出用O表示,输出集为:
[O]=[o1,o2]T=[1,0]T (11)
其中,o1为陡峭型类别,o2为平缓型类别。
CPSO-BP预测模型
在LVQ将样本数据分为两大类后,采用CPSO-BP预测,CPSO算法采用Logistic映射混沌***,迭代公式为:
zi+1=μzi(1-zi),i=0,1,…,μ∈(2,4] (12)
式中:μ为控制变量。
当μ=4,0≤z0≤1时,Logistic完全处于混沌状态,利用μ=4时的混沌运动特性改进标准粒子群算法,即
zi+1=4zi(1-zi),i=0,1,… (13)
粒子群初始化阶段利用公式(13)产生出优良的初始群体作为初始粒子群,以减轻随机初始化粒子群对算法优化性能的影响。
BP神经网络结构的层数为3层,设R、S1、S2分别为输入层、隐含层、输出层神经元个数,其中的输入层神经元为LVQ分类后的样本数据,隐含层为S函数,最佳隐含节点数为3个,输出层为煤体瓦斯渗透率。BP神经网络的权值和阈值与粒子群个***置向量P的维数N之间的映射关系如公式(14)所示:
N=S1×R+S1×S2+S1+S2 (14)
学习样本数据
经过影响因素筛选,本发明获得了主要的“四因素”,其部分样本数据来自汤国水,***,韩军,等.基于MABC-SVM的含瓦斯煤体渗透率预测模型[J].中国安全生产科学技术,2015,02:11-16.中的试验数据,见表1。
表1煤体瓦斯渗透率学习样本基本数据
组号 | 有效应力/MP<sub>a</sub> | 温度/℃ | 瓦斯压力/MP<sub>a</sub> | 抗压强度/MP<sub>a</sub> | 渗透率/(10<sup>-15</sup><sub>m</sub><sup>2</sup>) |
1 | 1.51 | 40 | 1.80 | 10.85 | 0.881 |
2 | 4.01 | 55 | 0.50 | 12.85 | 1.062 |
3 | 2.50 | 30 | 0.50 | 13.13 | 0.579 |
4 | 3.62 | 30 | 1.00 | 12.85 | 0.565 |
5 | 2.00 | 50 | 1.00 | 11.62 | 0.479 |
6 | 3.00 | 30 | 1.00 | 12.85 | 0.810 |
7 | 1.00 | 70 | 1.00 | 11.50 | 0.451 |
8 | 1.23 | 40 | 1.50 | 13.50 | 0.972 |
9 | 1.95 | 40 | 1.20 | 11.62 | 0.782 |
10 | 3.78 | 20 | 0.60 | 12.50 | 0.617 |
11 | 2.50 | 50 | 1.50 | 12.70 | 0.540 |
12 | 1.68 | 35 | 0.70 | 12.85 | 1.028 |
13 | 3.00 | 40 | 0.50 | 10.62 | 0.791 |
14 | 1.23 | 30 | 1.50 | 12.37 | 1.023 |
15 | 2.64 | 65 | 1.30 | 11.62 | 0.504 |
16 | 2.50 | 30 | 1.00 | 12.85 | 0.656 |
17 | 2.04 | 50 | 1.60 | 12.85 | 0.635 |
18 | 2.50 | 50 | 1.00 | 11.62 | 0.588 |
19 | 2.00 | 30 | 1.50 | 12.37 | 0.769 |
20 | 1.73 | 45 | 1.80 | 14.13 | 0.805 |
21 | 2.01 | 30 | 0.60 | 12.13 | 0.851 |
22 | 2.00 | 55 | 1.90 | 12.85 | 0.613 |
23 | 2.50 | 25 | 0.50 | 14.15 | 0.863 |
24 | 2.00 | 60 | 1.60 | 12.37 | 0.608 |
25 | 1.45 | 35 | 1.50 | 14.12 | 0.955 |
26 | 1.75 | 45 | 0.50 | 12.37 | 1.004 |
27 | 2.00 | 50 | 1.00 | 11.62 | 0.724 |
28 | 1.83 | 55 | 1.20 | 13.47 | 0.705 |
29 | 3.00 | 50 | 1.00 | 11.62 | 0.516 |
30 | 2.00 | 30 | 0.50 | 14.13 | 1.042 |
31 | 1.48 | 30 | 0.70 | 14.15 | 1.236 |
32 | 1.48 | 30 | 1.00 | 12.85 | 1.054 |
33 | 1.62 | 50 | 1.00 | 11.62 | 0.891 |
34 | 2.5 | 70 | 1.00 | 11.50 | 0.516 |
35 | 1.29 | 30 | 1.50 | 11.37 | 0.669 |
36 | 3.00 | 30 | 1.50 | 12.37 | 0.632 |
37 | 2.20 | 70 | 1.20 | 11.50 | 0.564 |
38 | 1.45 | 70 | 1.00 | 11.50 | 0.786 |
39 | 1.73 | 30 | 0.50 | 14.13 | 1.189 |
40 | 2.00 | 70 | 1.00 | 11.50 | 0.632 |
学习样本训练及预测结果
本发明采用Matlab2015神经网络工具箱结合表1数据编写了BP、GA-BP、PSO-BP、CPSO-BP四种组合方式渗透率预测模型代码,设定BP神经网络迭代次数1500次,学习率0.01,最小均方根误差10-5;CPSO优化神经网络的粒子取W=0.7298,c1=c2=1.4962,粒子的上界值为5,下界值为-5。
LVQ训练及分类结果
将表1样本数据输入LVQ进行训练,首先进行网络初始化:设定较小随机数作为初始权值W1和W2;确定初始学***均向量。通过期望值最小寻找获胜神经元。以胜出神经元小于本发明指定的有效应力平均值(约2.5MPa)的划分为正确分类,否则调整分类时的权值,更新学习效率,直到t=tm,获得分类结果。
以有效应力、温度、瓦斯压力、抗压强度的平均值作为输入样本,规定以有效应力平均值为拐点的临界值左边和右边的输出类型依次为1,0,其分类识别结果如表2所示。
表2识别正确率对比
从表2可以看出,采用LVQ模型对有效应力拐点划分具有较高的精度,分类识别正确率最高可达96.15%。
CPSO-BP训练及预测结果
将表1样本数据输入LVQ进行分类后的两组样本数据输入CPSO-BP进行训练,首先初始化BP神经网络3层结构,输入层神经元为4(有效应力、温度、瓦斯压力、抗压强度),隐含层神经元为5,输出层神经元为1(煤体瓦斯渗透率)。通过混沌初始化种群大小、维度D及粒子S速度和位置及均方误差粒子群适应度函数自寻优,将寻找出的最优解赋给BP网络的权值及阈值,继续BP网络的学习训练预测。
本发明采用四种预测模型进行了40组样本数据的学习训练,预测模型一为BP神经网络预测,预测模型二为GA-BP神经网络预测,预测模型三为PSO-BP神经网络预测,预测模型四为CPSO-BP。预测模型中的输入数据均为LVQ分类后的数据,其中类型1训练样本为26组;类型0训练样本为14组。四种模型两类预测曲线和汤国水,***,韩军,等.基于MABC-SVM的含瓦斯煤体渗透率预测模型[J].中国安全生产科学技术,2015,02:11-16.的试验测试曲线如图5所示。
图5的对比仿真结果充分表明:本发明建立的煤体瓦斯渗透率预测模型拟合度高,预测效果好。
类型1的煤体瓦斯渗透率高且变化量大,但预测精度高,能更有效避免安全隐患,因此按照有效应力拐点采用LVQ进行样本数据分类并结合CPSO-BP进行优化预测的方法,有利于煤矿开采过程中瓦斯的安全隐患预警。
评估分析
依据建立的误差评估方法对煤体瓦斯渗透率的四种预测结果与试验真实值之差进行误差计算,如表3所示。
表3四种预测模型误差评价指标结果
由表3可知,LVQ-CPSO-BP算法的煤体瓦斯渗透率预测精度为95.26%,优于其他三种模型,所建立的LVQ-CPSO-BP预测模型有效、适应性更强、预测精度更高。
本发明的工作原理为:
(1)提出了LVQ-CPSO-BP算法的煤体瓦斯渗透率预测方法,获得了有效应力小于拐点段煤体瓦斯渗透率预测精度更好的效果,有利于煤矿瓦斯安全隐患的预警。
(2)采用宏观和微观分析法,筛选出了影响煤体瓦斯渗透率的主要宏观因素为煤层埋深,主要微观因素为有效应力、温度、瓦斯压力、抗压强度。
(3)依据煤层埋深的分布特性,提出了临界点分层概念;依据煤体瓦斯渗透率与有效应力之间的关系存在拐点特性,提出了拐点分类概念。
(4)基于Matlab编程实现了LVQ-CPSO-BP算法的煤体瓦斯渗透率预测,选取“四要素”样本数据,获得了预测结果,验证了经LVQ分类后的CPSO-BP神经网络预测精度优于PSO-BP神经网络、GA-BP神经网络和BP神经网络。
(5)基于预测误差评价指标,计算了4种预测模型的误差值,其中基于LVQ分类后的CPSO-BP预测模型各种误差值分别为:MRE为6.36%,MAE为0.0852×10-15m2,RMSE为0.1741×10-15m2,FMA为95.26%,表明该模型具有较好地精确度和收敛性,预测结果更加逼近真实值。
(6)多组合算法预测精度明显优于单一算法,也优于双组合算法,验证了LVQ-CPSO-BP多组合算法在煤体瓦斯渗透率预测精度提高方面可行、有效。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何不经过创造性劳动想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求书所限定的保护范围为准。
Claims (10)
1.基于LVQ-CPSO-BP算法的煤体瓦斯渗透率预测方法,其特征在于,由如下步骤构成:
步骤一、煤层埋深划分及选择:对煤层埋深进行研究:煤层埋深分为浅埋层和非浅埋层,而浅埋层一般在150m以上,处于瓦斯风化带;而非浅埋层是指150m以下,最深达到1500m左右,大部分煤层处于200m-700m之间,将非浅埋层自定义分为两层,设临界值为700m,以浅为中埋层,以深为深埋层;
步骤二,主要影响因素选择:选择有效应力、温度、瓦斯压力、抗压强度作为主要影响因素;
有效应力与煤体瓦斯渗透率之间的关系出现拐点,有效应力小于拐点时,煤体瓦斯渗透率变化较快,大于拐点时,煤体瓦斯渗透率变化平缓,将2MPa设为有效应力拐点;
步骤三、在煤层埋深非浅埋层宏观划分基础上,采用LVQ算法对有效应力样本数据按拐点进行微观分类;
步骤四、建立基于CPSO-BP算法的煤体瓦斯渗透率预测子模型;
步骤五、通过LVQ-CPSO-BP预测***的训练,获得煤体瓦斯渗透率预测结果。
2.根据权利要求1所述的方法,其特征在于,为了评价所建立预测模型的可信度及精确度,对误差和精度进行量化,采用式(1)-式(4)预测误差和精度评价指标进行煤体瓦斯渗透率预测结果定量评估:
1)预测误差的可靠性用平均相对误差:
2)预测值偏离真值误差的大小用平均绝对误差:
3)样本数据离散程度用均方根误差:
4)预测模型精确度:
其中:Xi为真实值,X′i为与预测值,M为样本数量。
3.根据权利要求1所述的方法,其特征在于,所述步骤三中,采用LVQ以有效应力拐点为界将其分为煤体瓦斯渗透率变化平缓型和陡峭型两种类型,以样本数据中有效应力为1-5MPa的平均值温度平均值瓦斯压力平均值抗压强度平均值作为LVQ的输入,经过LVQ分类识别训练后,获得两类分类样本数据。
4.根据权利要求3所述的方法,其特征在于,所述步骤三中,煤体瓦斯渗透率预测学习样本向量用X表示,样本集为:
[X]=[x1,x2,x3,x4]T=[σ,T,P,σm]T (5)
其中σ为有效应力单位:MPa、T为温度单位:℃、P为瓦斯压力单位:MPa、σm为抗压强度单位:MPa;
输入层输入集为:
竞争层的样本用Y表示为:
Y=(y1,y2,…ym)T,yi∈{0,1},i=1,2,…,m (7)
输入层和竞争层之间的权值矩阵用W1:
竞争层到输出层之间的权值矩阵用W2:
网络的期望值输出用d表示为:
d=(d1,d2,…dl)T (10)
输出层的输出用O表示,输出集为:
[O]=[o1,o2]T=[1,0]T (11)
其中,o1为陡峭型类别,o2为平缓型类别。
5.根据权利要求4所述的方法,其特征在于,所述步骤三中、经过影响因素筛选,以有效应力,温度,瓦斯压力,抗压强度与渗透率的对应数据为样本进行学习;
采用Matlab2015神经网络工具箱结合样本数据编写了BP、GA-BP、PSO-BP、CPSO-BP四种组合方式渗透率预测模型代码,设定BP神经网络迭代次数1500次,学习率0.01,最小均方根误差10-5;CPSO优化神经网络的粒子取W=0.7298,c1=c2=1.4962,粒子的上界值为5,下界值为-5。
6.根据权利要求5所述的方法,其特征在于,所述步骤四中,在LVQ将样本数据分为两大类后,采用CPSO-BP预测,CPSO算法采用Logistic映射混沌***,迭代公式为:
zi+1=μzi(1-zi),i=0,1,…,μ∈(2,4] (12)
式中:μ为控制变量;
当μ=4,0≤z0≤1时,Logistic完全处于混沌状态,利用μ=4时的混沌运动特性改进标准粒子群算法,即
zi+1=4zi(1-zi),i=0,1,… (13)
粒子群初始化阶段利用公式(13)产生出优良的初始群体作为初始粒子群,以减轻随机初始化粒子群对算法优化性能的影响。
7.根据权利要求6所述的方法,其特征在于,所述步骤四中,将BP神经网络结构的层数设为3层,设R、S1、S2分别为输入层、隐含层、输出层神经元个数,其中的输入层神经元为LVQ分类后的样本数据,隐含层为S函数,最佳隐含节点数为3个,输出层为煤体瓦斯渗透率,BP神经网络的权值和阈值与粒子群个***置向量P的维数N之间的映射关系如公式(14)所示:
N=S1×R+S1×S2+S1+S2 (14)。
8.根据权利要求7所述的方法,其特征在于,所述步骤五中、将样本数据输入LVQ进行训练,首先进行网络初始化:设定较小随机数作为初始权值W1和W2;确定初始学***均向量,通过期望值最小寻找获胜神经元,以胜出神经元小于指定的有效应力平均值的划分为正确分类,否则调整分类时的权值,更新学习效率,直到t=tm,获得分类结果;
表1
9.根据权利要求8所述的方法,其特征在于,所述步骤五中,以有效应力、温度、瓦斯压力、抗压强度的平均值作为输入样本,规定以有效应力平均值为拐点的临界值左边和右边的输出类型依次为1,0,得到分类识别结果;
将样本数据输入LVQ进行分类后的两组样本数据输入CPSO-BP进行训练,首先初始化BP神经网络3层结构,输入层神经元为4:有效应力、温度、瓦斯压力、抗压强度;隐含层神经元为5,输出层神经元为1:煤体瓦斯渗透率;通过混沌初始化种群大小、维度D及粒子S速度和位置及均方误差粒子群适应度函数自寻优,将寻找出的最优解赋给BP网络的权值及阈值,继续BP网络的学习训练预测。
10.根据权利要求9所述的方法,其特征在于,采用四种预测模型进行了40组样本数据的学习训练,预测模型一为BP神经网络预测,预测模型二为GA-BP神经网络预测,预测模型三为PSO-BP神经网络预测,预测模型四为CPSO-BP;预测模型中的输入数据均为LVQ分类后的数据,其中类型1训练样本为26组;类型0训练样本为14组,采用LVQ模型对有效应力拐点划分具有较高的精度,分类识别正确率最高可达96.15%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710121485.3A CN106869990B (zh) | 2017-03-02 | 2017-03-02 | 基于lvq-cpso-bp算法的煤体瓦斯渗透率预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710121485.3A CN106869990B (zh) | 2017-03-02 | 2017-03-02 | 基于lvq-cpso-bp算法的煤体瓦斯渗透率预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106869990A CN106869990A (zh) | 2017-06-20 |
CN106869990B true CN106869990B (zh) | 2019-03-19 |
Family
ID=59169287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710121485.3A Expired - Fee Related CN106869990B (zh) | 2017-03-02 | 2017-03-02 | 基于lvq-cpso-bp算法的煤体瓦斯渗透率预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106869990B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109191423B (zh) * | 2018-07-18 | 2019-07-02 | 中国矿业大学 | 一种基于机器图像智能学习的多孔介质渗透率预测方法 |
CN109120961B (zh) * | 2018-07-20 | 2020-11-03 | 南京邮电大学 | 基于PNN-PSO算法的IPTV不平衡数据集的QoE的预测方法 |
CN109657361B (zh) * | 2018-12-21 | 2022-09-13 | 成都四方伟业软件股份有限公司 | 瓦斯含量预测方法、终端设备及计算机可读存储介质 |
CN110633504A (zh) * | 2019-08-21 | 2019-12-31 | 中联煤层气有限责任公司 | 一种煤层瓦斯渗透率预测方法 |
CN111461457A (zh) * | 2020-04-17 | 2020-07-28 | 福州大学 | 一种基于粒子群优化bp神经网络的基坑位移预测方法 |
CN111521537B (zh) * | 2020-04-30 | 2022-10-21 | 太原理工大学 | 一种块煤钻进过程的多维数据测量装置 |
CN111398130B (zh) * | 2020-04-30 | 2022-11-08 | 太原理工大学 | 多维度数据来源的块煤渗透率分析方法、测量装置及方法 |
CN111929285A (zh) * | 2020-08-11 | 2020-11-13 | 安徽理工大学 | 一种激光诱导荧光红酒掺杂酒精的光谱识别方法 |
CN113570165B (zh) * | 2021-09-03 | 2024-03-15 | 中国矿业大学 | 基于粒子群算法优化的煤储层渗透率智能预测方法 |
CN116930042B (zh) * | 2023-09-19 | 2023-12-01 | 常州建昊建筑鉴定检测有限公司 | 建筑防水材料性能检测设备及方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316448A (zh) * | 2014-11-03 | 2015-01-28 | 中国石油大学(华东) | 一种高阶煤岩气相渗透率动态变化的预测方法 |
CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
CN104612635A (zh) * | 2014-12-10 | 2015-05-13 | 中煤科工集团重庆研究院有限公司 | 煤层群瓦斯联合抽采达标预评判方法 |
CN104794550A (zh) * | 2015-05-13 | 2015-07-22 | 山东科技大学 | 基于wt-kpca-svr耦合模型的瓦斯涌出量预测方法 |
CN106202748A (zh) * | 2016-07-14 | 2016-12-07 | 中国海洋石油总公司 | 一种基于渗透率和井距的煤层气采收率预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6860147B2 (en) * | 2002-09-30 | 2005-03-01 | Alberta Research Council Inc. | Process for predicting porosity and permeability of a coal bed |
-
2017
- 2017-03-02 CN CN201710121485.3A patent/CN106869990B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104500050A (zh) * | 2014-10-31 | 2015-04-08 | 中国石油大学(华东) | 一种裂缝性储层渗透率张量及各向异性定量预测方法 |
CN104316448A (zh) * | 2014-11-03 | 2015-01-28 | 中国石油大学(华东) | 一种高阶煤岩气相渗透率动态变化的预测方法 |
CN104612635A (zh) * | 2014-12-10 | 2015-05-13 | 中煤科工集团重庆研究院有限公司 | 煤层群瓦斯联合抽采达标预评判方法 |
CN104794550A (zh) * | 2015-05-13 | 2015-07-22 | 山东科技大学 | 基于wt-kpca-svr耦合模型的瓦斯涌出量预测方法 |
CN106202748A (zh) * | 2016-07-14 | 2016-12-07 | 中国海洋石油总公司 | 一种基于渗透率和井距的煤层气采收率预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106869990A (zh) | 2017-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106869990B (zh) | 基于lvq-cpso-bp算法的煤体瓦斯渗透率预测方法 | |
CN110363344A (zh) | 基于miv-gp算法优化bp神经网络的概率积分参数预测方法 | |
CN110348713A (zh) | 一种基于关联性分析和数据挖掘的台区线损计算方法 | |
Wang et al. | Marcellus shale lithofacies prediction by multiclass neural network classification in the Appalachian Basin | |
Zhao et al. | A classification-based surrogate-assisted multiobjective evolutionary algorithm for production optimization under geological uncertainty | |
Sahoo et al. | Pattern recognition in lithology classification: modeling using neural networks, self-organizing maps and genetic algorithms | |
Aïfa et al. | Neuro-fuzzy system to predict permeability and porosity from well log data: A case study of Hassi R׳ Mel gas field, Algeria | |
Amiri et al. | Improving water saturation estimation in a tight shaly sandstone reservoir using artificial neural network optimized by imperialist competitive algorithm–A case study | |
CN105913450A (zh) | 基于神经网络图像处理的轮胎橡胶炭黑分散度评价方法及*** | |
Brown et al. | Use of fuzzy membership input layers to combine subjective geological knowledge and empirical data in a neural network method for mineral-potential mapping | |
Karkevandi-Talkhooncheh et al. | Application of hybrid adaptive neuro-fuzzy inference system in well placement optimization | |
CN108053052B (zh) | 一种油罐车油气泄漏速度智能监测*** | |
CN106295199A (zh) | 基于自动编码器和多目标优化的自动历史拟合方法及*** | |
CN114139458B (zh) | 基于机器学习的钻井参数优化方法 | |
Wang | Fuzzy clustering analysis by using genetic algorithm | |
CN101980298A (zh) | 基于多智能体遗传聚类算法的图像分割方法 | |
CN104680025A (zh) | 基于遗传算法极限学习机的抽油机参数优化方法 | |
Wang et al. | A critical review on intelligent optimization algorithms and surrogate models for conventional and unconventional reservoir production optimization | |
Shokir et al. | Permeability estimation from well log responses | |
Verma et al. | Quantification of sand fraction from seismic attributes using Neuro-Fuzzy approach | |
Olatunji et al. | Modeling permeability prediction using extreme learning machines | |
Chatterjee et al. | Ore grade estimation of a limestone deposit in India using an artificial neural network | |
CN113033637A (zh) | 一种基于岩石物理样本增强的集成学习岩相预测方法 | |
Sfidari et al. | Electrofacies clustering and a hybrid intelligent based method for porosity and permeability prediction in the South Pars Gas Field, Persian Gulf | |
CN115660221B (zh) | 基于混合神经网络的油气藏经济可采储量评估方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190319 |
|
CF01 | Termination of patent right due to non-payment of annual fee |