Abstract
In this article, we present a contribution to modeling, evaluation, and analysis of the inventory management systems performance and more generally stochastic discrete event systems with a batch behavior. For this contribution, we combine two models: the Supply Chain Operations Reference model, proposed by the Supply Chain Council, with the Batch Deterministic and Stochastic Petri Nets, which constitutes a very powerful dynamic modeling tool. To do this, we applied these tools on a typical model of inventory management system in order to show how the combination of these two tools can help us to model and analyze the performance of the inventory management system and to provide information on their behavior and the effects of their parameters. A resolution of the stochastic process associated with the warehouse management system will allow us to calculate the following performance indicators: average stock, average cost of stock, probability of an empty stock, and average supply and frequency. These indicators will help to monitor the activity of our stock management system, and therefore make the right decisions for the development of the organization.
Introduction
Modeling, evaluation, and performance analysis of discrete event systems including production systems, storage systems, and especially the supply chain (SC) remain a major concern for different scientific communities and especially manufacturing ones.
Historically, research in modeling and performance evaluation had always involved two closely linked issues. In addition, we can clearly identify in practice a combination of performance assessment methods to certain modeling formalisms. The first works are based on the direct analysis of stochastic processes and more generally of Markov chains. 1 –4 The work of Jackson 5 led, in the late of 50s, a new approach based on queuing networks. Other approaches have abandoned the formalism of queuing networks and adopt a new and more powerful formalism introduced in the early 60s. It is the Petri nets (PNs) 6 –8 formalism as basic form to which they added extensions ranging from simple timers constants 9 until mechanisms more sophisticated such as stochastic Petri nets (SPNs), 10 generalized stochastic Petri nets (GSPNs), 11,12 colored PNs, 13,14 and continuous PNs. 15
In this article, we introduce the Batch Deterministic and Stochastic Petri Nets (BDSPNs) as a tool for modeling and performance evaluation of SCs. This model is developed by enhancing DSPNs 16 with batch places and batch tokens.
The increasing complexity of logistics processes, the desire of a better configuration of various processes of SC, and the pursuit of multiple objectives 17 (minimizing the costs and delays, maximizing service levels and profit …) make the study of SCs a complex task. 18 Consequently, there is a need to use methods and effective modeling and analysis of these systems tools.
The objective of this article is to propose a modeling approach and evaluating the performance of a stock management system by combining two models: BDSPNs as a very powerful tool of modeling and the Supply Chain Operations Reference (SCOR) model, proposed by the Supply Chain Council (SCC).
SCOR is a modeling tool that defines an approach, processes, indicators, and best practices to represent, assess, and diagnose the SC. 19
In this article, we will first show the interest of using BDSPNs and SCOR model for modeling, assessment, and inventory management systems performance analysis, which plays a major role in business-organized SCs. Second, we will present our model and inventory management system. Finally, we will present the performance indicators, calculated by a resolution of stochastic process associated with this system.
Research background
Batch Deterministic and Stochastic Petri Nets
Introduction
There are several classes of SPN models proposed for modeling and performance evaluation of SCs, such as SPNs, GSPNs, 12 and DSPNs. 16 PNs introduced by Petri (1962), as a graphical and mathematical tool, but not be used for modeling and analyzing complex systems which can be characterized as deterministic and/or stochastic. 20 –22 SPNs emerged as a modeling formalism for performance analysis in the early 1980s. In an SPNs, it associates an exponentially distributed delay with the firing of each transition. 23,24 GSPN is an extension for SPN, with two additional modeling elements: immediate transitions and inhibition arcs. 12,25 A DSPN is an extension for GSPN, which adds deterministically timed transitions. 16
In this section, we will introduce the BDSPNs as a new class of SPNs, which is an extension of DSPNs by introducing batch places and batch tokens. Systems such as SCs, purchasing, production, and distribution are usually performed in a batch way with customer orders playing an important role, which justify the motivation for using the BDSPNs. In this last, we can distinguish between two types of places: discrete places and batch places. For the discrete places, tokens (which are the same as those in standard PNs) are viewed indifferently, while tokens in a batch places (which have sizes) are viewed as different individuals, called batch tokens. 26
BDSPNs are characterized by: They may have inhibitor arcs and marking-dependent weights of arcs. If a transition has an inhibitor arc, it is enabled only if the number of tokens in the input place of the arc is less than the weight of the arc and the firing of the transition does not remove any token from the place. They have three types of transitions: immediate transitions, deterministically timed transitions, and stochastic transitions with exponentially distributed firing time. For an immediate transition, tokens will be created in their output places as soon as it is fired, while tokens will be created either after a deterministic time or a time randomly generated according to an exponential distribution. The transition firing conditions and the state evaluation of BDSPNs are similar to those of DSPNs except in BDSPNs where the transitions may be fired in a “batch” way when batch tokens are involved.
Performance evaluation
The performance evaluation is based on the temporal evolution of the model μ-marking process.
26,27
Two main approaches may be possible in the case of deterministic and stochastic Petri networks batches: An analytical approach based on the graph of μ-markings used in conjunction with stochastic processes in order to benefit from the evaluation methods of stochastic processes. This approach is particularly applicable in the case where the network generates a finite number of states (μ-marking); An approach based on discrete event simulation. This technique is particularly useful in the case where the system is too complex for analytical study or when the number of states generated by the model is infinite.
In our study, we applied the analytical approach to evaluate analytically the performance of our system modeled by the BDSPN model. Indeed, the general procedure is given by the algorithm described in Table 1. 26
Algorithm of the general procedure of analysis for the evaluation of performance.
BDSPN: Batch Deterministic and Stochastic Petri Net.
To calculate the indices of performance of a BDSPN, it is necessary to calculate the state probabilities in permanent regime, that is to say, the probability of being in a specific μ-marking. For this calculation, the BDSPN must be limited, and its graph of μ-marking has a finite state space. The type of the resulting stochastic process must be identified in order to be resolved by an appropriate method appropriated to the nature of this process. Thus, the resolution of the latter determines the stationary distribution of probabilities πi of states si which allows the evaluation of performance indices of the modeled system.
In general, if r(µ) represents an “index” function, then the index mean is:
Unlike other class of stochastic Petri networks, BDSPNs can capture the performance metrics that characterize their two types of components: batch and discrete components.
Performance indicators relative to batches components are:
The average number NL of batch tokens in a batch place p ∊ Pb is given by:
where Card(µi ) is the number of batch tokens available in the p at µi state, and Pb means all batches places.
The average number NL(b, p) of batch tokens with size equal to b in a batch place p ∊ Pb is given by:
where Cardb(µi ) is the number of batches tokens with size equal to b available in the place p at µi state.
The average sum of the sizes of batches tokens of a batch place p ∊ Pb corresponding to the average M-marking M(p) Moy of this place is given by:
Let S(tj [q]) be set of μ-marking in the graph GµA(BDSPN, µ 0), where the batch transition tj is crossed by a batch crossing index q. The average frequency of crossing tj [q] is given by:
Modeling of inventory management systems
Whatever the activity of the company, its size, and its organization, stocks exist. Free production stock is almost inconceivable given the many functions performed by stocks.
28
The competitiveness of the company can be particularly affected by its inventory management,
29,30
which constitutes a sufficient reason to pay close attention. Indeed, stockpiling is necessary for many strategic and operational reasons. Supply means ensuring delivery of programming needs and inventory as part of the general planning of the company. Define a supply policy consists essentially on answering three questions: WHAT (what product) should you supply? WHEN should you supply it? HOW should you supply?
In theory, any combination of answers to these three questions determines an inventory management policy that could be adopted for the management of a particular item.
In practice, some combinations have advantages over others and are therefore needed more frequently. In particular, policies identified as the most used in practice inventory management are described in Table 2.
Most used policies in practice inventory management.
It is obvious that each policy is adapted to a product or product category. This often leads companies to use these policies simultaneously. The challenge is to choose the best policy for every product that prevents stock-outs without significant financial asset. Each procurement policy has its advantages and disadvantages. As advantage, we particularly appreciate the simplicity of the procedure, while for disadvantages, it will highlight the costs of the proceedings and especially the risk of possible break. Therefore, each method is ensuring that the disadvantages of this method are acceptable for the company. 28,31,32
Having chosen one of the policies described above, there has yet to set the parameter values that define the system (i.e. the values of s, S, R, and Q) so as to minimize the cost function in question. This optimization problem is generally complex, and thus, it is solved by using heuristic methods. 33
In this work, we are interested in modeling a stock management system in an automotive wiring business, with model based on a stock management policy with continuous review (s, S). The BDSPN model in Figure 1 models a stock managed by the policy (s, S), with ongoing review.

BDSPN model management policy (s, S). BDSPN: Batch Deterministic and Stochastic Petri Net.
In this model, the stock is modeled by the place p 1 and level of p 1 stock, M p 1 plus supply current orders, and M p 3 is modeled by the place p 2.
The output of the stock is modeled by the transition t 1 that we have seen exponential in these representations. It describes customer demand following a Poisson process. t 3 transition models the operation of the award of procurement orders from stock and t 2 transition models the supply operation of the stock. t 3 transition is immediate (crossing time is 0), which specifies that the revision of the stock is continuous.
In other words, the award of procurement orders is done immediately when the stock status is considered critical. Holding the stock is ensured thanks to the following: Place p
2 photographing with his M-marking M(p
2), the level of stock M(p
1), and supply current orders M(p
3). The arc inhibitor is combined to place p
2 whose weight corresponds to the threshold of a supply order. It thus controls the level of inventory, M(p
2), p
1 stock.
t
3 transition and its corresponding arcs modeling model triggering procurement order. The weight associated with these arcs corresponds to the quantity of the command to execute depending on the modeled policy.
SCOR model
Review of SCOR model
Every business process has different characteristics. Thus, a good understanding of the SCM processes is necessary before developing a simulation model. 34 In the past, different metrics were used to measure performance at different SC levels, and models for decision-making at strategic level were scarce. 35 It is not achievable to have a model that perfectly represents SCM but a closely adapted model 36,37 SCOR is achievable. In 1996, the SCC has developed and released a structured framework model SCOR for SCM systems and practices. 38 SCOR model is the first cross-industry framework for evaluating and improving SC performance and management. 35,39 The benefit of the SCOR model is the provision of a standard format and a comprehensive methodology to facilitate communication and to improve SC operations. 36,40 In order to achieve the desired performance in the internal and external levels, upper management of a company is highly recommended to use SCOR to design and reconfigure its SC, as it constitutes a common language and a flexible framework. 35 SCOR is widely used in academia and practice. For these reasons, we use the SCOR framework as a base of our analysis.
SCOR is a business process framework that spans from the supplier’s supplier to the customer’s customer. 19 SCOR 11.0 describes an SC by four levels of details. Level 1 is a top level that defines the scope and high-level configuration by six core processes, that is, plan, source, make, deliver, return, and enable. Level 2 is a configuration level and processes at this level, along with their positioning, determine the SC strategy. Level 3 is a process element level and describes the steps that need to be performed to execute all the processes of the level 2. Level 4 is the implementation level and describes industry-specific activities that are required to perform level 3 processes.
For the SCOR’s components, each one constitutes an important intraorganizational function and a critical interorganization process. 41
SCOR core processes are defined below and shown in Figure 2. 19 The process reference model utilized in our research is based on SCOR 11.0.

SCOR model.
SCOR sourcing process
The sourcing processes are critical to many industrial companies. The industrial company is looking for this process to expect deliveries (receiving, inspection, storage/transfer payment); to identify and select its sources, to define and organize the supply management rules (both in terms of evaluation of supplier performance as maintenance supplies data); and to manage its product arrivals, stocks, infrastructure network providers, export–import constraints, and supplier contracts.
The sourcing processes at level 2 include source stocked product (S1), source make-to-order product (S2), and source engineer-to-order product (S3). There are 17 process elements at level 3, including S1.1, S1.2, S1.3, S1.4, S1.5, S2.2, S2.3, S2.4, S2.5, S3.1, S3.2, S3.3, S3.4, S3.5, S3.6, and S3.7.
This research selects level 2 of the SCOR sourcing process “source make-to-order product (S2)” to modeling, evaluation, and analysis of the performance of our inventory management systems.
The detailed codes of the SCOR sourcing processes source make-to-order product (S2) at level 3 are as follows: Schedule product deliveries S2.1 Receive product S2.2 Verify product S2.3 Transfer product S2.4 Authorize supplier payment S2.4
Table 3 shows the performance metric codes of the SCOR sourcing processes source make-to-order product (S2) at level 2; there are 17 performance metrics for this level.
Performance metric definition of SCOR sourcing “make-to-order product” process at level 2.
SCOR: Supply Chain Operations Reference; GRC: Governance, Risk, Compliance.
After defining performance indicators sourcing processes source make-to-order product (S2) established by the model SCOR, we modeled, analyzed, and assessed some of these indicators based on the dynamic model BDSPNs through a case study presented in the next section.
Case study
In this section, an analytical study is presented on a stock management system modeled by a BDSPN. This study is based on the graph of BDSPN μ-markings used in conjunction with the associated stochastic process. The resolution of the stochastic process used to determine the probabilities of the various states and the use of these probabilities assesses the average performance of the system.
The BDSPN model of this system is limited (finite state) and shows no mathematical constraints for the characterization and study of the associated stochastic process.
System description
The storage system considered in this study is shown in Figure 3 and has the following characteristics:

