专利名称：一种基于分量热力学基因表达式编程的矿岩强度预测方法
技术领域：
本发明涉及一种矿岩強度预测方法，尤其是涉及ー种基于分量热力学基因表达式编程的矿岩强度预测方法。
背景技术：
矿岩强度问题是矿业工程中的基本问题，但是由于传统的矿岩强度预测方法需要进行大量的现场试验，从而导致人力、资金等急剧增加。在实际工程实践中受资金、设备等条件的制约，往往无法进行大量的现场试验。因此如何高效、准确地预测矿岩強度，是许多工程人员十分关注的重要问题。目前，用于矿岩強度的预测方法有很多，主要可以分为数学物理预测方法和智能预测方法。数学物理预测方法需要工程人员掌握比较多的专业知识，因此普适性不是很強。而智能预测方法则是融合计算机科学技木，尤其是智能计算技术而发展的一种简单易行，无需工程人员掌握很多的领域知识，并且具有灵活、通用性强，可 满足复杂多样的实际工程需求等优点。因此智能预测方法具有广阔的前景，是当前矿岩强度预测研究的热点。目前，矿岩强度的智能预测方法主要是神经网络，支持向量机等方法。但这些方法存在着很多不足，实际工程效果并不是很理想。神经网络和支持向量机的方法存在的主要问题是需要采集矿岩试件数量较多，当矿岩试件数量较少容易过度拟合，所获得的模型精度很有限，往往无法达到实际工程要求。而当矿岩试件数量较多时，一方面会导致人力、资金等増加；另一方面也増加了算法学习过程的复杂性，从而导致学习时间长。可以參见有关文献基于粒子群支持向量机的矿岩强度指标的超声预测.近年来，智能计算技术的迅速发展，出现了基因表达式编程算法。基因表达式编程算法是模拟达尔文生物进化论的自然选择和遗传学机理的生物进化过程的一种有效的智能计算方法，它采用固定长度的线性字符串表示解，因此它具有简洁、高性能和强鲁棒性等优点，已经成功应用到众多的实际工程领域中。然而在实际工程实践中发现，传统基因表达式编程算法应用于矿岩强度预测时，存在着收敛速度慢和容易陷入局部最优的问题。这两个问题的原因在很大程度上是由于在传统基因表达式编程算法的演化操作过程中，无法定量地协调选择压カ和种群多祥性之间的平衡，从而容易导致选择压カ和种群多祥性出现随意抖动，影响了传统基因表达式编程算法的收敛速度和求解精度。当传统基因表达式编程算法应用于矿岩强度预测时，选择压カ和种群多样性两者之间是相互冲突的。在通常情况下，当选择压力过大时，种群中接近于当前最优个体的个体就越多，种群的平均适应值就越优秀，可以使基因表达式编程算法收敛速度加快，但会导致种群中大部分的个体都趋向于当前最优个体的附近，种群多样性变差，増加算法陷入局部最优的概率；当选择压力过小时，虽然可以使得种群中个体分布趋于分散，种群多样性变好，増加算法收敛到全局最优解的概率，但这样会减慢基因表达式编程算法的收敛速度。目前，如何定量地协调基因表达式编程算法的选择压力和种群多祥性之间平衡的研究成果还很缺乏。

