Study on Leachate Recirculation Volume of Aerobic Bioreactor Landfills
-
摘要: 为了研究渗滤液回灌量对好氧填埋场稳定化进程的影响，通过模拟实验，探讨了不同回灌量时，好氧填埋场固相垃圾与渗滤液特性的变化趋势，填埋周期为122 d. 结果表明，不同的渗滤液回灌下各反应器中填埋垃圾的含水率始终保持在70%左右；填埋过程中，回灌量为20%的反应器填埋垃圾中总有机质含量、（半纤维素 + 纤维素）/木质素比（Q）值与沉降量变化速度最快，至填埋结束时，总有机质含量与Q值较其他反应器低8.9%～14.6%和9.9%～16.9%，沉降量较其他反应器提高了6.6%～13.3%；在20%回灌量下，整个填埋周期填埋场所产渗滤液中的还原性有机物的量较10%、15%、30%和全回灌分别低9.8%、12.5%、17.8%和14.9%，氨氮在52 d，即降至25 mg/L，达到GB 16889—2008所规定的垃圾填埋场渗滤液氨氮排放浓度值，较其他反应器提前7～21 d.Abstract: To study the influence of leachate recirculation on the stabilization process of aerobic landfills, an experiment was conducted to assess the variation trend of solid waste and leachate characteristics in an aerobic landfill with different recirculation rates. The landfill period was 122 days. Results showed that the water content of landfill refuse in each reactor remained at about 70% under different leachate recirculation rates. The total organic carbon content, （hemicellulose + cellulose）/lignin ratio（Q）, and sedimentation of landfill refuse with 20% recirculation changed fastest in the landfill process. At the end of the landfill period, the mass content and Q of landfill leachate were 8.9%–14.6% and 9.9%–16.9% lower, respectively, than that of other reactors, and the sedimentation was 6.6%–13.3% higher than that of other reactors. The amount of reductive organic matter in the leachate during the whole process with the recirculation volume as 20% was 9.8%, 12.5%, 17.8%, and 14.9% lower than that of 10%, 15%, 30%, and full recirculation, respectively, and the ammonia nitrogen was reduced to 25 mg/L, which met the concentration of ammonia nitrogen discharge in landfill leachate of GB 16889—2008, in 52 d, which was 7–21 d earlier than the other reactors.
-
Key words:
- aerobic landfill /
- recirculation volume /
- stabilization /
- biodegradation /
- leachate
-
截止2019年底，我国高速铁路营业里程已达到3.5万公里，位居世界第一. 我国已成为世界上高铁运营里程最长、运输密度最高、运营场景最复杂的国家. 铁路信号系统软件测试是保障高铁信号系统安全运行的重要方法，能够有效地预防列车追尾和出轨等行车事故，避免造成人员伤亡及重大财产损失[1]，对保障我国高铁的安全运行具有重要意义.
目前，高铁信号系统软件测试的研究热点之一为基于形式化模型的测试案例自动编制方法，其首先建立软件的测试模型，然后研究算法自动生成满足特定目标的测试案例. 铁路信号系统软件测试中常用的建模方法包括UML状态图（UML statechart）、时间自动机（timed automata，TA）和有色Petri网（colored Petri net，CPN）等. 文献[2]以CPN为测试建模理论基础，以无线闭塞中心切换功能为例，研究基于全路径覆盖和序列优选的测试案例自动编制方法，降低测试案例的重复度. 文献[3]以CTCS-3 （Chinese Train Control System 3）级ATP （automatic train protection）设备模式转换功能为例，基于CPN建模理论，解决状态空间爆炸等问题，生成满足全路径覆盖准则的测试序列. 文献[4]以TA为基础进行测试建模，提出一种满足全状态和全变迁覆盖准则的测试案例生成方法，提高测试案例的生成效率和重用性. 文献[5]研究CTCS-3故障特点，提出了一种基于输入输出时间自动机变异模型的测试案例生成方法. 文献[6]以轨道电路编码功能为研究对象，提出图覆盖、组合覆盖和文法分析相结合的方式自动从UML状态图模型生成测试案例. 文献[7]以UML顺序图为研究基础，利用深度优先算法和网络组合算法对计算机联锁系统测试案例进行组合优化，降低测试案例的重复度. 文献[8]针对铁路信号系统软件的测试度量方法，提出对信号系统软件进行风险分级，按照风险等级递减顺序进行安全性测试的测试策略，指出了以铁路信号系统软件风险等级为基础测试的必要性. 文献[9]对CTCS-3系统安全软件进行接口危险辨识，评估危险事件的风险等级，提出列控系统安全软件黑盒测试的危险分析方法，但未对软件测试建模理论进行研究. 文献[10-12]针对铁路信号系统的形式化描述和验证，分别从时钟约束、风险等级和安全约束等几个方面提出铁路信号系统的建模方法，但这些方法结构复杂，测试模型解析和测试案例搜索算法难以展开，不易直接应用于铁路信号系统软件测试建模.
目前，铁路信号系统软件测试的研究重点是如何以既有的测试建模方法为基础生成覆盖全面且高效的测试案例，研究提高了测试活动的自动化程度和测试效率. 但现有测试建模方法不能完整地描述软件所有的行为特征，如UML状态图和CPN等方法不能描述软件的时间约束，导致测试案例不能反映软件的时间特性. TA能够描述软件的时间约束，但不具有层次结构，对复杂功能逻辑的描述能力较差，容易出现模型难以理解，甚至错误的问题. 同时，UML状态图、CPN和TA等方法均不能描述软件测试模型中的风险特性. 本文通过分析铁路信号软件的领域特点，总结铁路信号软件测试的建模需求，结合Z规格说明语言（简称Z语言）提出风险时间状态机建模方法，满足铁路信号软件测试的建模需求.
1. 建模概述
铁路信号系统软件测试通常采用基于模型的测试（model-based testing，MBT）[13]理论，主要包括分析软件需求、建立测试模型、生成测试案例和执行测试4个阶段，如图1所示. 首先，基于被测软件的需求规格等文件对软件功能和结构进行抽象，采用形式化建模方法建立软件的测试模型；然后，依据测试模型，采用“平坦化”和状态组合等方法对测试模型进行解析，通过广度优先等搜索算法自动生成测试案例；最后，将测试案例加载到测试环境中执行测试，观察测试结果.
测试模型是编制测试案例和评判测试结果时的重要依据，直接关系到软件测试过程的全面性和有效性，铁路信号软件测试建模方法需要尽可能系统、全面地描述铁路信号系统软件的所有需求.
2. 建模需求
铁路信号系统软件测试建模包括功能建模需求和性能建模需求两个方面. 功能建模需求描述系统的功能逻辑（function logic），性能建模需求指与系统相关的性能约束或限制，包括时钟约束（clock constraint）和风险等级（risk level）两个方面.
1） 功能逻辑
随着铁路信号系统的规模越来越大，软件承担的功能越来越多、逻辑越来越复杂. 如CTCS-3级列控系统包括9种工作模式，14个主要运营场景，206个功能特征. 铁路信号系统已构成一个复杂的控制系统，软件存在着大量的并发、竞争、冲突等逻辑关系，控制状态转移条件复杂.
2） 时钟约束
铁路信号系统软件运算结果的正确性不仅取决于系统逻辑处理的正确性，还取决于运算过程中的时钟约束，要求软件功能需要在指定时间内完成，或在规定时限之后才能发生. 如《CTCS-3级列控系统总体技术方案》[14]要求“当列车前端距分相区还有10秒走行距离时，车载设备向司机发出提示”. 铁路信号系统软件的时钟约束通常为硬实时性（hard real-time）要求，如道岔转换时间、移动授权有效时间、目视模式确认时间等，时钟约束错误或缺失，可能造成重大安全事故.
3） 风险等级
铁路信号系统软件由安全功能（safety function）、安全相关功能（safety-related function）和非安全功能（non-safety function）组合构成. 如 《CTCS-3级列控系统测试案例（V3.0）》[15]的功能特征CTCS3-FT-71“由于超速触发紧急制动”属于安全功能，功能特征CTCS3-FT-116“记录数据的下载”属于非安全功能. 不同安全级别的功能对系统造成的安全影响差别较大，由风险等级进行描述. 风险等级的高低是铁路信号系统软件测试案例编制的重要基础，对组织测试案例优先级及分析测试需求具有重要意义.
根据以上分析，铁路信号系统软件测试建模方法需要描述系统复杂的功能逻辑、时钟约束和风险等级3个方面的特点，严密、准确地反应软件测试需求，且具有良好的可读性和可理解性.
3. 风险时间状态机
为定义风险时间状态机建模方法，首先介绍有限状态机（finite state machine, FSM），然后在FSM基础上，结合Z语言，在迁移元素中扩展出时钟约束，在状态元素中扩展出风险等级，最后给出风险时间状态机的定义.
3.1 有限状态机概述
FSM采用可视化的方式描述系统的功能逻辑，具有清晰、直观的优点，其是表示有限状态以及状态之间转移和动作等行为的数学模型，具有精确性、可推导性和可验证性，是基于模型的测试理论中的经典建模方法[13-16].
定义1 一个有限状态机M是一个六元组，如式（1）[17].
M=(S,s0,δ,λ,I,O), (1) 式中：
S为有限状态集合；
$ {s_0} $ 为初始状态，且$ {s_0} \in S $ ；$ \lambda $ 为状态转移函数，$ \lambda :S \times I \to S $ ；$\delta $ 为输出函数，$ \delta :S \times I \to O $ ；I 为有限输入符号集合；
O 为有限输出符号集合.
图2为一个FSM模型，包括s1和s2两个状态，其中s1为初始状态，A、B、C为输入，a、b、c为输出. s1接收到输入A时输出a，转移到s2状态. 模型对应的六元组如下：
$ S = \left\{ {{s_{\text{1}}},{s_{\text{2}}}} \right\} $ ，$ {s_0} = {s_{\text{1}}} $ ，$ I = \left\{ {A,B,C} \right\} $ ，$ O = \left\{ {a,{\kern 1pt} {\kern 1pt} b,{\kern 1pt} {\kern 1pt} c} \right\} $ ，$ \lambda \left( {{s_{\text{1}}},A} \right) = {s_{\text{2}}} $ ，$ \lambda \left( {{s_{\text{2}}},B} \right) = {s_{\text{1}}} $ ，$ \lambda \left( {{s_{\text{2}}},C} \right) = {s_{\text{2}}} $ ，$ \delta \left( {{s_{\text{1}}},A} \right) = a $ ，$ \delta \left( {{s_{\text{2}}},B} \right) = b $ ，$ \delta \left( {{s_{\text{2}}},C} \right) = c $ .接下来结合Z语言，在FSM描述功能逻辑基础上，增加层次结构，扩展出时钟约束和风险等级两个参数，提出风险时间状态机（risk timed statechart，RTSC）. 为描述方便，首先定义RTSC的基本元素，然后定义RTSC的层次结构，最后给出RTSC的形式化定义.
3.2 RTSC的基本元素
RTSC包括状态、时钟、信号表达式和迁移等基本元素.
1） 状态
状态是铁路信号系统软件在一定时期内的存在形式. 对RTSC的非空、有限状态集合S，Z语言描述如下：
[S]
∃m:N1•#S=m 状态包括简单状态（SIMPLE），或状态（OR）和与状态（AND）3种类型，Z语言描述如下：
P :: = SIMPLE | OR | AND
OR状态由同一层次的状态和迁移组成，包含且仅包含一个初始状态. 初始状态无源状态，且是任意状态的源状态. 多个并发的OR状态组成AND状态.
2） 时钟
RTSC采用时钟描述时间的流逝，有限时钟集合记为X，Z语言描述如下：
[X]
∃m:N•#X=m 对一个时钟变量集合X，时间约束集合定义如式（2）.
Σ::=x∝g|σ1∧σ2， (2) 式中：
$ {\sigma _1} $ 、$ {\sigma _{\text{2}}} $ 为时钟约束；$ x \in X $ ；$ g \in {\text{N}} $ ；$\propto \in \{ \leqslant , < , \geqslant , > , = \}$ .RTSC时钟约束Σ的Z语言描述如下：
[Σ]
∃m:N•#Σ=m 3） 信号
信号E具有产生和不产生两种状态，Z语言描述如下：
[E]
∃m:N•#E=m 在某一时刻，正交组件产生的信号称为动作.
4） 信号表达式
信号之间的与（and）关系、或（or）关系、非（not）关系构成信号表达式，记为ex如not e表示未产生信号e.
用Π为表示信号永远发生的常量.
5） 迁移
迁移是RTSC从源状态转移到目标状态的方式，包括源状态sc、时钟约束σ、信号表达式ex、转移时产生的信号e、重置的时钟x及目标状态st. 迁移的Z语言描述如下：
3.3 RTSC的层次结构
RTSC的状态层次H包括底状态γ、为状态分配子状态的有限状态层次函数ω，定义状态类型的有限状态类型函数κ，定义状态父状态的ρ函数. H满足以下性质：
1） ω为非SIMPLE类型的状态分配子状态；
2） 状态之间的层次关系不能形成循环结构；
3） γ为OR状态或AND状态，是唯一没有父状态的状态；
4） 任意非γ的状态都具有唯一的父状态，且均可以由γ通过ω传递到达.
H的Z模式定义如下：
其中：
$ \mathbb{F}$ S表示集合S的所有有限子集的集合；dom ω表示ω的定义域；ran ρ表示ρ的值域；$ \forall $ si: dom ρ; sj: dom ω表示ρ定义域中的任意状态si和ω定义域中的任意状态sj， i，j为不为0的自然数，且i ≠ j.RTSC由状态层次和迁移集构成，采用风险等级函数η为SIMPLE状态分配风险等级，满足以下性质：
1） η仅为SIMPLE状态分配风险等级；
2） OR状态和AND状态的风险等级为其包含的子状态的风险等级的最大值.
3） 对任意迁移，只能有一个同层次的源状态和目标状态，源状态和目标状态可以是AND或SIMPLE状态.
RTSC的Z模式定义如下：
其中：
$\mathbb{F}_1$ T表示集合T的所有非空有限子集的集合.4. RTSC的格局转移机制
4.1 格局
一个时钟集合X的时钟解释v是指每个时钟变量x到时间序列上的一个全映射. 某时刻，RTSC能够同时处于最大的状态集结合当前的时钟解释称为格局（configuration），记为U. 在任意时刻，RTSC只有一个活动的格局，满足以下规则：
1） U包含γ状态；
2） U包含AND状态s，则U 包含s的每一个子状态；
3） U包含OR状态s，则U包含s的某一个子状态；
4） U包含非γ的状态s，则U包含s的父状态；
5） U仅包含满足规则1） ~ 4）的所有状态.
格局公理描述如下：
其中：k表示不为0的自然数.
4.2 格局转移机制
RTSC格局转移过程记为 (us, v(X), E)，其中us表示当前格局包含的状态集，v(X)表示时钟解释，E表示转移过程中产生的动作集，具体过程如下：
1） 激励. 当RTSC接受到新的时钟解释或信号时开始格局转移：
$\left( {{u_{\rm{s}}},v\left( X \right),\varnothing } \right) \to \left( {{u_{\rm{s}}},v'\left( X \right),e} \right)$ .2） 计算. 格局转移过程以SIMPLE状态为基础，从γ状态，通过OR状态和AND状态向下分解，包括以下2个规则：
规则1 AND状态并行地计算其所包含的OR状态的转移.
规则2 OR状态计算当前激活的SIMPLE状态的变迁t上的信号，产生t上的动作，重置t的时钟，离开变迁源状态τ.sc到达变迁状态t.st，满足：
${u_{\rm{s}}}{\kern 1pt} {\kern 1pt} {\kern 1pt} = {\kern 1pt} {\kern 1pt} {\kern 1pt} {\kern 1pt} {u_{\rm{s}}}{\kern 1pt} {\kern 1pt} \backslash {\kern 1pt} {\kern 1pt} t.{s_{{\rm{c}}}} \cup t. {s_{\rm{t}}}$ ，$G = G \cup t.e$ ，$v\left( {t. x} \right) = 0$ .3） 结束. RTSC根据状态类型迭代相应规则直至转移稳定，保留时钟解释，清空转移期间的信号.
RTSC从初始格局开始，根据产生的信号或时钟的更新激活格局转移过程；依据层次结构，迭代AND状态和OR状态的规则，直到计算出下一组稳定的状态. 格局转移时，迁移的触发信号均发生，且时钟解释满足迁移上的时钟约束，同时产生迁移上的动作，重置迁移上的时钟变量.
4.3 RTSC的操作语义
RTSC的语义可以由其对应的时间标记迁移系统（timed transition system，TTS）进行定义.
定义2 一个时间标记迁移系统TTS是一个4元组：
Ω=(L,l0,→,Σ), (3) 式中：L为状态集合；
l0为初始状态，l0∈L；
$\to \subseteq L \times \left( {\varSigma \cup \mathbb{R}} \right) \times L$ 为转移关系，且满足：1） 如果$ l\xrightarrow{0}l' $ ，则$ l{\kern 1pt} {\kern 1pt} {\text{ = }}{\kern 1pt} {\kern 1pt} l' $ ； 2） 如果$ l\xrightarrow{d}l' $ ，且$ l'\xrightarrow{{d'}}l'' $ ，$ d,d' \in \mathbb{R} $ ，则$ l\xrightarrow{{d{\text{ + }}d'}}l'' $ ； 3） 如果$ l\xrightarrow{d}l' $ ，$ d \in \mathbb{R} $ ，则$ {\text{0}} \leqslant d' \leqslant d $ ，则有$ l'' \in L $ ，满足$ l\xrightarrow{{d'}}l'' $ 且$ l''\xrightarrow{{d - d'}}l' $ .在时钟解释v下，如果时钟约束σ为真，则称时钟解释v满足时钟约束σ，记为
$ v| = \sigma $ .定义3 对任意RTSC，其语义可以通过时间标记迁移系统
$\varOmega = {\kern 1pt} {\kern 1pt} \left( {L,{l_0}, \to ,\varSigma } \right)$ 进行定义，其中：$ L = U \times {\mathbb{R}^X} $ ；${l_0}{\text{ = }}\left( {{u_{{\rm{s}}0}},{v_0}} \right)$ ，${u_{{\rm{s}}0}}$ 为RTSC的初始状态集，且对任意$ x \in X $ ，$ {v_0}\left( x \right) = 0 $ ；转移关系
$ \to $ 满足：1） 对$\left( {{u_{\rm{s}}},v} \right)$ 和一个迁移${u_{\rm{s}}}\xrightarrow{{{e_x},\sigma ,e,x}}{u_{\rm{s}}}^\prime$ ，且$ v| = \sigma $ ，则$\left( {{u_{\rm{s}}},v} \right)\xrightarrow{{{e_x}}}\left( {{u_{\rm{s}}}^\prime ,v'} \right)$ ；2） 对$\left( {{u_{\rm{s}}},v} \right)$ 和时间增量$ d \in \mathbb{R} $ ，$\left( {{u_{\rm{s}}},v} \right)\xrightarrow{d}\left( {{u_{\rm{s}}},v + d} \right)$ .RTSC从初始格局开始，通过动作转移或时延转移实现格局变化.
1） 动作转移. 当前激活的信号和时钟解释满足变迁的信号表达式和时钟约束，触发动作迁移，激活变迁上的动作，重置变迁时钟为0.
2） 时延转移. 时延转移保持格局状态不变，对所有时钟增加相同的时间增量.
5. 案例分析
计算机联锁系统是铁路信号系统中典型的安全关键系统. 以计算机联锁系统的道岔定位选排子系统为例建立RTSC测试模型，解释RTSC的格局转移机制.
5.1 基于RTSC的道岔转换测试模型
道岔定位选排子系统根据道岔操作命令或进路请求，检查道岔当前是否处于定位. 当道岔状态与期望不一致时，选排道岔到需求位置[18]. 该系统涉及11个信号，主要信号名称及其含义见表1.
表 1 道岔子系统信号含义Table 1. Meaning of signals in switch subsystem信号 含义 信号 含义 信号 含义 e1 定位请求 e4 定位需求 e7 道岔定位 e2 道岔转换 e5 道岔锁闭 e8 转换成功 e3 转换超时 e6 反位请求 e9 选排一致 道岔定位选排子系统的RTSC模型如图3所示，记为DRTSC，其中ran η = {0, 1, 2, 3, 4}，分别对应EN 50128定义的5个安全完整度等级SIL 0～4. DRTSC包括1个时钟x1，最高风险等级为SIL 4，最低为SIL 2，底状态为s0，初始状态为
$ \left\{ {{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{6}}},{s_{\text{9}}},{s_{{\text{11}}}}} \right\} $ ，s1、s2、s3分别为SIMPLE状态、AND状态和OR状态，$ \eta \left( {{s_1}} \right) = {\text{4}} $ ，$ \eta \left( {{s_{\text{3}}}} \right) = \max \left( {\eta \left( {{s_{\text{6}}}} \right),\eta \left( {{s_{\text{7}}}} \right),\eta \left( {{s_{\text{8}}}} \right)} \right) = {\text{4}} $ ，$ \rho \left( {{s_2}} \right) = {s_0} $ ，$ {\kern 1pt} {\kern 1pt} \omega \left( {{s_2}} \right){\text{ = }}\left\{ {{s_3},{s_4},{s_5}} \right\} $ .DRTSC中：迁移的表现形式为[e] [σ] / [a] [x]；e为触发转移的信号；σ为时钟约束；a为转移时产生的动作；x为重置的时钟，记为
$ {s_1}\xrightarrow{{e,{\kern 1pt} \sigma {\kern 1pt} /{\kern 1pt} {\kern 1pt} a,{\kern 1pt} x}}{s_{\text{2}}} $ . 如$ {s_{14}}\xrightarrow{{{e_7},{\kern 1pt} {x_{\text{1}}} < 13{\kern 1pt} {\kern 1pt} /{\kern 1pt} {\kern 1pt} {e_8},{\kern 1pt} {x_{\text{1}}}}}{s_{13}} $ 表示状态s14激活时，当检测到信号e7，并满足时间约束x1<13，DRTSC迁移到s13，同时产生信号e8，并重置时钟x1.5.2 DRTSC格局转移
DRTSC模型中，s3、s4、s5分别表示道岔转换的“安全请求”，“道岔转换”和转换的“时钟约束”状态. DRTSC主要描述以下过程：
道岔“安全请求”状态检测到定位需求后，在道岔未锁闭、无反位请求、未处于定位时，进一步检查道岔转换状态，未处于转换过程中或转换超时时，判定定位请求安全. “道岔转换”状态接受到“定位请求”信号后，如果道岔未锁闭则重置时钟x1；条件保持1 s后，请求道岔转换. “时钟约束”状态接受道岔转换动作后，重置时钟x1并开始计时，13 s内检测到道岔转换到定位则转换成功，否则转换超时. 当“道岔转换”状态接受到“时钟约束”状态的“转换成功”信号后判定道岔转换成功. 当同时处于“定位请求”和“转换成功”状态时，判定“选排一致”. 转换过程中任意时刻，道岔子系统接受到反位请求信号立即重置时钟，返回初始状态.
DRTSC的初始格局为
$( ( {s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}}, {s_{\text{6}}},{s_{\text{9}}},{s_{{\text{13}}}} ), ( {\text{0}} ) )$ . 下面以格局$ \left( {\left( {{s_{\text{0}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}},{s_{{\text{11}}}},{s_{{\text{14}}}}} \right){\text{,}}\left( {\text{2}} \right)} \right) $ ，产生信号e7，$ v\left( x \right) = {\text{3}} $ 为例，分析DRTSC的格局转移过程.1） DRTSC接受激励信号和更新的时钟解释：
$\left( {\left( {{s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}},{s_{{\text{11}}}},{s_{{\text{14}}}}} \right),\left( {\text{2}} \right),\left( \varnothing \right)} \right) \to$ $( ( {s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}}, $ ${s_{\text{8}}},{s_{{\text{11}}}},{s_{{\text{14}}}} ){\text{,}}\left( {\text{3}} \right){\text{,}}( {{e_{\text{7}}}} ) ) $ .2） 由于s0为OR状态，根据规则 2，s0当前激活状态为s2. s2为AND状态，根据规则 1，s2并行计算
$ {s_{\text{3}}} $ ,$ {s_{\text{4}}} $ ,$ {s_{\text{5}}} $ ，计算过程为e7和
$ v\left( x \right) = {\text{3}} $ 满足$ {s_{{\text{14}}}} $ 迁移到$ {s_{{\text{13}}}} $ 的信号表达式和时钟约束，迁移时重置时钟x1，同时产生信号e8:$ \left( {\left( {{s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}},{s_{{\text{11}}}},{s_{{\text{14}}}}} \right){\text{,}}\left( {\text{3}} \right){\text{,}}\left( {{e_{\text{7}}}} \right)} \right) \to $ $( ( {s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}}, $ $ {s_{\text{8}}},{s_{{\text{11}}}},{s_{{\text{13}}}} ){\text{,}}\left( {\text{0}} \right){\text{,}}( {{e_{\text{7}}},{e_{\text{8}}}} ) ); $ e8信号触发s4的子状态s11，并转移到状态s12：
$ \left( {\left( {{s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}},{s_{{\text{11}}}},{s_{{\text{13}}}}} \right){\text{,}}\left( {\text{0}} \right){\text{,}}\left( {{e_{\text{7}}},{e_{\text{8}}}} \right)} \right) \to $ $( ( {s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}}$ ,${s_{\text{5}}},{s_{\text{8}}},{s_{{\text{12}}}},{s_{{\text{13}}}} ){\text{,}}\left( {\text{0}} \right){\text{,}} \left( {{e_{\text{7}}},{e_{\text{8}}},{e_{{\text{11}}}}} \right) ); $ s3的状态s8满足自迁移变迁：
$( ( {s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}}, {s_{{\text{12}}}},{s_{{\text{13}}}} ),\;\left( {\text{0}} \right),\;\left( {{e_{\text{7}}},{e_{\text{8}}},{e_{{\text{11}}}}} \right) ) \to$ $( ( {{s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}},{s_{{\text{12}}}},{s_{{\text{13}}}}} ) ,$ $( {\text{0}} ){\text{,}}( {{e_{\text{7}}},{e_{\text{8}}},{e_{{\text{11}}}},{e_{\text{1}}}} ) ). $ s3，s4和s5计算结束后，DRTSC激活的信号包括
$( {e_{\text{7}}}, {e_{\text{8}}},{e_{{\text{11}}}},{e_{\text{1}}} )$ ，激活s2至s1的迁移：$( ( {s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}}, {s_{{\text{12}}}}, {s_{{\text{13}}}} ){\text{,}}\left( {\text{0}} \right){\text{,}}\left( {{e_{\text{7}}},{e_{\text{8}}},{e_{{\text{11}}}},{e_{\text{1}}}} \right) ).$ $ \to \left( {\left( {{s_{\text{0}}},{s_{\text{1}}}} \right){\text{,}}\left( {\text{0}} \right){\text{,}}\left( {{e_{\text{7}}},{e_{\text{8}}},{e_{{\text{11}}}},{e_{\text{1}}}} \right)} \right) $ .3） 计算结束，DRTSC格局更新为
$ \left(\left({s}_{\text{0}}，{s}_{\text{1}}\right)\text{,}\left(\text{0}\right)\right) $ .根据以上分析，DRTSC格局转移过程记为
$\left( {\left( {{s_{\text{0}}},{s_{\text{2}}},{s_{\text{3}}},{s_{\text{4}}},{s_{\text{5}}},{s_{\text{8}}},{s_{{\text{11}}}},{s_{{\text{14}}}}} \right),\left( {\text{2}} \right)} \right)\xrightarrow{{{e_7},v\left( x \right)\; = \;3}}$ $\left(\left({s}_{\text{0}}， {s}_{\text{1}}\right),\left(\text{0}\right)\right),$ 表示道岔在处于“定位安全请求”，“道岔转换”，“道岔转换计时中”状态时，检测到道岔在第3 s转移到定位，判定道岔选排一致.5.3 DRTSC模型分析
DRTSC描述了道岔定位转移过程中的功能逻辑、时钟约束和风险等级.
5.3.1 测试需求覆盖
1） 功能逻辑
DRTSC采用16个状态，19条迁移和11个信号描述道岔转换过程的功能逻辑.
2） 时钟约束
DRTSC通过时钟x1描述道岔转换过程中信号保持时间1 s和转换超时判断时间13 s，包括时钟约束x1>1，x1<13和x1
$\geqslant $ 13. 如$\left({s}_{14}, \text{3}\right)\xrightarrow{{\rm{not}}\;{e}_{7},\;{x}_{\text{1}}\geqslant 13/{e}_{3}}\left({s}_{15}, \text{4}\right)$ 描述道岔转换超过13 s后进入超时状态s15.3） 风险等级
DRTSC描述了道岔转换过程中风险等级的变化情况，包括风险等级升高、降低以及保持不变，风险转移矩阵如图4所示. 其中，正数表示风险等级升高，负数表示风险等级降低，0表示风险等级不变，“−”表示状态之间无转移关系.
如R78 = 1描述了迁移
$\left( {{s_{\text{7}}},{\text{3}}} \right)\xrightarrow{{{\rm{not}}{\kern 1pt} {\kern 1pt} {e_{\text{3}}}{\kern 1pt} {\kern 1pt} {\rm{and}}{\kern 1pt} {\kern 1pt} {\rm{not}}{\kern 1pt} {\kern 1pt} {e_{10}}{\kern 1pt} {\kern 1pt} /{\kern 1pt} {\kern 1pt} {e_{\text{1}}}}}$ $ \left( {{s_{\text{8}}},{\text{4}}} \right) $ ，表示道岔由SIL 3的状态s7转移到SIL 4的s8，风险等级升高1个等级，软件测试时需要重点对此类变迁进行测试.5.3.2 与TA对比
为检测RTSC模型的描述能力，选取铁路信号系统软件测试中广泛采用的TA建模方法进行对比. 道岔定位选排的TA模型如图5所示，记为DTA. DTA由4个TA构件组成TA网络实现道岔定位选排的测试需求，其中，为描述DRTSC的否定信号，如未接收到信号e5，在DTA中增加en5. DTA和DRTSC均包括1个时钟，状态、变迁和信号数量方面的比较如表2所示，DRTSC比DTA节省62%的变迁数，节省54%的状态数. 其中，DRTSC的状态数为除去层次结构后的SIMPLE状态的数量.
进一步，随机选取计算机联锁系统软件中的4个子功能，分别采用RTSC和TA进行建模，RTSC模型和TA模型在构件、状态和变迁3个方面的对比如图6所示.
通过表2和图6，相对TA建模方法，RTSC模型具有更少的构件数、状态数、变迁数和信号数. 且随着功能逻辑复杂度的增加，RTSC节省的状态数和变迁数越多，RTSC模型更为简洁、清晰，表达能力更强，可读性更高，适合描述铁路信号系统软件复杂的功能逻辑.
RTSC继承了FSM可视化、直观和简介的特点，其状态转移、层次性和并发性能够描述软件的功能逻辑，时钟约束能够描述软件的时间关系，风险等级能够描述软件不同功能之间的安全性关系. 具体的，RTSC与铁路信号系统软件测试建模需求之间的关系如表3所示.
表 2 DRTSC与DTA 对比Table 2. Comparison between DRTSC and DTA模型 构件/个 状态/个 变迁/条 信号/个 DRTSC 1 11 19 11 DTA 4 24 49 18 表 3 RTSC与测试需求的关系Table 3. Relationship between RTSC and test requirements特性 功能性 实时性 风险 状态迁移 √ × × 层次性 √ × × 并发性 √ × × 时间约束 × √ × 风险等级 × × √ 注：“√”表示RTSC特性满足建模需求；“×”则相反. 6. 结 论
1） 铁路信号系统软件测试建模方法不能很好地描述铁路信号系统软件测试的建模需求，不能同时描述软件的功能逻辑、时钟约束和风险等级3个方面的特征.
2） 以有限状态机理论为基础，扩展出时钟约束和风险等级属性，提出风险时间状态机建模方法，采用Z语言给出了风险时间状态机的定义，阐述了风险时间状态机的格局转移规则，满足铁路信号系统软件测试在功能逻辑、时钟约束和风险等级方面的建模需求.
3） 以计算机联锁系统中的道岔子系统为例，利用风险时间状态机建立了道岔子系统软件的测试模型，并选取TA建模方法进行对比，结果表明风险时间状态机描述能力更强.
4） 文章提出的方法已在计算机联锁系统的软件测试中得到应用，实践表明该方法能够满足铁路信号系统软件测试的建模需求，为铁路信号系统软件测试提供了形式化建模的理论基础.
-
表 1 各反应器编号与运行参数
Table 1. Code of simulated landfill bioreactor
垃圾重量/kg 压实密度/（kg•m–3） 曝气频率 回灌量/% 35.0 635 1∶3 10 35.0 635 1∶3 15 35.0 635 1∶3 20 35.0 635 1∶3 30 35.0 635 1∶3 100（全回灌） 
下载: 导出CSV
-
HOORNWEG D, BHADA-TATA P, KENNEDY C. Environment:waste production must peak this century[J]. Nature, 2013, 502: 615-617 doi: 10.1038/502615a SLEZAK R, KRZYSTEK L, LEDAKOWICZ S. Degradation of municipal solid waste in simulated landfill bioreactors under aerobic conditions[J]. Waste Management, 2015, 43: 293-299 doi: 10.1016/j.wasman.2015.06.017 MORELLO L, RAGA R, CRISTINA L M, et al. The S.An.A.®;concept:semi-aerobic,anaerobic,aerated bioreactor landfill[J]. Waste Management, 2017, 67: 193-202 doi: 10.1016/j.wasman.2017.05.006 SUN Yingjie, WANG Yanan, SUN Xiaojie, et al. Production characteristics of N2O during stabilization of municipal solid waste in an intermittent aerated semi-aerobic bioreactor landfill[J]. Waste Management, 2013, 33(12): 2729-2736 doi: 10.1016/j.wasman.2013.08.013 韩智勇,刘丹,李启彬. 厌氧-准好氧联合型生物反应器填埋场产气规律的研究[J]. 环境科学,2012,33(6): 2118-2124HAN Zhiyong, LIU Dan, LI Qibin. Aerogenesis evolution of the anaerobic-semiaerobic bioreactor landfill[J]. Environmental Science, 2012, 33(6): 2118-2124 邱忠平,江海涛,王倩,等. 加速填埋场稳定化进程复合菌系的构建[J]. 中国环境科学,2012,32(3): 492-498 doi: 10.3969/j.issn.1000-6923.2012.03.017QIU Zhongping, JIANG Haitao, WANG Qian, et al. Construction of a multifunctional microbial community for accelerating the stabilization of landfill[J]. China Environmental Science, 2012, 32(3): 492-498 doi: 10.3969/j.issn.1000-6923.2012.03.017 GE Sai, LIU Lei, XUE Qiang, et al. Effects of exogenous aerobic bacteria on methane production and biodegradation of municipal solid waste in bioreactors[J]. Waste Management, 2016, 55: 93-98 doi: 10.1016/j.wasman.2015.11.024 HETTIARATCHI J P A, JAYASINGHE P A, BARTHOLAMEUZ E. M.Waste degradation and gas production with enzymatic enhancement in anaerobic and aerobic landfill bioreactors short communication[J]. Bioresource Technology, 2014, 159: 433-436 doi: 10.1016/j.biortech.2014.03.026 邱忠平. 生物强化技术加速填埋场稳定化进程实验研究[J]. 西南交通大学学报,2012,47(3): 533-537 doi: 10.3969/j.issn.0258-2724.2012.03.028QIU Zhongping. Experimental study on acceleration of landfill stabilization process using bioaugmentation technology[J]. Journal of Southwest Jiaotong University, 2012, 47(3): 533-537 doi: 10.3969/j.issn.0258-2724.2012.03.028 SUNA E A, ONAY T, YENIGUN O. Comparison of aerobic and anaerobic degradation of municipal solid waste in bioreactor landfills[J]. Bioresource Technology, 2008, 99(13): 5418-5426 doi: 10.1016/j.biortech.2007.11.008 HECHAM O, SOHRAB R. The mathematical model of the conversion of a landfill operation from anaerobic to aerobic[J]. Applied Mathematical Modelling, 2017, 50: 53-67 doi: 10.1016/j.apm.2017.05.020 ÖNCÜ G,REISER M,KRANERT M. Aerobic in situ stabilization of landfill konstanz dorfweiher:leachate quality after 1 year of operation[J]. Waste Management,2012,32(12): 2374-2384 SANG N N, SODA S, INOUE D, et al. Effects of intermittent and continuous aeration on accelerative stabilization and microbial population dynamics in landfill bioreactors[J]. Journal of Bioscience and Bioengineering, 2009, 108(4): 336-343 doi: 10.1016/j.jbiosc.2009.04.019 CLÉMENT R, OXARANGO L, DESCLOITRES M. Contribution of 3-D time-lapse ERT to the study of leachate recirculation in a landfill[J]. Waste Management, 2011, 31(3): 457-467 doi: 10.1016/j.wasman.2010.09.005 WHITE J K, BEAVEN R P, POWRIE W, et al. Leachate recirculation in a landfill:some insights obtained from the development of a simple 1-D model[J]. Waste Management, 2011, 31(6): 1210-1221 doi: 10.1016/j.wasman.2010.10.022 邱忠平, 刘洋, 刘源月, 等. 好氧生物反应器填埋场快速稳定的关键技术: 中国, 201510111362.2[P], 2015-03-13 水与废水监测分析方法编委会. 水与废水监测分析方法[M]. 4版. 北京: 中国环境科学出版社, 2003: 38, 210-213 周效志,桑树勋,程云环,等. 城市生活垃圾可生物降解有机质成分的测定[J]. 环境监测管理与技术,2007,19(2): 30-34 doi: 10.3969/j.issn.1006-2009.2007.02.009ZHOU Xiaozhi, SANG Shuxun, CHENG Yunhuan, et al. The method for determination of biodegradable organic substances in domestic garbage[J]. The Administration and Technique of Environmental Monitoring, 2007, 19(2): 30-34 doi: 10.3969/j.issn.1006-2009.2007.02.009 唐建,邱忠平,海维燕,等. 城市生活垃圾中纤维素含量测定方法优化[J]. 环境工程学报,2011,5(11): 2615-2618TANG Jian, QIU Zhongping, HAI Weiyan, et al. Optimization of determination conditions for cellulose content of municipal solid waste[J]. Chinese Journal of Environmental Engineering, 2011, 5(11): 2615-2618 蒋建国,张唱,黄云峰,等. 垃圾填埋场稳定化评价参数的中试实验研究[J]. 中国环境科学,2008(1): 58-62 doi: 10.3321/j.issn:1000-6923.2008.01.013JIANG Jianguo, ZHANG Chang, HUANG Yunfeng, et al. Pilot experiment on evaluation parameters of landfill stabilization process[J]. China Environmental Science, 2008(1): 58-62 doi: 10.3321/j.issn:1000-6923.2008.01.013 邱忠平,刘洋,连红民,等. 温度对好氧生物反应器填埋场稳定化进程的影响[J]. 西南交通大学学报,2013,48(3): 574-579QIU Zhongping, LIAN Hongmin, LIU Yang, et al. Effect of temperature on stabilizing process of aerobic bioreactor landfill[J]. Journal of Southwest Jiaotong University, 2013, 48(3): 574-579 -
下载:
下载:
点击查看大图
计量
- 文章访问数: 503
- HTML全文浏览量: 235
- PDF下载量: 6
- 被引次数: 0