Description of the storage system to study.
The stock is managed by the policy (s, S) with continuous review.
The stock of the raw material is modeled by p 1 instead.
The production hall is modeled by a set comprising: t 1: ready production orders; p 5: production orders in process; p 4: products required to produce a batch depending on the size of the production order; and t 4: production operation.
Orders productions are in a batch mode. They are generated by following different sizes.
The lead time of the stock is stochastic. It follows an exponential law.
Orders productions (batches) arrive via the transition t 4 and begin waiting in the lot instead p 6.
Modeling system
The storage system to be studied is modeled by the BDSPN in Figure 4. The input and output of the stock instead modeled by p 1 are represented by t 1 and t 2 transitions, respectively. The stock is managed by the policy (s, S) with continuous revision.

BDSPN model of the storage system to study.
Recall, it is to inspect the stock level continuously and when it falls to a predetermined level s (weight of inhibitor arc from p 2 to t 3), a supply control of the stock is on to try to recomplete the stock at a level set at S. 42
The size of the procurement order is determined by the variable weight, S − M(p 2), the arc connecting p 3 to t 3. The stock, p 1, is subject to random requests for uses by the production workshop and sizes. Requests for applications are randomly generated by the transition t 1 lot following a Poisson process of rate λ 1. The transition t 1 is exponential since in a Poisson process, the interarrival dates are distributed exponentially.
If the stock is available, the execution of a request of use is by crossing the batch transition t 1. Crossing at a time t 1 decreases the M-marking p 1 site (stock) and that of p 2 (representing the position of the stock, M(p 2) = M(p 1) + M(p 3)) up a lot amount of size equal to the size of the production order to produce.
Once the stock level (M-marking p 2) falls below the thresholds, monitored by the inhibitor arc connecting the transition t 3 in the place p 2, a supply order is immediately passed with the crossing of the immediate transition t 3 (continued inspection of stock). Crossing the transition t 3 adds in order with quantity (batch) of S − M(p 2) in the batch place p 3. This batch place represents the supply orders in execution waiting in order to attempt to bring back the level of stock to S.
The supply operation is executed by crossing the transition t 2 whose crossing delay times are exponentially distributed.
The parameters of the policy management of the stock (s, S) are fixed to (2, 4). We consider that there are only two possible sizes (1 or 2) concerning production orders, where the μ-marking of batch place p 6, modeling the production orders in expectations, is µ(p 6) = {1, 2}.
The initial μ-marking is: µ 0 = (4, 4, Ø, 0, Ø, {1, 2}).
Evolution of the modeled system
The evolution of the system can be expressed by using the graph of μ-markings BDSPN model that represents this system.
Each μ-marking of the graph represents the state of the system and each crossing between two μ-markings in the graph represents the execution of an operation that can be the award of a supply (crossing t 3) and supply stock (crossing t 2).
The graph of μ-markings of BDSPN is shown in Figure 5. This graph is built from the initial μ-marking µ 0 considering all possible crossings of a μ-marking to another. To recall, the immediate transitions are prior than timed transitions. The resulting graph contains 15 tangible states numbered from 0 to 14.