发明内容
本发明主要是解决现有技术所存在的技术问题；提供了ー种定量地协调基因表达式编程中选择压力和种群多祥性之间的平衡，以提高传统基因表达式编程的收敛速度、求解精度和算法的稳定性。利用提出的分量热力学基因表达式编程算法以矿岩试件的吸水率、干密度、波阻抗、动泊松比、动弹性模量为输入变量，抗压强度为输出变量，自动智能地挖掘出矿岩强度的数学模型，从而预测矿岩的抗压强度的一种基于分量热力学基因表达式编程的矿岩强度预测方法。本发明的上述技术问题主要是通过下述技术方案得以解决的一种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，包括以下步骤步骤1，采集N个矿岩样本，并针对N个矿岩样本的吸水率、干密度、波阻抗、动泊松比、动弹性模量以及抗压强度进行实验后获得每ー个矿岩样本的试验数据；并将N个矿 岩样本的试验数据记为矩阵A ；步骤2，用户自定义初始化參数，所述初始化參数包括种群大小PS，子种群大小M，最大评价次数MAX_FE，比例因子α，等级数K，Markov链长LK，初始温度T0，函数符和终结符，基因长度，基因个数，变异概率，插串概率，插串长度以及重组概率；步骤3，令当前演化代数t=0;温度下降因子k=0;温度T=TO;步骤4，产生初始化种群Pt，对每个个体的染色体解码成数学表达式，并利用矩阵A评价每个个体的适应值，其中吸水率、干密度、波阻抗、动泊松比、动弹性模量为输入变量，抗压强度为输出变量；然后保存最优个体，并计算第O代的活跃窗ロ Wtl ；步骤5，通过执行基因表达式编程的选择、变异、插串、重组算子对种群Pt中的PS个个体生成M个新个体，并将M个新个体组织成子种群0t，对Ot中的每个个体进行适应值评价；然后保存适应值最大的个体为最优个体；步骤6，将种群Pt中的PS个个体和子种群Ot中的M个个体组织成临时种群P’ t，然后计算第t+1代的活跃窗ロ Wt+1，并计算临时种群P’t的中的每个个体的自由能分量，再找出自由能分量值最大的前M个个体，然后删除这M个个体，得到由PS个个体组成的新一代种群Pt+1 ;步骤7，重复步骤5至步骤6直至评价次数达到MAX_FE后结束，将执行过程中得到的最优个体的染色体解码成数学表达式就可以得到以吸水率、干密度、波阻抗、动泊松比、动弹性模量为输入变量，抗压强度为输出变量的矿岩强度数学模型，利用获得的矿岩强度数学模型就可以预测矿岩的抗压强度。在上述的一种基于分量热力学基因表达式编程的矿岩强度预测方法，所述步骤4的具体操作步骤如下对每个个体的染色体解码成数学表达式，并利用矩阵A评价每个体的适应值，并保存适应值最大的个体为最优个体，设任意个体j的染色体解码成数学表达式记为&，其
适应值计算公式为i .Oc -10 + X (f} (An, An, Ai3, A14, Jis) - Ai6)其中，N 为矿岩样本总数，Ail, Ai2, Ai3, Ai4, Ai5, Ai6,分别为矩阵A中存储的第i矿岩样本的吸水率、干密度、波阻抗、动泊松比、动弹性模量和抗压强度，最优个体的判断标准是在所有个体中适应值最大的个体。在上述的一种基于分量热力学基因表达式编程的矿岩强度预测方法，所述步骤5中的具体操作步骤如下步骤5. 1，判断当前评价次数是否大于等最大评价次数MAX_FE，若是则转到步骤7，否则执行以下步骤；步骤5. 2，i=0 ;其中，i为临时变量，为记数器；步骤5. 3，判断i是否大于等于Markov链长LK，若是则执行温度下降因子k=k+l ；当前温度T=T/(l+k)，后执行步骤5.1，否则执行以下步骤；步骤5. 4，通过执行基因表达式编程的选择、变异、插串、重组等算子生成M个新个 体；步骤5. 5，将M个新个体组织成子种群0t，然后按公式
F(X) ニ---ΓI .Oe —10 + Σ (/ノ (A11, A12，A13，Al4，Al5) — Α,6)计算Ot中的每个个体的适应值，并且保存适应值最大的个体为最优个体。在上述的一种基于分量热力学基因表达式编程的矿岩强度预测方法，所述步骤6的具体操作步骤如下步骤6. 1，计算第t+Ι代的活跃窗ロ fft+1 ;并令临时种群P，t=Pt U Ot ;，其中Pt是当前种群，Ot是子种群步骤6. 2，计算P’ t中的PS+M个个体的相对能量；步骤6. 3，根据步骤6. 2中计算的相对能量，计算P’t中的PS+M个个体的自由能
分量;步骤6. 4，从P’ t中找出自由能分量最大的前M个个体；步骤6. 5，将步骤6. 4中找出的M个具有最大自由能分量的个体从P’ t中删除；步骤6. 6,当前演化代数t=t+l ；,下一代新种群Pt=P’ t ；i=i+l,后执行步骤5. 3。在上述的一种基于分量热力学基因表达式编程的矿岩强度预测方法，所述步骤6. 2中，计算每个个体的相对能量方法如下e'(Wt X ) = , ' )，其中F(X)为个体X的适应值函数，且其特性为个体越优，则适应值越大。在上述的一种基于分量热力学基因表达式编程的矿岩强度预测方法，所述步骤6. I中，计算第t代的活跃窗ロ Wt的方法如下当t=0 时，W。= [I0, U。],其中I0 = min{_F(X) X e P。}，u0 = max{_F(X) X e P0}；当 t>0 时,若 Wt^1 = [It^1, ut_J ,且第 t 代子种群0t (XN+1, XN+2,...，XN+M) e S，则 Wt = [lt, Ut]，其中It = min (I^1, min {_F (X) | X e Oj), Ut = max (U^1, max {_F (X) | X e Oj)。在上述的一种基于分量热力学基因表达式编程的矿岩强度预测方法，所述步骤6. 3中，计算个体的自由能分量方法如下
权利要求
1.一种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，包括以下步骤 步骤1，采集N个矿岩样本，并针对N个矿岩样本的吸水率、干密度、波阻抗、动泊松比、动弹性模量以及抗压强度进行实验后获得每ー个矿岩样本的试验数据；并将N个矿岩样本的试验数据记为矩阵A ； 步骤2，用户自定义初始化參数，所述初始化參数包括种群大小PS，子种群大小M，最大评价次数MAX_FE，比例因子α，等级数K，Markov链长LK，初始温度T0，函数符和终结符，基因长度，基因个数，变异概率，插串概率，插串长度以及重组概率； 步骤3，令当前演化代数t=0;温度下降因子k=0;温度T=TO; 步骤4，产生初始化种群Pt，对每个个体的染色体解码成数学表达式，并利用矩阵A评价每个个体的适应值，其中吸水率、干密度、波阻抗、动泊松比、动弹性模量为输入变量，抗压强度为输出变量；然后保存最优个体，并计算第O代的活跃窗ロ Wtl ； 步骤5，通过执行基因表达式编程的选择、变异、插串、重组算子对种群Pt中的PS个个体生成M个新个体，并将M个新个体组织成子种群0t，对Ot中的每个个体进行适应值评价；然后保存适应值最大的个体为最优个体； 步骤6，将种群Pt中的PS个个体和子种群Ot中的M个个体组织成临时种群P’t，然后计算第t+Ι代的活跃窗ロ Wt+1，并计算临时种群P’ t的中的每个个体的自由能分量，再找出自由能分量值最大的前M个个体，然后删除这M个个体，得到由PS个个体组成的新一代种群 Pt+i ； 步骤7，重复步骤5至步骤6直至评价次数达到MAX_FE后结束，将执行过程中得到的最优个体的染色体解码成数学表达式就可以得到以吸水率、干密度、波阻抗、动泊松比、动弹性模量为输入变量，抗压强度为输出变量的矿岩强度数学模型，利用获得的矿岩强度数学模型就可以预测矿岩的抗压强度。
2.根据权利要求I所述的ー种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，所述步骤4的具体操作步骤如下 对每个个体的染色体解码成数学表达式，并利用矩阵A评价每个体的适应值，并保存适应值最大的个体为最优个体，设任意个体j的染色体解码成数学表达式记为ち，其适应值计算公式为LOe —10 + ^(/；(為"為“ふも)—べ6)其中，N为矿岩样本总数，Ail, Ai2, Ai3, Ai4, Ai5, Ai6,分别为矩阵A中存储的第i矿岩样本的吸水率、干密度、波阻抗、动泊松比、动弹性模量和抗压强度，最优个体的判断标准是在所有个体中适应值最大的个体。
3.根据权利要求I所述的ー种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，所述步骤5中的具体操作步骤如下 步骤5. 1，判断当前评价次数是否大于等最大评价次数MAX_FE，若是则转到步骤7，否则执行以下步骤； 步骤5. 2i=0 ;其中，i为临时变量，为记数器； 步骤5. 3，判断i是否大于等于Markov链长LK，若是则执行温度下降因子k=k+l ;当前温度T=T/(l+k)，后执行步骤5.1，否则执行以下步骤； 步骤5. 4，通过执行基因表达式编程的选择、变异、插串、重组等算子生成M个新个体； 步骤5. 5，将M个新个体组织成子种群Ot，然后按公式 计算Ot中的每个个体的适应值，并且保存适应值最大的个体为最优个体。
4.根据权利要求I所述的ー种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，所述步骤6的具体操作步骤如下 步骤6. 1，计算第t+Ι代的活跃窗ロ Wt+1 ;并令临时种群P’ t=Pt U Ot ；，其中Pt是当前种群，Ot是子种群 步骤6. 2，计算P’ t中的PS+M个个体的相对能量； 步骤6. 3，根据步骤6. 2中计算的相对能量，计算P’t中的PS+M个个体的自由能分量； 步骤6. 4，从P’ t中找出自由能分量最大的前M个个体； 步骤6. 5，将步骤6. 4中找出的M个具有最大自由能分量的个体从P’t中删除； 步骤6. 6,当前演化代数t=t+l ；,下一代新种群Pt=P’ t ；i=i+l,后执行步骤5. 3。
5.根据权利要求4所述的ー种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，所述步骤6. 2中，计算每个个体的相对能量方法如下 其中F(X)为个体X的适应值函数，且其特性为个体越优，则适应值越大。
6.根据权利要求4所述的ー种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，所述步骤6. I中，计算第t代的活跃窗ロ Wt的方法如下 当 t=0 时，W。= [I0, U。],其中I0 = min{-F(X) |X e P0}，u0 = max{-F(X) |X e PcJ ;当 t>0 时，若 Ww = [It^1, ut_J ,且第 t 代子种群0t (XN+1, XN+2,…，XN+M) e S，则 Wt = [lt, ut]，其中 It = min (Ih, min {-F (X) | X e Oj), ut = max Ov1, max {_F (X) | X e Oj)。
7.根据权利要求4所述的ー种基于分量热力学基因表达式编程的矿岩强度预测方法，其特征在于，所述步骤6.3中，计算个体的自由能分量方法如下 其中，T为当前温度，个体X在种群Pt的等级为b，种群Pt中落入Wt的第b个等级尤中的个体数目为nb，且 中，a为比例因子，K为等级数。
全文摘要
本发明涉及一种基于分量热力学基因表达式编程的矿岩强度预测方法。本发明利用提出的分量热力学基因表达式编程算法以矿岩试件的吸水率、干密度、波阻抗、动泊松比、动弹性模量为输入变量，抗压强度为输出变量，自动智能地挖掘出矿岩强度的数学模型，从而预测矿岩的抗压强度。本发明能定量地协调基因表达式编程中选择压力和种群多样性之间的平衡，从而提高了传统基因表达式编程应用于矿岩强度预测的收敛速度、求解精度和算法的稳定性。
文档编号G01N33/24GK102854299SQ201210261469
公开日2013年1月2日 申请日期2012年7月26日 优先权日2012年7月26日
发明者郭肇禄, 吴志健, 董晓健, 李元香, 张勇 申请人:武汉大学