CN102692491B - 基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 - Google Patents
基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 Download PDFInfo
- Publication number
- CN102692491B CN102692491B CN201210208585.7A CN201210208585A CN102692491B CN 102692491 B CN102692491 B CN 102692491B CN 201210208585 A CN201210208585 A CN 201210208585A CN 102692491 B CN102692491 B CN 102692491B
- Authority
- CN
- China
- Prior art keywords
- optimizing
- soil moisture
- soil
- stage
- algorithm
- 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
Landscapes
- General Factory Administration (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供一种基于分阶段禁忌寻优算法的土壤水分特征参数计算方法,包括:S1:获取土壤饱和导水率和饱和含水量;S2:根据S1:计算土壤非饱和导水率和土壤水分扩散率;S3:建立一维垂直土壤水分运动数学模型;S4:得到土壤含水率理论值计算公式θ(a,b,m);S5:确定土壤含水率实际测量值θ;S6:根据S4和S5中的θ(a,b,m)和θ,并根据改进的分段禁忌算法进行优化,确定出a、b、m的最优值;S7:将S6中的最优值带入S3中计算出土壤水分特征参数。本发明公开的方法能够确定在负压地下灌溉条件下,一维土壤水分的运动规律,从而指导并进一步改进负压地下灌溉系统,提高土壤水分利用率,进而提高土壤水分利用率。
Description
技术领域
本发明涉及农业技术领域,特别涉及一种土壤水分特征参数的计算方法,尤其涉及一种基于分阶段禁忌寻优算法的土壤水分特征参数计算方法。
背景技术
土壤水分运动是一个复杂的过程,与灌溉方式、土壤条件和外界环境等众多因素密切相关。在特定的灌溉条件下,研究土壤水分运动,对理解土壤水分分布,进而改进灌溉系统,提高水分利用效率具有重要意义。
土壤水分运动模拟研究中,通常是在特定的灌溉条件、土壤组成和外界环境条件下,利用实验方法建立经验公式来模拟土壤水分运动规律,或者利用已有的成熟公式,利用参数优化算法,结合实验数据,计算出优化参数,进而确定出土壤水分运动规律。
目前应用于土壤水分运动参数优化的算法主要有遗传算法,粒子群算法,最小二乘法等。
遗传算法是目前应用最广泛的研究土壤水分运动参数的算法。遗传算法从串集开始搜索,覆盖面大,利于全局择优;采用概率的变迁规则来指导搜索方向;具有自组织、自适应和自学习性。但是遗传算法较容易陷入局部最优点,并且迭代次数较多。
粒子群算法是近年来发展起来的一种新的进化算法。从随机解出发,通过迭代寻找最优解,也是通过适应度来评价解的品质,但比遗传算法规则更为简单,算法在多维空间函数寻优、动态目标寻优等方面有着收敛速度快、解质量高、鲁棒性好、需调整的参数较少等优点。但是也存在着精度较低,易发散等缺点。
最小二乘法通过最小化误差的平方和寻找参数的最佳匹配函数。目前应用最小二乘法主要是针对Van Genuchten方程描述的土壤水分特征曲线进行参数寻优。利用PS最小二乘法求解Van Genuchten方程参数算法实现简单、运算速度快、计算精度高。但是对其他土壤水分运动方程应用的较少。
禁忌搜索(Tabu Search)算法是一对局部邻域搜索的扩展,为一种全局迭代寻优算法。通过模拟人的经验,利用一种灵活的存储结构(禁忌表)和相应的禁忌准则来避免反复搜索,达到跳出局部最优解的目的。并通过特赦规则来赦免一些被禁忌的优良解,从而保证多样化的有效搜索,从而实现全局优化。禁忌算法具有很强的通用性,并且对问题无特殊要求,因而其应用领域非常广泛。尽管禁忌算法通过建立禁忌表来避免陷入局部最优,但是基本禁忌算法的收敛速度和普适性较差。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是,针对现有技术的不足,提出一种土壤水分特征参数的计算方法,对一维土壤水分垂直渗入模型参数进行寻优。通过数值模拟和实测实验交互验证,该算法计算出来的寻优结果很好的模拟了土壤水分垂直运动规律。
(二)技术方案
本发明提供一种基于分阶段禁忌寻优算法的土壤水分特征参数的计算方法,包括:S1:通过实验测量获取土壤饱和导水率和饱和含水量;S2:根据S1中的土壤饱和导水率和饱和含水量,计算出土壤非饱和导水率和土壤水分扩散率;S3:建立一维垂直土壤水分运动数学模型;S4:将S2中的土壤非饱和导水率和土壤水分扩散率带入到S3的模型中,转化为土壤含水率理论值计算公式θ(a,b,m);S5:根据入渗实验,确定土壤含水率的实际测量值θ;S6:建立最优化目标函数,并将S4和S5中的θ(a,b,m)和θ带入目标函数中,根据改进的分段禁忌算法进行优化,以确定出a、b、m的最优值;所述改进的分段禁忌算法为以基本禁忌算法为基础,分别在三个阶段中的每个阶段的参数进行对应的设置和选取;S7:将S6中的a、b、m最优值带入S3中的一维垂直土壤水分运动数学模型中计算出土壤水分特征参数。
更好地,在S6中利用改进的分段禁忌算法对目标函数进行参数寻优,所述改进的分段禁忌算法具体流程步骤包括:S6.1:对目标函数进行初步寻优;S6.2:对目标函数以S6.1寻优结果为基础进行全局寻优;S6.3:对目标函数以S6.2的寻优结果为基础进行快速搜索。
更好地,所述S6.1具体包括:S6.1.1:初始值设置,初始值以随机方式在整个取值空间内取得;禁忌表为空,并采用单向列表方式加以管理;S6.1.2:邻域点的选取,利用目标函数的差分来指导邻域点的选取。
更好地,所述S6.2具体包括:S6.2.1:初始值设置,本步骤的初始值以S6.1中的初步寻优输出结果作为起点;S6.2.2:邻域点设置,本步骤的邻域点设置依然利用目标函数的差分来指导邻域点的选取。
更好地,所述S6.3具体包括:S6.3.1:初始值设置:本阶段的初始值以S6.2中的全局寻优输出结果作为起点。禁忌表初值则为S6.2中的全局寻优得出的最优解S6.3.2:邻域点设置,在上述S6.1初步寻优阶段和S6.2全局寻优阶段的寻优基础之上,寻优的范围缩小为S6.2全局寻优阶段中最小取值范围。
更好地,所述S6.1.2具体包括:S6.1.2.1计算目标函数对每一个变量的差分及其在当前解的差分值;S6.1.2.2在整个取值范围内,利用多维正态分布函数来选取邻域点。
(三)有益效果
本发明公开的方法能够确定在负压地下灌溉条件下,一维土壤水分的运动规律,从而指导并进一步改进负压地下灌溉系统,提高土壤水利用率,进而提高土壤水分利用率。利用分段禁忌算法对土壤水分多参数寻优方法,能够有效确定土壤水分运移规律,对其他灌溉方式下,土壤水分运动规律的模拟具有重要指导意义。
附图说明
图1为本发明土壤水分特征参数的计算方法流程图;
图2为本发明改进的分段禁忌算法初步寻优阶段的方法流程图;
图3为本发明改进的分段禁忌算法全局寻优阶段的方法流程图;
图4为本发明改进的分段禁忌算法加速寻优阶段的方法流程图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
优化算法一般用于解决最优化问题。所谓最优化问题是指设计一个目标函数,通过选取参数,从而使得目标函数达到最小值。
优化算法一般分为经典优化算法和启发式优化算法。其中:经典优化算法是以一个可行解作为初始值,按照确定的搜索策略,来寻找最优值。经典算法是利用局部信息(如单个初始点及其导数等)进行寻优,使得经典算法无法避免局部极小问题。启发式优化算法是指受大自然的运行规律或者人类积累的工作经验启发而来的方法。启发式优化算法均是从随机可行初始解出发,按照一定搜索策略来不断逼近问题的最优解。由于启发式优化算法对目标函数性质无特别要求,因此可以很好解决实际应用问题。但是启发式优化算法的计算量一般很大,且算法的普适性较差。
禁忌算法属于启发式优化算法。通过对局部邻域搜索的扩展,为一种全局迭代寻优算法。通过模拟人的经验,利用一种灵活的存储结构(禁忌表)和相应的禁忌准则来避免反复搜索,达到跳出局部最优解的目的。并通过特赦规则来赦免一些被禁忌的优良解,从而保证多样化的有效搜索,从而实现全局优化。
在禁忌算法中,邻域、禁忌表、禁忌长度、特赦规则和终止规则为算法的关键所在,也是影响算法的求解精度、速度和普适性的重要因素。
本发明提出利用一种改进的分阶段禁忌多变量寻优算法,来求解一维垂直入渗土壤水分运动参数的方法。
如图1所示,具体步骤如下:
S1:通过实验测量获取土壤饱和导水率Ks和饱和含水量θs;
S2:根据S1中的土壤饱和导水率Ks和饱和含水量θs,计算出土壤非饱和导水率K(θ)和土壤水分扩散率D(θ);计算公式如下:
其中:θ为当前土壤含水量,θs为饱和土壤含水量,Ks为饱和导水率通过实验测得。
S3:建立一维垂直土壤水分运动数学模型
θ=θ1 t=0 h≥0
θ=θ2 t>0 h=0
θ=θ1 t>0 h→∞
其中,θ1为土壤中均匀分布的含水率,θ2为地下因接近储水布而维持不变的含水率,h:垂直距离,t:时间(min)。
S4:将S2中的土壤非饱和导水率和土壤水分扩散率带入到S3的模型(3)中,转化为土壤含水率理论值计算公式θ(a,b,m);
S5:根据入渗实验,确定土壤含水率的实际测量值θ;
S6:建立最优化目标函数(4),并将S4和S5中的θ(a,b,m)和θ带入(4)中;并根据改进的分段禁忌算法进行优化,以确定出a、b、m的最优值;所述改进的分段禁忌算法为以基本禁忌算法为基础,分别在三个阶段中的每个阶段的参数进行对应的设置和选取。
al≤a≤au
b1≤b≤bu
待求经验参数约束条件为ml≤m≤mu,其中i为采样点;。
S7:将S6中的a、b、m最优值带入S3中的一维垂直土壤水分运动数学模型中计算出土壤水分特征参数。
其中,在S6中利用改进的分段禁忌算法对目标函数进行参数寻优,所述改进的分段禁忌算法,以基本禁忌算法为基础,分别在三个阶段加以实现。对于每阶段来说,本质上均为一个基本的禁忌算法,但是各个阶段在初始值设置、邻域点选取、取值范围设置等方面均有所不同,从而有效提高算法的全局最优解的精度和算法的收敛速度。
具体流程步骤包括:
如图2所示,初步寻优阶段:对目标函数进行初步寻优包括以下步骤:
S6.1.1:初始值设置。对于全局优化算法来说,初始值都是无法预测的,因此在初步寻优阶段,基本禁忌算法的初始值以随机方式在整个取值空间内取得。禁忌表为空,并采用单向列表方式加以管理;
S6.1.2:邻域点的选取。为了保证邻域点选取的更为有针对性,本发明利用目标函数的差分来指导邻域点的选取。具体的做法:
S6.1.2.1计算目标函数对每一个变量的差分及其在当前解的差分值;
S6.1.2.2在整个取值范围内,利用多维正态分布函数来选取邻域点。其中该多维正态分布函数的标准差是根据差分值来设置的。差分值越大,则方差越小。利用数学变换和正态分布的原则来将标准差设置在合理的范围。本邻域选取方法有效地避免了陷入局部最优。
如图3所示,全局寻优阶段:对目标函数以S6.1寻优结果为基础进行全局寻优包括以下步骤:
S6.2.1:初始值设置,本步骤的初始值以S6.1中的初步寻优输出结果作为起点。经过初步寻优阶段的预处理,获得了一个较为可信的初始值,将加快本阶段算法的收敛速度。禁忌表的初值则全部初始化为预处理阶段输出的最优解;
S6.2.2:邻域点设置:本步骤的邻域点设置依然利用目标函数的差分来指导邻域点的选取。但是有别于预处理阶段。具体体现在取值范围不同。在本阶段,每次得到未被禁忌的伪最优解时,取值空间就会缩小为原空间的1/2。利用多维正态分布函数来选取邻域点。其中该多维正态分布函数的标准差设置方法同预处理阶段。通过本阶段的搜索,基本可以确定最优解的大致范围。
如图4所示加速寻优阶段:对目标函数以S6.2的寻优结果为基础进行快速搜索包括如下步骤:
S6.3.1:初始值设置:本阶段的初始值以S6.2中的全局寻优输出结果作为起点。禁忌表初值则为S6.2中的全局寻优得出的最优解;
S6.3.2:邻域点设置,在上述S6.1初步寻优阶段和S6.2全局寻优阶段的寻优基础之上,寻优的范围缩小为S6.2全局寻优阶段中最小取值范围。在本取值范围内,采用定步长搜索策略。步长初始值为整个取值域的1/n,其中n是设定值,n越大,步长越小。定步长搜索的耗时短,并且在所限制的范围内小步搜索,可以提高解的精度。
优化过程如下:
邻域设置:取5个邻域点。其选取过程为计算每个变量在当前最优解处的差分值,以该差分为基础,经过数学变换产生一个适当的标准差。以标准差为方差,当前最优解为均值确定一个多维正态分布函数,通过该函数来选取5个邻域点。
禁忌表内容:禁忌表内容为前5次所找到的最优解的坐标、目标函数值和禁忌长度。
特赦规则:本算法以禁忌对象的禁忌长度为特赦规则。即如果某一个禁忌对象其禁忌长度超过5,则将该禁忌对象特赦。
终止规则:
1.算法执行到规定的最大迭代次数
2.在一定迭代次数内解无明显变化
3.在一定迭代次数内未更新禁忌表
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (3)
1.一种基于分阶段禁忌寻优算法的土壤水分特征参数的计算方法,其特征在于,包括:
S1:通过实验测量获取土壤饱和导水率和饱和含水量;
S2:根据S1中的土壤饱和导水率和饱和含水量,计算出土壤非饱和导水率和土壤水分扩散率;
S3:建立一维垂直土壤水分运动数学模型;
S4:将S2中的土壤非饱和导水率和土壤水分扩散率代入到S3的模型中,转化为土壤含水率理论值计算公式θ(a,b,m);
S5:根据入渗实验,确定土壤含水率的实际测量值θ;
S6:建立最优化目标函数,并将S4和S5中的θ(a,b,m)和θ代入目标函数中,根据改进的分段禁忌算法进行优化,以确定出a、b、m的最优值;所述改进的分段禁忌算法为以基本禁忌算法为基础,分别在三个阶段中的每个阶段的参数进行对应的设置和选取;
S7:将S6中的a、b、m最优值代入S3中的一维垂直土壤水分运动数学模型中计算出土壤水分特征参数;
其中,在S6中利用改进的分段禁忌算法对目标函数进行参数寻优,所述改进的分段禁忌算法具体流程步骤包括:
S6.1:初步寻优阶段:对目标函数进行初步寻优;
S6.2:全局寻优阶段:对目标函数以S6.1寻优结果为基础进行全局寻优;
S6.3:加速寻优阶段:对目标函数以S6.2的寻优结果为基础进行快速搜索;
其中,所述S6.1具体包括:
S6.1.1:初始值设置,初始值以随机方式在整个取值空间内取得;禁忌表为空,并采用单向列表方式加以管理;
S6.1.2:邻域点的选取,利用目标函数的差分来指导邻域点的选取;
其中,所述S6.1.2具体包括:
S6.1.2.1计算目标函数对每一个变量的差分及其在当前解的差分值;
S6.1.2.2在整个取值范围内,利用多维正态分布函数来选取邻域点。
2.如权利要求1所述计算方法,其特征在于,所述S6.2具体包括:
S6.2.1:初始值设置,本步骤的初始值以S6.1中的初步寻优输出结果作为起点;
S6.2.2:邻域点设置,本步骤的邻域点设置依然利用目标函数的差分来指导邻域点的选取。
3.如权利要求1所述计算方法,其特征在于,所述S6.3具体包括:
S6.3.1:初始值设置:本阶段的初始值以S6.2中的全局寻优输出结果作为起点,禁忌表初值则为S6.2中的全局寻优得出的最优解;
S6.3.2:邻域点设置,在上述S6.1初步寻优阶段和S6.2全局寻优阶段的寻优基础之上,寻优的范围缩小为S6.2全局寻优阶段中最小取值范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210208585.7A CN102692491B (zh) | 2012-06-19 | 2012-06-19 | 基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210208585.7A CN102692491B (zh) | 2012-06-19 | 2012-06-19 | 基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102692491A CN102692491A (zh) | 2012-09-26 |
CN102692491B true CN102692491B (zh) | 2014-07-30 |
Family
ID=46858086
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210208585.7A Expired - Fee Related CN102692491B (zh) | 2012-06-19 | 2012-06-19 | 基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102692491B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5858080B2 (ja) * | 2013-08-23 | 2016-02-10 | 横河電機株式会社 | 運転計画策定方法および運転計画策定システム |
CN104198317B (zh) * | 2014-09-03 | 2016-06-01 | 塔里木大学 | 测定土壤物理参数s的方法 |
CN104764934B (zh) * | 2015-04-27 | 2019-04-26 | 塔里木大学 | 土壤饱和浸提液电导率的测量方法 |
CN105022922B (zh) * | 2015-07-16 | 2018-01-16 | 中国水利水电科学研究院 | 基于水循环全过程非饱和带土壤水分消耗特性的评价方法 |
CN109857937A (zh) * | 2019-01-25 | 2019-06-07 | 同济大学 | 一种基于Sine映射和分段式并行连续禁忌搜索数据处理方法 |
CN110235756A (zh) * | 2019-07-11 | 2019-09-17 | 中国水利水电科学研究院 | 一种水稻全生育期灌溉水量及时间的确定方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763601A (zh) * | 2010-01-12 | 2010-06-30 | 武汉大学 | 基于禁忌搜索算法的土地利用分区方法 |
CN102289562A (zh) * | 2011-06-08 | 2011-12-21 | 南京航空航天大学 | 基于禁忌搜索算法的压电作动器方向优化配置方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7085690B2 (en) * | 2000-06-10 | 2006-08-01 | Mark Edward Sale | Unsupervised machine learning-based mathematical model selection |
-
2012
- 2012-06-19 CN CN201210208585.7A patent/CN102692491B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101763601A (zh) * | 2010-01-12 | 2010-06-30 | 武汉大学 | 基于禁忌搜索算法的土地利用分区方法 |
CN102289562A (zh) * | 2011-06-08 | 2011-12-21 | 南京航空航天大学 | 基于禁忌搜索算法的压电作动器方向优化配置方法 |
Non-Patent Citations (11)
Title |
---|
A Staged Continuous Tabu Search Algorithm for the Global Optimization and its Applications to the Design of Fiber Bragg Gratings;Zheng R.T.等;《Computational Optimization and Applications》;20051231;第30卷(第3期);第319-335页 * |
Estimating unsaturated soil hydraulic parameters using ant colony optimization;K.C. Abbaspour等;《Advances in Water Resources》;20010831;第24卷(第8期);第827-841页 * |
K.C. Abbaspour等.Estimating unsaturated soil hydraulic parameters using ant colony optimization.《Advances in Water Resources》.2001,第24卷(第8期),第827-841页. |
Zheng R.T.等.A Staged Continuous Tabu Search Algorithm for the Global Optimization and its Applications to the Design of Fiber Bragg Gratings.《Computational Optimization and Applications》.2005,第30卷(第3期),第319-335页. |
周明丽等.土壤水分运动参数研究.《现代农业科技》.2009,(第4期),第136-138页. |
唐亚莉等.优化算法确定土壤水分特征曲线的分析.《新疆大学学报(自然科学版)》.2006,第23卷(第2期),第240-243页. * |
土壤水分运动参数研究;周明丽等;《现代农业科技》;20090220(第4期);第136-138页 * |
基于连续函数优化的禁忌搜素算法;张晓菲等;《中国计量学院学报》;20100915;第21卷(第3期);第251-256页 * |
基于遗传算法的土壤水分运动参数识别;杨坤等;《农业工程学报》;20090930;第25卷(第9期);第32-35页 * |
张晓菲等.基于连续函数优化的禁忌搜素算法.《中国计量学院学报》.2010,第21卷(第3期),第251-256页. |
杨坤等.基于遗传算法的土壤水分运动参数识别.《农业工程学报》.2009,第25卷(第9期),第32-35页. |
Also Published As
Publication number | Publication date |
---|---|
CN102692491A (zh) | 2012-09-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102692491B (zh) | 基于分阶段禁忌寻优算法的土壤水分特征参数计算方法 | |
CN109829217B (zh) | 压裂性裂缝油藏产能模拟方法及装置 | |
CN105353620A (zh) | 一种基于地表水深信息的地面灌溉控制方法 | |
CN104598565B (zh) | 一种基于随机梯度下降算法的k均值大规模数据聚类方法 | |
CN110941890B (zh) | 基于最优控制理论的海上油藏动态实时生产优化方法 | |
CN102262692B (zh) | 飞机翼面蒙皮亚音速颤振优化方法 | |
CN104179499A (zh) | 考虑油藏参数时变的数值模拟方法 | |
CN107607692B (zh) | 基于土壤最大储水容量的土壤水分监测优化布点方法 | |
CN103106256A (zh) | 一种基于Markov链的正交化插值的GM(1,1)模型预测方法 | |
CN111222271A (zh) | 基于基质-裂缝非稳态窜流油藏裂缝数值模拟方法及系统 | |
CN103149600A (zh) | 一种基于优化控制点确定水文地质参数的自动配线方法 | |
CN103310284B (zh) | 经济负荷和降温负荷的确定方法及负荷确定装置 | |
CN102682172A (zh) | 基于参数分类的超临界机翼极多参数优化设计方法 | |
CN103530700B (zh) | 城区配电网饱和负荷综合预测方法 | |
CN101739718B (zh) | 一种基于参数模板的玉米叶片虚拟仿真建模方法 | |
CN114676522B (zh) | 融合gan和迁移学习的气动形状优化设计方法及系统及设备 | |
CN109800454A (zh) | 一种基于卡尔曼滤波的煤层气采集方法、系统及存储介质 | |
CN110765677B (zh) | 一种高精度、快速三维地质模型有限元模型的建模方法 | |
CN106096183B (zh) | 一种基于特征线法的多重并行方法 | |
CN102937946B (zh) | 一种基于约束正则模式的复杂函数极小值搜索方法 | |
CN105373846A (zh) | 基于分级策略的油气集输管网拓扑结构智能优化方法 | |
CN103593534A (zh) | 基于工程地质因素关联性的盾构机智能选型方法及装置 | |
CN103984832A (zh) | 一种铝电解槽电场仿真分析方法 | |
CN105490836B (zh) | 一种复杂网络可靠度的蒙特卡罗评估方法 | |
CN103970610A (zh) | 一种供水管网中节点流量的监控方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140730 Termination date: 20150619 |
|
EXPY | Termination of patent right or utility model |