μ-Markings graph model.
The μ-marking graph generated by the BDSPN of the modeled system is isomorphic to a continuous time Markov chain shown in Figure 6. 2,3 The nature of the stochastic process generated is verified by the nature of distribution laws associated with the transitions used (without memory law).

The Markov chain associated.
In other words, during the evolution of the model, only the new μ-marking determines the future behavior of the stochastic process.
Resolution of the associated stochastic process
For the rest of the study of this system, we will be limited to the following case: λ 1[1] = λ 1[2] = λ 1 and λ 3[3] = λ 3[4] = λ 3.
λ 1[1] = λ 1[2] = λ 1: means that the arrival of the two types of production orders follows two exponential laws of identical parameters.
λ 3[3] = λ3[4] = λ 3: means that the supply delay times don’t depend on the size of the supply orders.
t 4 transition is exponentially distributed with a parameter #λ 2. The symbol # indicates that λ 2 parameter of the transition t 4 depends on the size of the batch that crosses this transition: λ 2[2] = λ 2 means that λ 2[1] = 2λ 2.
In this case, the generator of the Markov process Q (transition matrix) associated with the stochastic μ-markings process is given in Table 4.
Transition matrix.
The probability distribution of the states (μ-marking), depending on the parameters λ
1, λ
2, and λ
3; the permanent regime of system is obtained by solving the following system:
System performance evaluation
By knowing the probability distribution of the states (μ-marking), it is possible to assess several performance of the modeled system using the performance indices associated with BDSPN model. In the following, the system performances are given as functions that depend on the three parameters λ 1, λ 2, and λ 3. This will allow us to analyze the effect of these parameters on the variation in system performance.
Average stock
The function of the average stock Sa based on parameters λ 1, λ 2, and λ 3 is given by the M-marking means and the place p 1. It is calculated by the following formula:
Average cost of stock
Knowing the unit stock price Cs (price of storing a product unit in a unit of time in the stock), the function of the average storage costs Cs a based on parameters λ 1, λ 2, and λ 3 is given by multiplying the average stock by the unit cost:
Probability of an empty stock
On the μ-markings graph of BDSPN, μ-markings of the place p 1 is 0 (µ = 0) in the four μ-markings μ 8, μ 12, μ 13, and μ 14. Thus, the function of the probability of an empty stock is given by:
Average supply frequency
As the supply of stock is done with the clearing of the batch transition t
2, we observe on the graph of μ-markings that this transition is in the μ-markings:
µ
6, µ
7, µ
9, µ
13, and µ
14 with a batch of size equal to 3 (crossing t
2[3]) corresponding to the size of a supply order.
µ
8 and µ
12 with a batch of size equal to 4 (crossing t
2[4]) corresponding to the size of a supply order.
The average frequency of supply stock based on parameters λ 1, λ 2, and λ 3, noted FSa (λ 1, λ 2, λ 3), corresponds to the average frequency of crossing the transition t 2 and thus, the sum of the average frequency crossing, F(t 2[3]) and F(t 2[4]), and since we assumed in the statement of this study that λ 3[3] = λ 3[4] = λ 3, we get:
For our industrial application, the settings of the stock system that we have chosen are S = 10; s = 2; λ 1 = 10; λ 2 = 4; λ 3 = 1; Cs = 2.
The probability distribution of the states (µ 0, µ 1, µ 2, µ 3, µ 4, µ 5, µ 6, µ 7, µ 8, µ 9, µ 10, µ 11, µ 12, µ 13, µ 14) is the vector [π 0, π 1, π 2, π 3, π 4, π 5, π 6, π 7, π 8, π 9, π 10, π 11, π 12, π 13, π 14] = [0.0083, 0.0127, 0.0166, 0.0346, 0.0178, 0.0384, 0.0576, 0.0178, 0.0296, 0.0311, 0.0210, 0, 0.0593, 0.0311, 0.6242].
We also presented the settings of the stock system for our industrial application to calculate these indicators, we obtained: Average stock = 0.5112 Average cost of stock = 1.0224 Probability of an empty stock = 0.7442 Average supply frequency = 0.8507
According to these results, some indicators need improvement to achieve the objectives dictated by the company. For example, the “probability of an empty stock” is very high , thus, there is a necessity to search the best combination of the three parameters (S, s, and λ 3) to decrease this indicator.
Conclusion
In this article, we achieved our objective of modeling and evaluation of the performance of a stock management system by combining two models: BDSPNs model and the SCOR model.
The main idea of this work is to start with the definition and measurement of performance indicators based on a reference static model, appropriate and dedicated to the field; it is the SCOR model that offers a large number of indicators. These indicators provide a range of measures to evaluate the performance of an inventory management system. Then, we modeled and simulated our stock management system through a case study by a high-performance dynamic model, namely BDSPNs. The rich BDSPNs model allowed us to easily assess important performance indicators of the stock management system that we had already defined through the SCOR model.
In the second section, we presented a research background about the two models. The first part of this section was about the BDSPNs. First, we introduced this model and its characteristics. In addition, we presented the approaches used for performance evaluation. In our study, we adopted the analytical approach to evaluate analytically the performance of our system and we presented the algorithm of the general procedure. We also described the most used policies in practice inventory management. As we are interested in modeling a stock management system in an automotive wiring business, we selected the stock management policy with continuous review (s, S). Then, we presented the BDSPN model management policy (s, S) and its characteristics. The second part concerns the SCOR model in which we presented a review of this model, its important contribution to the SCM, and the reasons of its choice as a base for our analysis. We also described the SCOR sourcing process. In our study, we selected the level 2 of the SCOR sourcing process source make-to-order product (s2) to modeling, evaluation, and analysis of the performance of our inventory management systems.
The third section of this article concerns the case study. An analytical study is presented on a stock management system modeled by the BDSPNs. This study is based on the graph of BDSPNs μ-markings used in conjunction with the associated stochastic process. First, we described the storage system and its characteristics. Then, we modeled this system by the BDSPNs. We also expressed the evolution of the system by using the graph of μ-markings BDSPNs model that represents this system. After this, we presented the resolution of the associated stochastic process; this resolution is used to determine the probabilities of the various states and the use of these probabilities assesses the average performance of the system. Finally, we evaluated the performances of our system by the following indicators: average stock, average cost of stock, probability of an empty stock, and average supply frequency.
These results revealed that some indicators need improvements to achieve the objective values required by the company.
In this work, we focused on the analytical approach for the modeling and performance evaluation of a stock management system. As future work, we will use techniques of performance evaluation based on simulation in the case of too complex systems.
Footnotes
Author’s Note
Jamal Fattah is now affiliated to Modeling, Control Systems and Telecommunications team, High School of Technology of Meknes (ESTM), Moulay Ismail University P.O. Box 3103, Toulal, Meknes 50000, Morocco.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
