Research article Special Issues

A framework for multi-perspective process mining into a BPMN process model

  • Process mining is mainly focused on process discovery from control perspective. It is further applied to mine the other perspectives such as time, data, and resources by replaying the events in event logs over the initial process model. When process mining is extended far beyond discovering the control-flow models to capture additional perspectives; roles, bottlenecks, amounts of time passed, guards, and routing probabilities in the process can be identified. This is a such extensions are considered under the topic of multi-perspective process mining, which makes the discovered process model more understandable. In this study, a framework for applying multi-perspective process mining and creating a Business Process Modelling Notation (BPMN) process model as the output is introduced. The framework, which uses a recently developed application programming interface (API) for storing the BPMN Data Model which keeps what is produced from each perspective as an asset into a private blockchain in a secure and immutable way, has been developed as a plugin to the ProM tool. In doing so, it integrates a number of techniques for multi-perspective process mining in literature, for the perspectives of control-flow, data, and resource; and represents a holistic process model by combining the outputs of these in the BPMN Data Model. In this article, we explain technical details of the framework and also demonstrate its usage over a case in medical domain.

    Citation: Merve Nur TİFTİK, Tugba GURGEN ERDOGAN, Ayça KOLUKISA TARHAN. A framework for multi-perspective process mining into a BPMN process model[J]. Mathematical Biosciences and Engineering, 2022, 19(11): 11800-11820. doi: 10.3934/mbe.2022550

    Related Papers:

    [1] Zhongcai Zhu, Xiaomei Feng, Xue He, Hongpeng Guo . Mirrored dynamics of a wild mosquito population suppression model with Ricker-type survival probability and time delay. Mathematical Biosciences and Engineering, 2024, 21(2): 1884-1898. doi: 10.3934/mbe.2024083
    [2] Chen Liang, Hai-Feng Huo, Hong Xiang . Modelling mosquito population suppression based on competition system with strong and weak Allee effect. Mathematical Biosciences and Engineering, 2024, 21(4): 5227-5249. doi: 10.3934/mbe.2024231
    [3] Rajivganthi Chinnathambi, Fathalla A. Rihan . Analysis and control of Aedes Aegypti mosquitoes using sterile-insect techniques with Wolbachia. Mathematical Biosciences and Engineering, 2022, 19(11): 11154-11171. doi: 10.3934/mbe.2022520
    [4] Liming Cai, Shangbing Ai, Guihong Fan . Dynamics of delayed mosquitoes populations models with two different strategies of releasing sterile mosquitoes. Mathematical Biosciences and Engineering, 2018, 15(5): 1181-1202. doi: 10.3934/mbe.2018054
    [5] Bo Zheng, Lihong Chen, Qiwen Sun . Analyzing the control of dengue by releasing Wolbachia-infected male mosquitoes through a delay differential equation model. Mathematical Biosciences and Engineering, 2019, 16(5): 5531-5550. doi: 10.3934/mbe.2019275
    [6] Mugen Huang, Moxun Tang, Jianshe Yu, Bo Zheng . The impact of mating competitiveness and incomplete cytoplasmic incompatibility on Wolbachia-driven mosquito population suppressio. Mathematical Biosciences and Engineering, 2019, 16(5): 4741-4757. doi: 10.3934/mbe.2019238
    [7] Yang Li, Jia Li . Stage-structured discrete-time models for interacting wild and sterile mosquitoes with beverton-holt survivability. Mathematical Biosciences and Engineering, 2019, 16(2): 572-602. doi: 10.3934/mbe.2019028
    [8] Mingzhan Huang, Shouzong Liu, Xinyu Song . Study of the sterile insect release technique for a two-sex mosquito population model. Mathematical Biosciences and Engineering, 2021, 18(2): 1314-1339. doi: 10.3934/mbe.2021069
    [9] Hui Wan, Huaiping Zhu . A new model with delay for mosquito population dynamics. Mathematical Biosciences and Engineering, 2014, 11(6): 1395-1410. doi: 10.3934/mbe.2014.11.1395
    [10] Xiaoli Wang, Junping Shi, Guohong Zhang . Bifurcation analysis of a wild and sterile mosquito model. Mathematical Biosciences and Engineering, 2019, 16(5): 3215-3234. doi: 10.3934/mbe.2019160
  • Process mining is mainly focused on process discovery from control perspective. It is further applied to mine the other perspectives such as time, data, and resources by replaying the events in event logs over the initial process model. When process mining is extended far beyond discovering the control-flow models to capture additional perspectives; roles, bottlenecks, amounts of time passed, guards, and routing probabilities in the process can be identified. This is a such extensions are considered under the topic of multi-perspective process mining, which makes the discovered process model more understandable. In this study, a framework for applying multi-perspective process mining and creating a Business Process Modelling Notation (BPMN) process model as the output is introduced. The framework, which uses a recently developed application programming interface (API) for storing the BPMN Data Model which keeps what is produced from each perspective as an asset into a private blockchain in a secure and immutable way, has been developed as a plugin to the ProM tool. In doing so, it integrates a number of techniques for multi-perspective process mining in literature, for the perspectives of control-flow, data, and resource; and represents a holistic process model by combining the outputs of these in the BPMN Data Model. In this article, we explain technical details of the framework and also demonstrate its usage over a case in medical domain.



    Mosquito-borne diseases such as dengue, yellow fever, and Zika are threatening more than half the world's population. Due to the lack of vaccines, the primary traditional method in control of these mosquito-borne diseases to suppress the mosquito population density by spraying insecticides. However, this method failed to achieve a sustainable effect on keeping mosquito population density below the critical level of epidemic risk. Even worse, heavy applications of pesticides have led to insecticide resistance and environmental pollution. In recent years, releasing sterile mosquitoes has provided an effective and biologically safe control method for eliminating or reducing mosquito populations and thus to control mosquito-borne diseases. In such a method, male mosquitoes are first sterilized using radiological or chemical techniques, and then released into the field to sterilize wild females. A wild female mosquito that mates with a sterile male will either be non-reproductive, or lay eggs that are not hatchable. Repeating releases of sterile mosquitoes could eventually wipe out wild mosquitoes, or, more realistically, suppress the wild mosquito population [1,2].

    Various mathematical models have been developed to study the interactive dynamics of wild and sterile or Wolbachia-infected mosquitoes, including ordinary differential equation models [3,4,5,6,7,8,9], delay differential equation models [10,11,12,13,14], partial differential equation models [15,16,17,18], and stochastic dynamical equation models [19,20], to cite only a few. Recently, Li in [7] formulated a simple model with constant release rate of sterile mosquitoes

    {dSν(t)dt=aSν(t)Sν(t)+Sg(t)(1ξνSν(t))Sν(t)μ1Sν(t),dSg(t)dt=bμ1Sg(t), (1.1)

    where Sν(t) and Sg(t) are the numbers of wild and sterile mosquitoes at time t respectively, a is the number of offspring produced per individual female adult per unit of time, b is the constant release rate of sterile mosquitoes, ξν is the carrying capacity parameter such that 1ξνSν describes the effect of density dependence, and μ1 is the death rate of wild or sterile mosquitoes. In (1.1), it was assumed implicitly that the mating of sterile male mosquitoes with wild females has an instant impact on the reproduction such that the model is based on ordinary differential equations without time-delay. It does not incorporate the development or maturation of mosquitoes that undergo four distinct life stages.

    Including the maturation process, we let τ be the average waiting duration from the eggs to the eclosion of adults in the next generation. We then extend model (1.1) to the following model with time delay

    {dSν(t)dt=aeμ0τSν(tτ)Sν(tτ)+Sg(tτ)(1ξνSν(tτ))Sν(tτ)μ1Sν(t),dSg(t)dt=bμ2Sg(t). (1.2)

    In this new model eμ0τ is the survival rate of the immature mosquitoes that were born at time tτ and are still alive at the time t, and μ1,μ2 are the death rates of wild and sterile mosquitoes, respectively.

    The initial condition for (1.2) is given as

    ϕ(t)=(ϕν(t),ϕg(t))C1([t0τ,t0],(0,1ξν)×(0,)), (1.3)

    where t0τ. For convenience, we write

    M=maxt[t0τ,t0]ϕg(t),M0=max{b/μ2,M},m=mint[t0τ,t0]ϕg(t),m0=min{b/μ2,m}. (1.4)

    It is well-known that the oscillation phenomena may appear when time delays are included in differential equations, as shown, for example, in recent studies of mosquito population models with time delays [10,12]. We aim to offer an accurate description of the oscillation phenomena of (1.2), and study how the oscillatory properties change with respect to the release of sterile mosquitoes. Specifically to our model (1.2)-(1.3), we define the oscillations as follows.

    Definition 1.1. Let (Sν(t),Sg(t)) be a non-constant positive solution of (1.2)-(1.3). We say that (Sν(t),Sg(t)) is non-oscillatory about its equilibrium (¯Sν,¯Sg) if Sν(t)¯Sν and Sg(t)¯Sg are positive or negative eventually. Otherwise, it is said to be oscillatory about (¯Sν,¯Sg). We say that (1.2) is oscillatory about (¯Sν,¯Sg) if every non-constant solution of (1.2) is oscillatory. Otherwise, (1.2) is said to be non-oscillatory about (¯Sν,¯Sg).

    In this paper, we study the oscillatory properties of (1.2). The paper is organized as follows. In Section 2, we establish useful lemmas and find that, (1.2) is oscillatory about the unique non-zero equilibrium when no sterile mosquitoes are released. In Section 3, we determine an oscillation threshold, denoted by ˆb, for the constant release rate b of sterile mosquitoes, such that all non-trivial positive solutions oscillate when b<ˆb, and the oscillation disappears when b>ˆb. In Section 4, numerical simulations are provided to demonstrate our new findings. Concluding remarks are finally given in Section 5.

    In this section, we first establish useful lemmas that help us to prove the main results of this paper. In addition, we show that every non-trivial solution of (1.2) oscillates with respect to its unique non-zero equilibrium when no sterile mosquitoes are released. We first determine the monotonicity of the following birth-progression function.

    Lemma 2.1. Consider the birth-progression function

    g(x,z)=aeμ0τ(1ξνx)x2x+z,0<x<1/ξν,z>0,

    and let

    C(z)=13ξνz+(13ξνz)2+16ξνz4ξν.

    Then

    gx(x,m0){>0,0<x<C(m0),=0,x=C(m0),<0,C(m0)<x<1ξν.

    Proof. We omit the detail of the proof since the result can be obtained by directly taking the derivative of g(x,z) with respect to x.

    We next discuss the boundedness and positivity of solutions of (1.2).

    Lemma 2.2. Suppose μ1>ξνg(C(m0),m0). Then system (1.2)-(1.3) has a unique solution (Sν,Sg) which is bounded and positive for all t[t0,).

    Proof. The existence and uniqueness of the solution of (1.2) follows from the standard results in the theory of delay differential equations [21]. From the second equation of (1.2), it is easy to see that 0<m0Sg(t)M0 for all tt0, where m0 and M0 are defined in (1.4). It remains to verify the positivity and boundedness of Sν in (1.2).

    We first confirm that Sν(t)<1/ξν holds for all tt0 when μ1>ξνg(C(m0),m0). Otherwise, there exists t1>t0 such that Sν(t1)=1/ξν, and Sν(t)<1/ξν for all t[t0τ,t1). Hence Sν(t1)0. However, from (1.4) and the first equation of (1.2), we have

    Sν(t1)μ1Sν(t1)+g(Sν(t1τ),m0)1ξν(μ1ξνg(C(m0),m0))<0,

    which is a contradiction to Sν(t1)0.

    Next, we prove the positivity of Sν(t). If not, there exists t2>t0 such that Sν(t2)=0 and 0<Sν(t)<1/ξν for all t[t0τ,t2). Hence Sν(t2)0, which contradicts

    Sν(t2)=aeμ0τ(1ξνSν(t2τ))S2ν(t2τ)Sν(t2τ))+Sg(t2τ))>0.

    The proof is complete.

    Remark 1. We assume that the condition μ1>ξνg(C(m0),m0) holds throughout the rest of this paper.

    Define the intrinsic growth rate of wild mosquito population by

    r0=aeμ0τμ1.

    The following lemma clarifies the existence of equilibria of (1.2).

    Lemma 2.3. Define

    b:=μ2(r01)24ξνr0. (2.1)

    If 0<b<b, then (1.2) has three nonnegative equilibria: N0:=(0,bμ2) and N=(Sν,bμ2), where

    Sν:=(r01)μ2(r01)2μ224bξνr0μ22ξνr0μ2;

    When b=b, (1.2) has two nonnegative equilibria: N0 and N:=(r012r0ξν,bμ2); When b>b, N0 is the only nonnegative equilibrium of (1.2).

    We next investigate the oscillation of (1.2) with respect to the unique non-zero equilibrium when b=0. In this case, (1.2) becomes

    {dSν(t)dt=aeμ0τSν(tτ)Sν(tτ)+Sg(tτ)(1ξνSν(tτ))Sν(tτ)μ1Sν(t),dSg(t)dt=μ2Sg(t). (2.2)

    It is clear that limtSg(t)=0. Therefore, the oscillation of (2.2) is exactly the same as that of the next equation

    dSν(t)dt=aeμ0τSν(tτ)(1ξνSν(tτ))μ1Sν(t), (2.3)

    which has a unique positive equilibrium S(1)ν=(r01)/r0ξν.

    We now study the oscillation of (2.3) with respect to S(1)ν when r0>2. Let y(t)=Sν(t)S(1)ν. Then instead of considering the oscillation of Sν(t) with respect to S(1)ν in (2.3), we consider the oscillation of y(t) about y(t)=0 for the following equation

    y(t)=μ1y(t)qy(tτ)μ1r0ξνy2(tτ), (2.4)

    where q=μ1(r02)>0.

    To proceed, we first establish the following three lemmas to reach the conclusion of Theorem 1.

    Lemma 2.4. If r0>2, then every non-oscillatory solution of (2.4) about y(t)=0 converges to zero as t.

    Proof. Let y(t) be an arbitrary non-oscillatory solution of (2.4) about y(t)=0. Then it is positive or negative eventually. It suffices to prove that

    limty(t)=0. (2.5)

    (ⅰ) Suppose that y(t) is eventually positive. Then there exists ¯t1>t0 large enough such that y(t)>0 for all t>¯t1+τ. In the meantime, based on (2.4), we have

    y(t)=μ1y(t)qy(tτ)μ1r0ξνy2(tτ)<0.

    Thus limty(t)=L holds and we claim that L is non-negative. If L>0, then

    limty(t)=μ1LqLμ1r0ξνL2:=ϱ<0.

    Thus we have y(t)ϱ<0 for sufficiently large t, and hence limty(t)=, which is a contradiction. Therefore, limty(t)=L=0.

    (ⅱ) Suppose that y(t) is eventually negative. To prove (2.5), we let

    ¯y=lim supty(t),y_=lim infty(t).

    It is easy to see that (2.5) is true if and only if y_=0. Otherwise, y_<0. There are only two cases to consider: y_=¯y, and y_<¯y.

    If y_=¯y<0, then limny(t)=y_=¯y exists. Taking limits in (2.4) on both sides gives

    limty(t)=y_(μ1+q+μ1r0ξνy_).

    Notice that y(t) is bounded, then we derive μ1+q+μ1r0ξνy_=0 and y_=S(1)ν. Hence there exists a sequence {tn} large enough such that y(tn)0, y(tn)S(1)ν, and y(tn)=minttny(t). Then it follows from (2.4) that

    0y(tn)=μ1y(tn)qy(tnτ)μ1r0ξνy2(tnτ),

    which leads to

    μ1y(tn)y(tnτ)[q+μ1r0ξνy(tnτ)],

    and

    μ1y(tnτ)y(tn)[q+μ1r0ξνy(tnτ)][q+μ1r0ξνy(tnτ)],

    Solving the inequality above, we obtain y(tnτ)S(1)ν, which is a contradiction to y(tnτ)>S(1)ν.

    If y_<¯y<0, then there exists a sequence {sn} large enough such that y(sn)=0, y(sn)y_. From (2.4), we have

    μ1y(sn)=qy(snτ)μ1r0ξνy2(snτ).

    Notice that y(snτ) is bounded in the interval (S(1)ν,0). Then there exists a convergent subsequence, denoted by {sn} again, such that

    μ1y_=qy1μ1r0ξνy21,

    where y_y1:=limny(snτ)¯y. Hence

    μ1=y1y_[q+μ1r0ξνy1][q+μ1r0ξνy1],

    and y1S(1)ν. Meanwhile, y1y_S(1)ν, thus y1=y_=S(1)ν.

    Take t>t0 large enough such that y(t)=0, and y(t)=mintty(t). Then it follows from (2.4) that

    μ1y(t)=qy(tτ)μ1r0ξνy2(tτ),

    which yields

    μ1=y(tτ)y(t)[q+μ1r0ξνy1][q+μ1r0ξνy(tτ)],

    and y(tτ)S(1)ν, a contradiction to y(tτ)>S(1)ν. The proof is complete.

    To show the oscillatory behavior of solutions of (2.4) about y(t)=0, we linearize (2.4) at y(t)=0 which leads to

    y(t)=μ1y(t)qy(tτ). (2.6)

    Based on the results in [22], we immediately have the following necessary and sufficient condition for the oscillation of solutions of (2.6) with respect to y(t)=0.

    Lemma 2.5. [22] Every non-trivial solution of (2.6) is oscillatory about y(t)=0 if and only if

    qτeμ1τ>1e, (2.7)

    where q=μ1(r02)>0.

    We then show the equivalence of the oscillations between (2.4) and (2.6) about y=0.

    Lemma 2.6. Equation (2.4) is oscillatory about y=0 if and only if (2.6) is oscillatory about y=0.

    Proof. Assume that (2.4) is oscillatory about y=0. We confirm that (2.6) is also oscillatory about y=0. Otherwise, there exists a solution y(t) of (2.6) that is non-oscillatory about y(t)=0. Without loss of generality, suppose that y(t) is eventually negative. Then there exist δ>0 and t1>t0 such that δ<y(t)<0, for all t>t1. Let Y(t) be a solution of (2.4) with the same initial value condition as that of y(t) in (2.6). Notice that

    μ1y(t)qy(tτ)μ1ξνr0y2(tτ)μ1y(t)qy(tτ).

    By the comparison principle, we have Y(t)y(t)<0, for t>t1, which is a contradiction to the assumption that (2.4) is oscillatory about y(t)=0.

    On the other hand, we show that (2.4) oscillates about y=0 if (2.6) oscillates about y=0. If not, then there exists a solution y1(t) of (2.4), which is either eventually positive or eventually negative. There are two cases to consider.

    (ⅰ) If y1(t) is eventually positive, then there exists t2>t0 such that y1(t)>0, for all t>t2. Let Y1(t) be the solution of (2.6) with the same initial condition as that of y1(t) in (2.3). Again, by the comparison principle, we have Y1(t)y1(t)>0, for all t>t2, a contradiction to the oscillation of (2.6) about y(t)=0.

    (ⅱ) If y1(t) is eventually negative, then there exists t3>t0 such that y1(tτ)<0, y1(t)<0, for all t>t3. From Lemma 4, we have limty1(t)=0. Since (2.7) holds, there exists a positive constant ε0<1 such that

    ε0qτeμ1τ>1e,

    which suggests that the equation

    y(t)+μ1y(t)+ε0qy(tτ)=0

    is oscillatory about y(t)=0. Moreover, since

    limtq+μ1ξνr0y(tτ)q=1,

    there exists t4>t3 such that

    q+μ1ξνr0y(tτ)>ε0q,for allt>t4.

    Then

    (q+μ1ξνr0y(tτ))(y(tτ))>ε0q(y(tτ)),for allt>t4,

    and thus

    μ1y(t)qy(tτ)μ1ξνr0y2(tτ)>μ1y(t)ε0qy(tτ).

    Let ¯Y(t) be a solution of the equation y(t)+μ1y(t)+ε0qy(tτ)=0 with the same initial condition as the solution y1(t) of (2.4). Then ¯Y(t)<y1(t)<0 for all t>t4, which is a contradiction to the oscillation of the equation

    y(t)+μ1y(t)+ε0qy(tτ)=0

    about y(t)=0. Therefore, equation (2.4) is oscillatory about y=0 if and only if (2.6) is oscillatory about y=0. The proof is complete.

    We note that the oscillations of Sν about S(1)ν between system (2.2) and (2.3) are exactly the same. Based on Lemmas 2.3, 2.4, and 2.5 above, a sufficient and necessary condition for the oscillations of solutions of equation (2.2) can be summarized as follows.

    Theorem 1. If r0>2, then every non-trivial positive solution of (2.2) is oscillatory about S(1)ν if and only if (2.7) holds.

    Remark 2. Theorem 1 provides a sufficient and necessary condition for the oscillations of non-trivial positive solutions about the unique non-zero equilibrium in equation (1.2) when no sterile mosquitoes are released. In fact, from the estimation of parameters of (1.2) in Table 1, we find that (2.7) is always true and so (2.2) is oscillatory about S(1)ν, which is consistent with the oscillation phenomenon of the annual abundance of wild mosquitoes in the wild. These results are important for the discussion of the main results in Section 3.

    Table 1.  Parameter values adapted to Aedes albopictus population suppression in subtropical monsoon climate, especially in Guangzhou.
    Para. Definitions Ranges References
    τ The average waiting duration from eggs to the eclosion [16, 66] [25]
    of adults in the next generation (day)
    a The number of offsprings produced per individual, [3.15, 12.81] [26,28]
    per unit of time
    eμ0τ The survival rate of the immature mosquitoes (day1) 0.05 [27]
    ξν The carrying capacity parameter of wild mosquitoes 0.0025 Given
    μ1 The death rate of wild mosquitoes (day1) [0.0231, 0.0693] [26]
    μ2 The death rate of sterile mosquitoes (day1) 1/7 [23,24]

     | Show Table
    DownLoad: CSV

    Given that (1.2) is oscillatory about the unique non-zero equilibrium in the absence of sterile mosquitoes, it becomes interesting to ask how the oscillatory property of (1.2) with the increase of the releases of sterile mosquitoes. Interestingly, we find an oscillation threshold, denoted by ˆb, for the release rate of the sterile mosquitoes, and show that this oscillation phenomenon will be maintained when the release rate b<ˆb, whereas it will disappear when b>ˆb. The result of oscillation threshold is described as follows.

    Theorem 2. Let ˆb be the unique solution of

    q(b)τeμ1τ=1e,

    where

    q(b)=μ1μ2S+ν(3r0ξνS+ν+12r0)μ2S+ν+b>0.

    Then ˆb is the oscillation threshold of the release rate b, below which the oscillation phenomenon maintain, and above which it will disappear.

    Proof. By using the similar argument as that in Theorem 1, we conclude that every non-trivial positive solution of (1.2) is oscillatory about the non-zero equilibrium N+ if and only if

    q(b)τeμ1τ>1e. (3.1)

    The proof is nearly the same as that of Theorem 1, so we omit it.

    Next, we claim that q(b) decreases with the increase of b. It follows from the definition of S+ν in Lemma 2.3 that S+ν(b) is monotonously decreasing. Hence 3r0ξνS+ν2r0+1 is positive and monotonously decreasing with respect to b. Moreover, By taking derivative, we obtain

    ddb(S+νμ2S+ν+b)<0,

    thus S+ν/(μ2S+ν+b) is monotonously decreasing. Consequently, q(b) is monotonously decreasing.

    Notice that (1.2) is oscillatory about N+ when b=0. According to Theorem 1, we have q(0)τeμ1τ>1/e. Meanwhile, when b=b, we have

    q(b)=μ1μ2S+ν(3r0ξνS+ν+12r0)μ2S+ν+b=(12r012)μ1μ2S+νμ2S+ν+b<0.

    Thus, there exists unique ˆb(0,b) such that q(ˆb)τeμ1τ=1/e. Notice that (1.2) is oscillatory about N+ if and only if (3.1) holds, we conclude that system (1.2) is oscillatory about N+ when 0b<ˆb, and non-oscillatory when ˆbbb. The proof is complete.

    In this section, numerical simulations are given to validate the main theoretical results in Section 3.

    In Table 1 below, we list important parameter values for Aedes albopictus and sterile mosquitoes, most of which are taken from earlier experimental data in the literature [23,24,25,26,27,28,29]. Since the parameter values are affected by many factors such as temperature and rainfall, we can not determine their exact values but take reasonable estimations of the ranges. Using a similar method for the parameter estimations as in [12], we estimate τ[16,66]. By using the measured half-lives, we have μ10.0277 from [29], and μ1[0.0231,0.0693] from [26], hence the life span of adult mosquitoes is T1=1/μ136.10, and T1[14.43,43.29]. Furthermore, it takes 5 days on average for a female mosquito to lay eggs after the eclosion, and the total number of eggs laid by each female is 120.76 on average [26]. Consequently, the egg-laying phase of a female mosquito is T2[9.43,38.29], and a[3.15,12.81]. From [23,24], we find that sterile mosquitoes will die in about 7 days after releasing to the wild field. Thus we estimate that the death rate of sterile mosquitoes is about 1/7. Besides, according to [27], we obtain that the survival probability of the average maturation time of wild mosquitoes is about 0.05. We note that ξν is usually proportional to the area size. From simulations, it shows that different values of ξν produce similar dynamics, which allows us to take ξν=0.0025 as a representative value.

    To make the oscillation behavior be clearer, we employ simulation examples with different release rates to testify the result of oscillation threshold given in Theorem 2.

    Suppose that the model parameters of (1.2) are given as

    a=4,τ=36,eμ0τ=0.05,μ1=0.05,μ2=1/7,ξν=0.0025. (4.1)

    All the parameter values in (4.1) are within the estimation ranges given in Table 1. By substituting these parameter values into (2.1), we derive b32.143. Notice that the oscillation threshold ˆb is the unique solution of the equation q(b)τeμ1τ=1/e, we have ˆb21.7754.

    We select three different release rates b=0, b=3 and b=6 from the interval [0,ˆb). For these three cases, the number of wild mosquitoes Sν(t) oscillates around S+ν, as shown in Figure 1. This testifies the first part of Theorem 2, that is, the non-trivial positive solutions of (1.2) are oscillatory about N+ when b<ˆb. When b=0, there is no sterile mosquitoes released into the wild, the number of mosquitoes oscillates with a larger amplitude, as shown in Figure 1 (A). As the release rate increases, the amplitude is reduced. When b=6, the number of mosquitoes oscillates in a narrow range around S+ν, as shown in Figure 1 (C).

    Figure 1.  The oscillatory behavior of (1.2) with respect to N+ when the release rates b=0,3,6 are smaller than the oscillation threshold ˆb21.7754.

    We increase the release rate of sterile mosquitoes further, and let ˆb<b<b. For instance, we take b=24, b=27 and b=30 separately from the interval [ˆb,b]. As shown in Figure 2, we find that the number Sν(t) is less than S+ν, and then gradually increases and approaches to S+ν. This phenomenon testifies the second part of Theorem 2, that is, the solutions (Sν,Sg) are non-oscillatory about N+.

    Figure 2.  The release rates b=24, b=27 and b=30, respectively, are larger than the oscillation threshold ˆb21.7754. The corresponding solutions (Sν,Sg) of (1.2) are non-oscillatory about N+.

    Figure 2 also shows that the wild mosquito population is first suppressed to a low level by releasing sterile mosquitoes with higher rates, then gradually tends to a stable level S+ν. Figure 1 and 2 both show that the stable level of wild mosquitoes is reduced when the release rate increases from 0 to b.

    To control the life-threatening mosquito-borne diseases such as dengue, yellow fever, and Zika, one of the biologically safe methods is to release male sterile mosquitoes into the field to suppress the wild mosquito population. We, in this paper, consider the oscillatory properties of a delayed mosquito population suppression model with a constant release rate of sterile mosquitoes. We find that every non-trivial positive solution of the model oscillates with respect to its unique non-zero equilibrium when no sterile mosquitoes are released, i.e., b=0. We then study the oscillatory behavior of the model with the releases of sterile mosquitoes, and establish an oscillation threshold, denoted by ˆb, for the constant release rate b of sterile mosquitoes such that the oscillation of solutions maintains when b<ˆb, whereas the oscillation disappears when b>ˆb. Furthermore, through numerical examples, we show that the amplitudes of the oscillation become smaller and smaller and then the oscillation disappears as release rates of sterile mosquitoes increase. Oscillatory phenomenon is common for solutions of delay differential equations. However, to the best of our knowledge, an establishment of the oscillation threshold which determines the existence or disappearance of solution oscillations is one of the first in the field of mosquito population dynamics.

    This work was supported by the National Natural Science Foundation of China (No. 11631005), the Program for Changjiang Scholars and Innovative Research Team in University (No. IRT-16R16), the Guangzhou Postdoctoral International Training Program Funding Project, the Science and Technology Program of Guangzhou (No. 201707010337) and the Foundation of Yunnan Educational Committee (No. 2017ZDX027). The authors thank Professor Jianshe Yu for his help in improving the original Lemma 2.4.

    The authors have declared that no competing interests exist.



    [1] M. Dumas, M. La Rosa, J. Mendling, H. A. Reijers, Fundamentals of business process management, Springer, (2013). https://doi.org/10.1007/978-3-642-33143-5
    [2] A. J. M. M Weijters, W. M. van Der Aalst, A. A. De Medeiros, Process mining with the heuristics miner-algorithm, Technische Universiteit Eindhoven Tech. Rep. WP, 166 (2006), 1–34.
    [3] M. De Leoni, W. M. van der Aalst, M. Dees, A general process mining framework for correlating, predicting and clustering dynamic behavior based on event logs, Inform. Syst., 56 (2016), 235–257. https://doi.org/10.1016/j.is.2015.07.003 doi: 10.1016/j.is.2015.07.003
    [4] F. Mannhardt, Multi-perspective process mining, BPM (Dissertation/Demos/Industry), (2018).
    [5] J. L. Peterson, Petri net theory and the modeling of systems, Prentice Hall PTR, (1981).
    [6] R. M. Dijkman, M. Dumas, C. Ouyan, Semantics and analysis of business process models in BPMN, Inform. Software Technol., 50 (2008), 1281–1294. https://doi.org/10.1016/j.infsof.2008.02.006 doi: 10.1016/j.infsof.2008.02.006
    [7] J. M. Colom, J. Desel, Application and theory of petri nets and concurrency, Springer, Italy, 2013. https://doi.org/10.1007/978-3-319-19488-2
    [8] W. M. Van der Aalst, B. F. van Dongen, C. W. Günther, A. Rozinat, E. Verbeek, T. Weijters, ProM: The process mining toolkit, BPM (Demos), 489 (2009).
    [9] B. Ekici, A BPMN data model to keep a multi-perspective process model on the blockchain, Fen Bilimleri Enstitüsü, Hacettepe University in Turkey, 2021.
    [10] B. Ekici, T. G. Erdogan, A. K. Tarhan, BPMN data model for multi-perspective process mining on blockchain, Int. J. Software Eng. Knowl. Eng., (2022), 1–29. https://doi.org/10.1142/S0218194022500115
    [11] T. G. Erdogan, A. K. Tarhan, Multi-perspective process mining for emergency process, Health Inform. J., 28 (2022), 14604582221077195. https://doi.org/10.1177/14604582221077195 doi: 10.1177/14604582221077195
    [12] W. Van Der Aalst, Process mining: Data science in action, Springer, (2016). https://doi.org/10.1007/978-3-662-49851-4_1
    [13] W. V. D. Aalst, A. Adriansyah, A. K. A. D. Medeiros, F. Arcieri, T. Baier, T. Blickle, et al., Process mining manifesto, in Springer, International conference on business process management, (2011), 169–194. https://doi.org/10.1007/978-3-642-28108-2_19
    [14] M. Rovani, F. M. Maggiand, M. De Leoni, W. M. Van Der Aalst, Declarative process mining in healthcare, Expert Syst. Appl., 42 (2015), 9236–9251. https://doi.org/10.1016/j.eswa.2015.07.040 doi: 10.1016/j.eswa.2015.07.040
    [15] W. Van der Aalst, T. Weijters, L. Maruster, Workflow mining: Discovering process models from event logs, IEEE Transact. Knowl. Data Eng., 16 (2004), 1128–1142. https://doi.org/10.1109/TKDE.2004.47 doi: 10.1109/TKDE.2004.47
    [16] M. De Leoni, W. M. van der Aalst, Data-aware process mining: discovering decisions in processes using alignments, in Proceedings of the 28th annual ACM symposium on applied computing, (2013), 1454–1461. https://doi.org/10.1145/2480362.2480633
    [17] S. J. Leemans, D. Fahland, W. M. Van Der Aalst, Process and deviation exploration with inductive visual miner, BPM (demos), 1295 (2014).
    [18] A. A. F. G. Mohamed, Process mining application considering the organizational perspective using social network analysis, University of Porto, 2016.
    [19] M. Song, W. M. van der Aalst, Supporting process mining by showing events at a glance, in Proceedings of the 17th Annual Workshop on Information Technologies and Systems (WITS), (2007), 139–145.
    [20] F. Mannhardt, M. De Leoni, H. A. Reijers, The multi-perspective process explorer, BPM (Demos), 1418 (2015), 130–134.
    [21] T. G. Erdogan, A. Tarhan, Systematic mapping of process mining studies in healthcare, IEEE Access, 6 (2018), 24543–24567. https://doi.org/10.1109/ACCESS.2018.2831244 doi: 10.1109/ACCESS.2018.2831244
    [22] A. Rozinat, C. W. Günther, R. Niks, Process mining and automated process discovery software for professionals-Fluxicon Disco, (2017).
    [23] E. Rojas, C. Fernández-Llatas, V. Traver, V. Munoz-Gama, M. Sepúlveda, V. Herskovic et al., PALIA-ER: Bringing question-driven process mining closer to the emergency room, in BPM (Demos), (2017).
    [24] Process Mining and Execution Management Software, Celonis, 2022. Available from: https://www.celonis.com/.
    [25] R. Gatta, J. Lenkowicz, M. Vallati, E. Rojas, A. Damiani, L. Sacchi, et al., pMineR: An innovative R library for performing process mining in medicine, in Conference on artificial intelligence in medicine in europe, (2017), 351–355. https://doi.org/10.1007/978-3-319-59758-4_42
    [26] G. Janssenswillen, B. Depaire, M. Swennen, M. Jans, K. Vanhoof, bupaR: Enabling reproducible business process analysis, Knowledge-Based Syst., 163 (2019), 927–930. https://doi.org/10.1016/j.knosys.2018.10.018 doi: 10.1016/j.knosys.2018.10.018
    [27] U. Celik, E. Akcetin Surec madenciligi araclari karsilastirmasi AJIT-e: Bilişim Teknolojileri Online Dergisi, 9 (2018), 97–104. https://doi.org/10.5824/1309-1581.2018.4.007.x
    [28] D. Dakic, D. Stefanovic, I. Cosic, T. Lolic, M. Medojevic, Business process mining application: A literature review, Ann. DAAAM Proceed., 29 (2018). https://doi.org/10.2507/29th.daaam.proceedings.125 doi: 10.2507/29th.daaam.proceedings.125
    [29] A. Kalenkova, A. Burattin, M. de Leoni, W. van der Aalst, A. Sperduti, Discovering high-level BPMN process models from event data, Business Process Manag. J., (2019). https://doi.org/10.1108/BPMJ-02-2018-0051
    [30] D. Bell, UML basics: An introduction to the Unified Modeling Language, Rational Edge, (2003).
    [31] J. Mendling, M. Nüttgens, EPC markup language (EPML): An XML-based interchange format for event-driven process chains (EPC), Inform. Syst. E-business Manag., 4 (2006), 245–263. https://doi.org/10.1007/s10257-005-0026-1 doi: 10.1007/s10257-005-0026-1
    [32] M. Geiger, S. Harrer, J. Lenhard, G. Wirtz, BPMN 2.0: The state of support and implementation, Future Gener. Computer Syst., 80 (2018), 250–262. https://doi.org/10.1016/j.future.2017.01.006 doi: 10.1016/j.future.2017.01.006
    [33] J. C. Buijs, B. F. V. Dongen, W.M van Der Aalst, On the role of fitness, precision, generalization and simplicity in process discovery, in OTM Confederated International Conferences" On the Move to Meaningful Internet Systems", (2012), 305–322. https://doi.org/10.1007/978-3-642-33606-5_19
    [34] J. Carmona Vargas, J. Cortadella, M. Kishinevsky, Region-based algorithms for process mining and synthesis of Petri nets, Polytechnic University of Catalonia, (2009).
    [35] R. Bergenthum, J. Desel, R. Lorenz, S. Mauser, Process mining based on regions of languages, in International Conference on Business Process Management, (2007), 375–383. https://doi.org/10.1007/978-3-540-75183-0_27
    [36] J. M. E. Van der Werf, B. F. van Dongen, C. A. Hurkens, A. Serebrenik, Process discovery using integer linear programming, in International conference on applications and theory of petri nets, (2008), 368–387 https://doi.org/10.1007/978-3-540-68746-7_24
    [37] C. W. Günther, W. M. Van Der Aalst, Fuzzy mining–adaptive process simplification based on multi-perspective metrics, in International conference on business process management, (2007), 328–343. https://doi.org/10.1007/978-3-540-75183-0_24
    [38] W. M. van der Aalst, On the representational bias in process mining, in 2011 IEEE 20th International Workshops on Enabling Technologies: Infrastructure for Collaborative Enterprises, (2011), 2–7. https://doi.org/10.1109/WETICE.2011.64
    [39] W. M. Van der Aalst, A. K. De Medeiros, A. J. Weijters, Genetic process mining, in International conference on application and theory of petri nets, (2005), 48–69. https://doi.org/10.1007/11494744_5
    [40] P. Weber, B. Bordbar, P. Tiňo, A principled approach to the analysis of process mining algorithms, in International Conference on Intelligent Data Engineering and Automated Learning, (2011), 474–481. https://doi.org/10.1007/978-3-642-23878-9_56
    [41] R. Ghawi, Process discovery using inductive miner and decomposition, arXiv preprint arXiv: 1610.07989, (2016).
    [42] A. Rozinat, W. M. van der Aalst, Decision mining in ProM, in International Conference on Business Process Management, (2006), 420–425. https://doi.org/10.1007/11841760_33
    [43] R. S. Mans, M. H. Schonenberg, M. Song, W. M. van der Aalst, P. J. Bakker, Application of process mining in healthcare–a case study in a dutch hospital, in International joint conference on biomedical engineering systems and technologies, (2008), 425–438. https://doi.org/10.1007/978-3-540-92219-3_32
    [44] M. Bozkaya, J. Gabriels, J. M. van der Werf, Process diagnostics: A method based on process mining, in International Conference on Information, Process, and Knowledge Management, 1 (2009), 22–27. https://doi.org/10.1109/eKNOW.2009.29
    [45] W. Van Der Aalst, Service mining: Using process mining to discover, check, and improve service behavior, IEEE Transact. Services Comput., 6 (2012), 525–535. https://doi.org/10.1109/TSC.2012.25 doi: 10.1109/TSC.2012.25
    [46] E. Gupta, Process mining algorithms, Int. J. Adv. Res. Sci. Eng., 3 (2014), 401–412.
    [47] F. Folino, G. Greco, A. Guzzo, L. Pontieri, Discovering multi-perspective process models: The case of loosely-structured processes, in International Conference on Enterprise Information Systems, 19 (2008). https://doi.org/10.1007/978-3-642-00670-8_10
    [48] A. Pini, R. Brown, M. T. Wynn, Process visualization techniques for multi-perspective process comparisons, in Asia-Pacific Conference on Business Process Management, (2015), 183–197. https://doi.org/10.1007/978-3-319-19509-4_14
    [49] S. Schönig, C. D. Ciccio, F. M. Maggi, J. Mendling, Discovery of multi-perspective declarative process models, International Conference on Service-Oriented Computing, 9936 (2016), 87–103. https://doi.org/10.1007/978-3-319-46295-0_6 doi: 10.1007/978-3-319-46295-0_6
    [50] S. Jablonski, M. Röglinger, S. Schönig, K. M. Wyrtki, Multi-perspective clustering of process execution traces, Enterprise Model. Inform. Syst. Archit. (EMISAJ), 14 (2019), 2–11. https://doi.org/10.18417/emisa.14.2 doi: 10.18417/emisa.14.2
    [51] R. Sikal, H. Sbai, L. Kjiri, Promoting resource discovery in business process variability, in Proceedings of the 2nd International Conference on Networking, Information Systems & Security, (2019), 1–7. https://doi.org/10.1145/3320326.3320380
    [52] C. Sturm, S. Schönig, C. Di Ciccio, Distributed multi-perspective declare discovery, in BPM (Demos), (2017).
    [53] Tijs Slaats, Flexible process notations for cross-organizational case management systems, IT University of Copenhagen, Theoretical computer Science section, 2015.
    [54] A. J. M. M. Weijters, J. T. S. Ribeiro, Flexible heuristics miner (FHM), in Flexible heuristics miner (FHM), (2011). https://doi.org/10.1109/CIDM.2011.5949453
    [55] W. V. D. Aalst, A. Adriansyah, B. V. Dongen, Causal nets: A modeling language tailored towards process discovery, in International conference on concurrency theory, 6901 (2011), 28–42. https://doi.org/10.1007/978-3-642-23217-6_3
    [56] V. M. Van der Aalst, V. Rubin, H. M. W. Verbeek, B. F. van Dongen, E. Kindler, C. W. Günther, Process mining: a two-step approach to balance between underfitting and overfitting, Software Syst. Model., 9 (2010), 87–111. https://doi.org/10.1007/s10270-008-0106-z doi: 10.1007/s10270-008-0106-z
    [57] F. Mannhardt, D. Blinde, Analyzing the trajectories of patients with sepsis using process mining, in RADAR+ EMISA@ CAiSE, (2017).
    [58] Á. Rebuge, D. R. Ferreira, Business process analysis in healthcare environments: A methodology based on process mining, Inform. Syst., 37 (2012), 99–116. https://doi.org/10.1016/j.is.2011.01.003 doi: 10.1016/j.is.2011.01.003
    [59] ProM 6.10, Accessed on 31.12.2021. Available from: http://www.promtools.org/doku.php?id=prom610.
  • This article has been cited by:

    1. Genghong Lin, Yuanxian Hui, Stability analysis in a mosquito population suppression model, 2020, 14, 1751-3758, 578, 10.1080/17513758.2020.1792565
    2. Manuel De la Sen, Asier Ibeas, Aitor J. Garrido, Stage-Dependent Structured Discrete-Time Models for Mosquito Population Evolution with Survivability: Solution Properties, Equilibrium Points, Oscillations, and Population Feedback Controls, 2019, 7, 2227-7390, 1181, 10.3390/math7121181
    3. Liping Luo, Zhenguo Luo, Yunhui Zeng, New results for oscillation of fractional partial differential equations with damping term, 2020, 0, 1937-1179, 0, 10.3934/dcdss.2020336
    4. Yuanxian Hui, Zhong Zhao, Qiuying Li, Liuyong Pang, Asymptotic stability in a mosquito population suppression model with time delay, 2023, 16, 1793-5245, 10.1142/S1793524522500929
  • Reader Comments
  • © 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3614) PDF downloads(199) Cited by(8)

Figures and Tables

Figures(10)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog