CN115841165A - Crude oil scheduling optimization method and device applied to oil refinery - Google Patents
Crude oil scheduling optimization method and device applied to oil refinery Download PDFInfo
- Publication number
- CN115841165A CN115841165A CN202111102934.2A CN202111102934A CN115841165A CN 115841165 A CN115841165 A CN 115841165A CN 202111102934 A CN202111102934 A CN 202111102934A CN 115841165 A CN115841165 A CN 115841165A
- Authority
- CN
- China
- Prior art keywords
- oil
- tank
- crude oil
- long
- pipeline
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a crude oil scheduling optimization method and a crude oil scheduling optimization device applied to an oil refinery, which belong to the technical field of process industrial production scheduling optimization and comprise the following steps: step 1, establishing a crude oil scheduling optimization model based on a discrete time strategy; step 2, based on the established mathematical model, flexibly defining model input parameters and manual intervention events; step 3, according to a divide-and-conquer solving strategy, a crude oil scheduling model is optimized and solved to obtain a crude oil scheduling scheme and crude oil physical property values at all times in an event period; step 4, checking the physical property values of the crude oil in the step 3: checking the physical properties of the crude oil stored in the tank at each moment and the physical properties of the atmospheric and vacuum processed crude oil, judging whether the physical properties pass the checking, if the physical properties pass the checking, executing the step 5, and otherwise, returning to the step 2 to check and adjust the input parameters and the manual intervention event; and 5, forming a complete solution in the scheduling event period and outputting a report. The method comprehensively considers each decision point of the crude oil scheduling problem, has short solving time and is suitable for popularization and application.
Description
Technical Field
The invention belongs to the technical field of scheduling optimization of process industrial production, relates to a crude oil scheduling optimization method and device applied to an oil refinery, and particularly relates to a crude oil scheduling optimization method and device based on a mixed integer programming and divide-and-conquer solving strategy.
Background
Crude oil scheduling is the front end of refinery units and occupies an important position in the whole refinery scheduling, but research on crude oil scheduling is just started and a large blank exists. The main problems currently studied for crude oil scheduling include arrival of crude oil at a port (reservoir), oil discharge operation, pipeline transportation, oil recovery from a tank in a factory, atmospheric and vacuum equipment, and the like. The optimization objective is to minimize operating costs. In the aspect of model solution, because the existing solution software has great difficulty in solving mixed integer nonlinear programming, a simplified processing method is usually adopted to linearize a nonlinear problem in the solution, and the nonlinear problem is changed into a mixed integer programming problem (MILP) to solve.
In the prior art, a Chinese patent with a publication number of CN101286065 discloses a crude oil blending scheduling method for crude oil blending multi-period optimization, which summarizes a mathematical model of a problem; aiming at a nonlinear mixed integer programming model formed by problem regression, a satisfactory solution of the problem is obtained by utilizing an efficient mixed solving strategy. The scheme mainly aims at improving the stability of the properties of the blended crude oil, and simultaneously aims at improving the yield of the expected product, so that the influence on the stability of a distillation device and subsequent devices thereof caused by overlarge property fluctuation of the blended crude oil can be reduced, and the yield of the expected product can be increased to a certain extent. The algorithm structure for solving the discrete decision variables and the continuous decision variables in a layered mode is adopted, and the problems of low algorithm efficiency, poor search performance and the like caused by simultaneous optimization of two groups of variables are solved.
Chinese patent publication No. CN103984990A discloses a modeling method based on discrete time scheduling for the whole refinery, which divides the whole refinery system into three parts of crude oil supply, refinery production, and finished oil blending delivery, based on discrete time, modeling is performed from the perspective of the operation mode of the production apparatus and the transition process of the operation mode of the production apparatus, based on mode switching in the production scheduling for various finished oils of the refinery enterprise and discrete time optimization operation control of the transition process, the scheduling control for the whole refinery is given, and a scheduling model capable of minimizing the production cost and material storage cost of the production process and violating order penalties is constructed.
Although the existing method solves certain technical problems, the scale of the model is limited to a certain extent due to the problem of algorithm efficiency, and the accuracy of solution is also influenced to a certain extent due to the inevitable introduction of errors caused by linearization. In addition, the modeling and solving of the whole segment of the current crude oil scheduling often simplifies certain important business scenarios and even does not consider the scenarios (such as the transmission of a long-distance pipeline oil head, the implementation of rolling scheduling and the like). In the prior art, the technical problems of large solving scale and long solving time of crude oil scheduling in an oil refinery exist, and a crude oil scheduling optimization method and a crude oil scheduling optimization device which are efficient, close to a service scene and capable of realizing rolling scheduling and are applied to the oil refinery are urgently needed.
Disclosure of Invention
The invention aims to solve the technical problems and provides a crude oil scheduling optimization method applied to an oil refinery, which has short solving time and meets the requirements of industrial practical application.
In order to achieve the purpose, the invention mainly provides the following technical scheme:
in one aspect, the embodiment of the invention provides a crude oil scheduling optimization method applied to an oil refinery,
the method comprises the following steps:
step (1), establishing a crude oil scheduling optimization model based on a discrete time strategy;
step (2), based on the mathematical model established in the step (1), flexibly defining model input parameters and manual intervention events;
step (3) according to a divide-and-conquer solving strategy, a crude oil scheduling model is optimized and solved to obtain a crude oil scheduling scheme and crude oil physical property values at all moments in an event period;
and (4) checking the physical property values of the crude oil in the step (3): checking the physical properties of the crude oil stored in the tank at each moment and the physical properties of the atmospheric and vacuum processed crude oil, judging whether the checking is passed, if the checking is passed, executing the step (5), and otherwise, returning to the step (2) to check and adjust the input parameters and the manual intervention event;
and (5) forming a complete solution in the scheduling event period, and outputting a report, wherein the scheduling event summary table can be used as the input of the scheduling of the next period to realize rolling scheduling.
Specifically, in step (1), the crude oil scheduling optimization model of the discrete time strategy comprises:
material balance model of crude oil transportation tanker: the oil unloading amount of the crude oil tanker in the period is equal to the carrying amount of the crude oil tanker in the period;
tank field storage tank material restraint model: the tank oil storage amount at the end of each time in the storage tank period is equal to the tank storage amount at the end of the last time plus the tank entering amount at the moment minus the tank exiting amount at the moment, the tank entering amount at the moment comprises the oil tanker oil unloading and tank entering amount and the other storage tank pipeline transmission and tank entering amount, and the tank exiting amount at the moment comprises the tank exiting amount transmitted to other storage tanks or an atmospheric and vacuum device by pipelines;
selecting a model for an atmospheric and vacuum device scheme: and selecting a processing scheme at each moment of the atmospheric and vacuum distillation unit, and processing the crude oil according to the processing scheme, wherein the total processing amount of the oil seeds in a period is more than or equal to the required amount of the oil seeds.
Further, the storage amount of the storage tank is kept within the range of the upper limit and the lower limit of the storage amount of the safety tank;
the pipeline transmission follows a topological structure, the transmission quantity is within the range of the upper limit and the lower limit of the pipeline transmission, and the transmission oil seeds are matched with the oil seeds of the starting tank and the target tank;
the oil head is considered in the long-distance pipeline transmission, and the oil type and the oil quantity at the inlet of the long-distance pipeline at the current moment are equal to the oil type and the oil quantity at the outlet of the long-distance pipeline at the current moment plus the corresponding moment after the long-distance pipeline transmission time.
The optimization goal of the crude scheduling optimization model is to minimize the total cost (minimum production fluctuations) throughout the crude scheduling process.
The expression of the total cost obtained by adding the overdue cost of the oil tanker, the storage cost of the tank area, the switching cost of the atmospheric and vacuum oil mixing scheme, the switching cost of the oil type of the tank area and the fluctuation cost of the pipeline is as follows:
the total cost = tanker overdue cost + tank farm storage cost + atmospheric and vacuum oil mixing scheme switching cost + tank farm tank oil species switching cost + pipeline fluctuation cost.
Introducing human intervention events and flexible parameters to improve the practicability of the step (1) model.
Specifically, in step (2), the parameters and the manual intervention event are specifically:
the parameters include:
oil tanker input parameters: arrival time, estimated departure time, type and quantity of transported crude oil, arrival dock, unit overdue cost;
storage tank input parameters: the storage tank is located in a tank area, the upper limit and the lower limit of the physical tank capacity, the upper limit and the lower limit of the safe tank capacity, the current oil type and the stock, the initial state (oil receiving, oil paying and standing), the standing time of the initial state, the inventory cost, the physical property and whether the oil type is allowed to be switched or not;
pipeline input parameters: pipeline starting point, pipeline finishing point, pipeline transmission capacity upper and lower limits, pipeline transmission cost and long-distance pipeline oil head information (pipe section, oil type and oil quantity);
the device inputs parameters: the crude oil proportioning is selected according to the initial scheme of the atmospheric and vacuum distillation unit;
the requirement input parameters are as follows: the periodic demand of each oil;
the human intervention event comprises: crude oil enters a tank (oil is unloaded to a storage tank by an oil tanker), is transferred to the tank (the tank is transferred to the storage tank without a long-distance pipeline), is transferred by a crude oil pipe (a wharf tank field is transferred to the long-distance pipeline), and is switched according to the oil type (the storage tank stores the oil type); the human intervention event may also be derived from a scheduled event output scheduled in the last cycle.
The crude oil scheduling model in the step (1) adopts a traditional solving method to solve for a long time, so that the practical application is difficult, and the solving time can be greatly shortened by introducing a divide-and-conquer solving strategy to carry out optimization solving.
Specifically, in step (3), the divide-and-conquer solution strategy includes the following sub-steps:
step 31), module segmentation: the model is divided into three modules by taking a long transmission pipeline as a node: before the long-distance pipeline (oil tanker-dock tank field-long-distance pipeline inlet), after the long-distance pipeline (long-distance pipeline inlet-long-distance pipeline outlet), after the long-distance pipeline (long-distance pipeline outlet-atmospheric pressure relief device);
step 32), solving before the long-distance pipeline (oil tanker-wharf tank field-long-distance pipeline inlet): adopting a strategy of integral modeling and solution emphasizing, neglecting the oil head transmission of the long-distance pipeline, assuming that the transmission time of the long-distance pipeline is fixed time, establishing an integral crude oil scheduling model, wherein an objective function comprises global important cost, but the cost more emphasizing the front of the long-distance pipeline (oil tanker-wharf tank area-long-distance pipeline inlet) is the minimum, and obtaining global feasible solution and optimal solution/near optimal solution of the front of the long-distance pipeline (oil tanker-wharf tank area-long-distance pipeline inlet);
step 33), solving the long-distance pipeline (inlet of the long-distance pipeline-outlet of the long-distance pipeline): simulating the inlet oil type and the oil quantity of the long-distance pipeline obtained in the step 32) according to the transmission rate obtained in the step 32) by adopting a simulation strategy and considering oil head transmission to obtain the outlet oil type and the oil quantity of the long-distance pipeline at each moment;
step 34), solving after long-distance pipeline transportation (outlet of the long-distance pipeline-atmospheric and vacuum device): modeling and solving are carried out only for the long-distance pipeline outlet-atmospheric and vacuum device by adopting a module modeling and solving strategy, the cost of the long-distance pipeline outlet-atmospheric and vacuum device is the minimum by an objective function, and the optimal solution/near-optimal solution behind the long-distance pipeline (the long-distance pipeline outlet-atmospheric and vacuum device) is obtained;
and step 35) combining the solution of the front part (oil tanker-wharf tank area-long conveying pipeline inlet) of the long conveying pipeline in the global variable in the step 32) and all the solutions of the rear part (long conveying pipeline outlet-atmospheric and vacuum unit) of the long conveying pipeline in the step 34) to form the solution of the whole crude oil scheduling.
And (4) the solution obtained in the step (3) needs to be verified through physical property to be identified as an effective scheduling instruction which can be applied.
Further, in step (3), slack is also introduced into the divide-and-conquer solution strategy, which includes: tank storage relaxation, crude oil production plan relaxation, and storage reduction relaxation (tank storage reduction level in the tank area at the end of the term breaks through a set value).
And (4) inputting the solution obtained in the step (3) to form a report after the verification in the step (4).
Specifically, in step (5), the report includes:
tank farm oil seed balance sheet: the balance table of the inlet and outlet of each oil type in each tank area in the whole period;
tank field storage tank balance table: the inlet and outlet balance tables of all storage tanks of all tank areas in the whole period;
tank trend table: storing oil seeds and oil mass in each storage tank at each moment;
scheduling event summary table: the scheduling event summary table comprises a crude oil tank inlet table, a crude oil pipe transmission table, a crude oil tank transfer table, an overland transportation table, a crude oil processing table, a tank oil type switching table, a physical property output table and the like, and can be used as an input (a manual intervention event) of the next periodic scheduling, so that the rolling scheduling is realized, and the following specific contents of each event type are as follows:
crude oil enters a tank: refers to a tanker that transports crude oil to a port tank, including start time, end time, oil type, oil quantity, start point (tanker lot) and end point (port tank number);
crude oil pipe transportation: the method comprises the steps that a wharf tank conveys crude oil to an in-plant tank through a long-distance pipeline, and the in-plant tank comprises starting time, ending time, oil type, oil quantity, a starting point (the wharf tank) and a finishing point (the in-plant tank);
transferring crude oil to a tank: refers to the mutual transmission between other crude oil storage tanks besides the crude oil pipeline transmission, including the start time, the end time, the oil type, the oil quantity, the starting point (oil storage tank) and the end point (oil storage tank);
land transportation: the method comprises the steps that a land transportation starting point conveys crude oil to an in-plant tank through a land transportation pipeline, and the land transportation starting point comprises starting time, ending time, oil type, oil quantity, a starting point (land transportation starting point) and an end point (in-plant tank);
crude oil processing: the method is characterized in that the in-plant tank conveys crude oil to an atmospheric and vacuum distillation unit for processing, and comprises an atmospheric and vacuum distillation unit processing scheme, a processing amount, a crude oil type and crude oil quantity conveyed by the in-plant tank, a starting time, an ending time, a starting point (in-plant tank) and an end point (atmospheric and vacuum distillation unit);
switching of the oil types of the tanks: the switching of the oil storage tank in the oil storage is carried out, and comprises switching time, oil types after switching and tank stock after switching;
physical property output table: comprises an atmospheric and vacuum distillation unit processing crude oil physical tracking list and a storage physical tracking list.
In another aspect, an embodiment of the present invention further provides a crude oil scheduling optimization apparatus applied in an oil refinery, including:
the model establishing unit is used for establishing a crude oil scheduling optimization model based on a discrete time strategy;
the input unit is used for supporting flexible definition of model input parameters and manual intervention events based on the mathematical model established by the model establishing unit;
the model optimization unit is used for optimizing and solving the crude oil scheduling model according to the divide-and-conquer solving strategy to obtain a crude oil scheduling scheme and crude oil physical property values at all times in an event period;
the physical property checking unit is used for checking the physical property value of the crude oil obtained by the model optimizing unit;
and the output unit is used for forming a complete solution in the scheduling event period and outputting a report.
In another aspect, an embodiment of the present invention further provides a crude oil scheduling optimization system applied in an oil refinery, including:
the method comprises the following steps: one or more processors; a memory for storing one or more programs; the processor is configured to execute program instructions stored in the memory that when executed perform the above-described crude oil scheduling optimization method applied to a refinery.
In another aspect, the present invention also provides a computer-readable storage medium, on which a computer program is stored, wherein the computer program, when executed by one or more processors, implements the above-mentioned crude oil scheduling optimization method applied to a refinery.
The crude oil scheduling optimization method and device applied to the oil refinery comprehensively consider all decision points of the crude oil scheduling problem, are high in conformity with actual services, and are good in performability; the technical scheme provided by the invention provides an optimized solving strategy, which greatly improves the operation time, has short solving time, meets the actual production execution requirement and is suitable for popularization and application.
Drawings
FIG. 1 is a flow chart of a crude scheduling optimization method of the present invention as applied to a refinery;
FIG. 2 is a view showing the topology of a tank farm pipeline according to embodiments 1 and 2 of the present invention;
FIG. 3 is a tank farm pipeline topology structure diagram of embodiment 3 of the present invention;
FIG. 4 is a schematic flow chart of a divide-and-conquer solution strategy according to an embodiment of the invention;
FIG. 5 is a schematic diagram of a crude scheduling optimization unit for use in a refinery according to the present invention;
FIG. 6 is a graph of the density trend for example 1 (examples CDU2 and H116);
FIG. 7 is a trend trace of sulfur content for example 1 (examples CDU1 and H111);
FIG. 8 is an overview of the dispatch Gantt chart of example 1 (example atmospheric and vacuum device);
FIG. 9 is a Gantt chart section (e.g., atmospheric and vacuum distillation apparatus) of example 1;
FIG. 10 is a plot of the density trend for example 2 (examples CDU2 and H116);
FIG. 11 is a trend trace of sulfur content for example 2 (examples CDU1 and H111);
FIG. 12 is a dispatch Gantt overview (example atmospheric and vacuum device) of example 2;
FIG. 13 is a part of a dispatching Gantt chart (e.g., atmospheric and vacuum relief device) of example 2;
FIG. 14 is a plot of the density trend for example 3 (examples CDU2 and H116);
FIG. 15 is a trend trace of sulfur content for example 3 (examples CDU1 and H111);
FIG. 16 is a Schedule overview of example 3 (example atmospheric and vacuum;
fig. 17 is a schematic diagram of the schedule chart (example atmospheric and vacuum apparatus) in example 3.
Detailed Description
Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the invention are shown in the drawings, it should be understood that the invention can be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Referring to fig. 1, the invention provides a crude oil scheduling optimization method applied to an oil refinery, which specifically comprises the following steps:
step (1), establishing a crude oil scheduling optimization model based on a discrete time strategy;
step (2), based on the mathematical model established in the step (1), flexibly defining model input parameters and manual intervention events;
step (3), according to a divide-and-conquer solving strategy, a crude oil scheduling model is optimized and solved to obtain a crude oil scheduling scheme and crude oil physical property values at all times in an event period;
and (4) checking the physical property values of the crude oil in the step (3): checking the physical properties of the crude oil stored in the tank at each moment and the physical properties of the atmospheric and vacuum processed crude oil, judging whether the checking is passed, if the checking is passed, executing the step (5), and otherwise, returning to the step (2) to check and adjust the input parameters and the manual intervention event;
and (5) forming a complete solution in the scheduling event period, and outputting a report, wherein the scheduling event summary table can be used as the input of the scheduling of the next period to realize rolling scheduling.
In the specific application process, the method specifically comprises the following steps:
the method comprises the following steps:
step (1) establishing a crude oil scheduling discrete time model, specifically operating as follows:
the material balance model of the oil tanker is appointed as follows:
formula (1) represents that the actual departure time of the tanker is out of date if the expected latest departure time of the tanker is exceeded; equation (2) indicates that the tanker can only transmit to the port tanks of the matched oil species and can only transmit within the valid time window (arrival time — latest departure time); formula (3) shows that the sea transportation volume of the oil tanker in the whole period is equal to the total crude oil carried by the oil tanker; the formula (4) represents the material balance of the oil tanker: the surplus at the end of the oil tanker time = the surplus at the end of the time on the oil tanker minus the amount of the oil tanker entering the plant at the time on the oil tanker.
The tank field storage tank material constraint model is given as follows:
the formula (5) shows that the storage tank in the tank area is in one of three states of oil receiving, oil applying and standing at any time, wherein the oil applying shows that the oil is output and the oil tank is in an available state; formula (6) shows that when the oil tank is used as a crude oil transfer destination tank, the oil tank state is 'oil receiving'; formula (7) shows that the oil tank can enter the oil supply (available) state after the oil collection is finished and the standing time requirement is met, namely the oil tank after oil collectionThe oil supply state can not be entered in the standing time, and if the secondary oil receiving occurs in the standing period, the standing time needs to be recalculated after the secondary oil receiving is finished; the formula (8) and the formula (9) show that the oil tank can be output when in the oil supply state; equation (10) represents the tank equilibrium constraint: the storage tank end-of-period stock = last-time end stock of the storage tank + the amount of the tank entering at the moment-the amount of the tank leaving at the moment; formula (11) shows that when the oil type of the oil tank changes at adjacent time, the switching is determined to be one time; equation (12) indicates that the tank can only be switched between the same type of crude oil type, where k 0 Representing the type set of the oil seeds at the beginning of the oil tank; equation (13) shows that the can cutting operation is only allowed if the amount of cans reaches the lower limit of can storage, and since this expression is non-linear, the linearization yields equation (14). />
The pipeline transport constraint model is given as follows:
formula (15) indicates that the transmission of each pipeline can only occur on the active pipeline and meets the transmission capability limit; the formula (16) shows that the transmission rate of each pipeline meets the limitation of the upper and lower transmission capacity limits; equations (17) - (21) represent the pipeline topology constraints, which are, in order from top to bottom: at most one effective transmission pipeline exists at each moment in the wharf tank area and the commercial tank area, at most one effective transmission pipeline exists at each moment in the wharf tank area, at most one effective transmission pipeline exists at each moment in the commercial tank area, at most one effective transmission pipeline must exist at each moment in each land transportation node, and at most one effective transmission pipeline exists at each moment from the tank area in the factory to the atmospheric and vacuum distillation unit; equations (22) - (24) represent the pipeline oil transport restrictions, which are, in order from top to bottom: the land transportation of crude oil must match the oil seeds of oil tanks collected in a factory, the transportation of a wharf tank field needs to match the oil seeds of auxiliary oil tanks and oil tanks, and the pipeline transportation from an in-factory tank to an atmospheric and vacuum tank field needs to match the oil seeds of the in-factory tank;
the long delivery pipeline transmission constraint model is given as follows:
equations (25) - (26) show that a valid input pipeline and a valid output pipeline must exist in the long-distance pipeline at any time; equations (27) - (28) represent long transfer line head transfers: in-plant tank bay front TL mt-cn The outlet oil type (oil head transmission time) is equal to the pipe section oil type at the moment, and the outlet oil quantity is equal to the pipe section oil quantity at the moment; equations (29) - (30) represent the long haul pipeline terminal tank yard incoming transport: represents TL mt-cn (oil head transfer time) after the time, the outlet oil type of the tank area in the factory at the time t is equal to t-TL mt-cn The oil quantity at the inlet of the long-distance pipeline at the moment is equal to the oil quantity at the inlet; formula (31) shows that when the source tank of the long-distance pipeline at the adjacent moment is changed, the long-distance pipeline is determined to be switched; equations (32) - (33) define the fluctuation amount of the long-distance pipeline: the fluctuation amount is an absolute value of the transmission amount of the pipeline at the adjacent time.
The selection model of the oil mixing scheme of the atmospheric and vacuum device is given as follows:
the formula (34) shows that each atmospheric and vacuum equipment can only adopt one oil mixing scheme at each moment; formula (35) shows that the crude oil feeding amount from the tank area in the plant to atmospheric and vacuum pressure matches the atmospheric and vacuum pressure processing requirement; the formula (36) shows that the total processing amount of various oils of the atmospheric and vacuum equipment meets the requirements of a whole-cycle crude oil production plan; the formula (37) shows that the atmospheric and vacuum equipment adopts different schemes at adjacent moments and the oil mixing scheme switching is determined to occur.
The overall scheduling aims at minimizing the total cost, and the specific expression is as follows:
the first term in the formula (38) is the overdue cost of the tanker, the second term is the storage cost of the tank field, the third term is the switching cost of the atmospheric and vacuum oil mixing scheme, the fourth term is the switching cost of the tank oil type of the tank field, and the fifth term is the pipeline fluctuation cost.
Step (2): and (3) flexibly defining model input parameters and manual intervention events based on the mathematical model established in the step (1).
And (3): and according to a divide-and-conquer solving strategy, optimizing and solving the crude oil scheduling model to obtain crude oil scheduling schemes and crude oil physical property values at all times in an event period.
The divide-and-conquer solving strategy of the embodiment takes the long-distance pipeline as a division point, firstly, the global modeling is carried out, the objective function emphasizes and optimizes the front (oil tanker-long-distance pipeline inlet) part of the long-distance pipeline, and the method comprises the following solving steps:
step 2, operation: considering the tank is loose, the oil seeds can be switched at any time when the tank is cut, and the tank is not allowed to be cut when the tank is not cut. Minimizing the target content: the relaxation cost, the tank oil seed switching cost and the tanker overdue cost;
step 3, operation: and locking the oil types stored in the oil tanks in the wharf tank field and the switching time point. Minimizing the objective function: switching the inlets of the long-distance pipelines and the number of the oil delivery tanks;
step 4, operation: and locking the state of the wharf tank field and the oil conveying tank at the inlet of the long conveying pipeline. Minimizing the objective function: fluctuation cost + mutual transmission cost;
and 5, no operation is performed, and fluctuation is optimized. Minimizing the objective function: long transfer lines fluctuate costs.
And 6, obtaining inlet information of the long-distance pipeline, including time, oil type, oil quantity, physical property and the like.
And for the inlet information of the long-distance pipeline, combining the oil head to perform pipeline section combination and long-distance pipeline simulation to obtain a long-distance pipeline outlet information table, and obtaining the outlet information of the long-distance pipeline after simulation, wherein the outlet information comprises time, oil type, oil quantity, physical property and the like.
And for the outlet information of the long-distance pipeline, establishing a model after the long-distance pipeline (the outlet of the long-distance pipeline-the atmospheric and vacuum device) and solving to obtain a scheduling scheme after the long-distance pipeline. The solving steps are as follows:
step 2, operation: the oil seeds can be switched at any time when the can is cut, and the can is not allowed to be cut when the can is not cut. Minimizing the objective function: relaxation cost + tank cutting cost + scheme switching cost;
step 3, operation: and locking the state of the tank area in the factory and the oil conveying tank at the inlet of the long conveying pipeline. Minimizing the objective function: fluctuation cost + number of oil supply tanks;
step 4, operation: and locking the state of the tank in the factory. Minimizing the objective function: fluctuating cost + scheme switching cost.
And combining the scheduling schemes before and after the long-distance pipeline to obtain a scheduling scheme which is preliminary in the whole crude oil scheduling and is not subjected to physical property verification and a physical property tracking trend chart of the tank storage and atmospheric and vacuum devices.
And (4): and (3) checking the physical property values of the crude oil in the oil tank and the atmospheric and vacuum distillation unit in the step (3).
And (5): and (5) outputting a report if the physical properties in the step (4) meet the requirements.
The parameters, variables, and introduced relaxation in this example are defined as follows:
(1) Parameters and sets
I ship : set of marine crude oil categories, I ∈ I ship 。
I land : set of land-transported crude oil types, I belongs to I land 。
I: crude oil class set, I = I ship ∪I land 。
K: the set of crude types, K ∈ K.
t: the time scale is set, T ∈ T.
H: wharf set, H ∈ H.
F: and F is the same as F.
N land : and (4) carrying out land transport on the crude oil starting point set, wherein N belongs to N.
n long : long distance pipeline node.
C: tank field tank assembly, C = C mt ∪C sc ∪C cn Where mt represents a terminal tank farm, sc represents a commercial storage tank farm, and cn represents an in-plant tank farm.
e: and E belongs to the set of the atmospheric and vacuum equipment.
L: pipeline assembly for crude oil dispatching, comprising tanker (dock) → dock tank field,Dock tank field → long-distance pipeline, long-distance pipeline → in-plant tank field, dock tank field, business storage tank field and land transportation crude oil pipeline, L (o, d) belongs to L.
L distill : in-plant → set of atmospheric and vacuum pipelines, L (o, d) E is L distill 。
L land : the land transportation crude oil pipeline set is formed by o belonging to N and d belonging to C cn ,l(o,d)∈L lanu ,
L mt-long : wharf tank field → pipeline set of long-distance pipeline, o belongs to C mt ,l(o,n long )∈L mt-long ,
L long-cn : long distance pipeline → pipeline set in tank field in plant, d ∈ C cn ,l(n long ,d)∈L long-cn ,
M e : oil mixing scheme set of atmospheric and vacuum equipment, wherein M belongs to M e ,e∈E。
Δ T : the length of time per time instant.
T desalt : standing for desalting time.
TL mt-cn : and calculating the oil head speed or event intervention at the oil head transmission moment of the long transmission pipeline.
V long : lower limit of total transmission capacity of long transmission pipeline.
tail mt : the wharf tank area descending proportion and the final tank area total tank storage descending coefficient upper limit.
tail cn : the in-plant tank area descending proportion and the final tank area total tank storage descending coefficient upper limit.
γ switch : the oil blending scheme is cost-effective for a single switch.
γ transit : and the cost of cutting the crude oil tank once.
γ entrance : single switch cost of pipeline oil.
Υ variation : the unit transmission of the pipeline fluctuates the cost.
(2) Decision variables
θ o,d,t,i : and a variable of 0-1, wherein the variable is 1 when the oil seeds i are conveyed from the tank area in the plant to the atmospheric and vacuum device pipelines l (o, d) at the moment t, and the variable is 0 otherwise.
λ o,d,t,i : the period t is the transmission quantity of the oil seeds i from the tank area in the plant to the pipelines l (o, d) of the atmospheric and vacuum distillation unit,
A variable 0-1, 1 means that the tank c is in an oil-chargeable state at the time t, otherwise 0, and/or>
A variable 0-1, 1 means that the tank c is in a stationary desalted state at time t, otherwise 0, and/or>
x e,t,m : the variable is 0 to 1, the period t is 1 when the atmospheric and vacuum device e adopts the mixing scheme m, otherwise, the period t is 0,
A variable 0-1, which is 1 when the entrance of the long transport pipeline l at the time t is switched, or is 0 when the entrance of the long transport pipeline l is not switched, or is combined with the long transport pipeline l at the time t>
A variable of 0 to 1 is 1 when the oil mixing scheme switching occurs in the atmospheric and vacuum device e at the moment t, otherwise, the variable is 0,
(3) Relaxation of the introduction
sc c,t : the continuous variable indicates the tank storage lower limit slack amount of the tank c at time t.
As shown in fig. 5, a crude oil scheduling optimization apparatus applied to an oil refinery includes:
the model establishing unit is used for establishing a crude oil scheduling optimization model based on a discrete time strategy;
the input unit is used for supporting flexible definition of model input parameters and manual intervention events based on the mathematical model established by the model establishing unit;
the model optimization unit is used for optimizing and solving the crude oil scheduling model according to the divide-and-conquer solving strategy to obtain a crude oil scheduling scheme and crude oil physical property values at all times in an event period;
the physical property checking unit is used for checking the physical property value of the crude oil obtained by the model optimizing unit;
and the output unit is used for forming a complete solution in the scheduling event period and outputting a report.
The embodiment of the invention also provides a crude oil scheduling optimization system applied to the oil refinery, which comprises:
the method comprises the following steps: one or more processors; a memory for storing one or more programs; the processor is configured to execute program instructions stored in the memory that when executed perform the above-described crude oil scheduling optimization method applied to a refinery.
Embodiments also provide a computer readable storage medium having a computer program stored thereon, wherein the computer program, when executed by one or more processors, implements the above-described crude oil scheduling optimization method applied to a refinery.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above described systems, systems and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
The algorithms and displays presented herein are not inherently related to any particular computer, virtual machine, or other apparatus. Various general purpose systems may also be used with the teachings herein. The required structure for constructing such a system will be apparent from the description above. Moreover, the present invention is not directed to any particular programming language. It is appreciated that a variety of programming languages may be used to implement the teachings of the present invention as described herein, and any descriptions of specific languages are provided above to disclose the best mode of the invention.
In addition, the memory may include volatile memory in a computer readable medium, random Access Memory (RAM) and/or nonvolatile memory such as Read Only Memory (ROM) or flash memory (flash RAM), and the memory includes at least one memory chip.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create a system for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including an instruction system which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
In a typical configuration, a computing device includes one or more processors (CPUs), input/output interfaces, network interfaces, and memory.
The memory may include forms of volatile memory in a computer readable medium, random Access Memory (RAM) and/or non-volatile memory, such as Read Only Memory (ROM) or flash memory (flash RAM). The memory is an example of a computer-readable medium.
Computer-readable media, including both non-transitory and non-transitory, removable and non-removable media, may implement information storage by any method or technology. The information may be computer readable instructions, data structures, modules of a program, or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static Random Access Memory (SRAM), dynamic Random Access Memory (DRAM), other types of Random Access Memory (RAM), read Only Memory (ROM), electrically Erasable Programmable Read Only Memory (EEPROM), flash memory or other memory technology, compact disc read only memory (CD-ROM), digital Versatile Discs (DVD) or other optical storage, magnetic cassettes, magnetic tape magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information that can be accessed by a computing device. As defined herein, a computer readable medium does not include a transitory computer readable medium such as a modulated data signal and a carrier wave.
It should also be noted that the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or apparatus. Without further limitation, an element defined by the phrases "comprising a," "8230," "8230," or "comprising" does not exclude the presence of additional identical elements in the process, method, article, or apparatus comprising the element.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
To further illustrate the technical solution of the present invention, the following 3 examples are provided, and the present invention is applicable to the topology structure of the dock tank field and the commercial tank field for crude oil transportation through the common pipeline, and the dock tank field for crude oil transportation through the long transportation pipeline to the in-plant tank field. The topological structures of the embodiment 1 and the embodiment 2 are shown in fig. 2, and the topological structures comprise the following components: wharf tank areas, commercial storage tank areas, in-plant tank areas, long-distance pipelines and oil storage tanks; the wharf tank area carries out one-way crude oil transmission to the in-plant tank area through the long transmission pipeline, the commercial storage tank area and the wharf tank area carry out two-way crude oil transmission through the common pipeline, and the oil head needs to be considered in the transmission of the long transmission pipeline. The following examples 1 and 2 are based on actual production data of the same refinery in the same time period, and only different settings are made for human intervention events. This refinery includes 1 oil ship, 3 tank districts, total 23 storage tanks, 2 sets of atmospheric and vacuum distillation unit CDU1 and CDU2, total scheduling duration is 7 days, the time scale is 2 hours, consider density and sulphur content as main rerum natura, oil ship arrival time is the scheduling start time, long distance pipeline all will incessant transport crude oil in whole cycle, atmospheric and vacuum distillation unit will process crude oil according to the oil mixing scheme in whole cycle, the artificial intervention incident can appoint crude oil to enter the jar, the emergence time of the incident of crude oil pipe transmission and jar oil type switching, the end time, the oil type, the oil mass, the tank number, the oil wheel is in batches etc.. The common parameters referred to above are:
basic information: the codes and names associated with the refinery resource allocation are shown in table 1.
TABLE 1
Setting parameters: the model calculates the settings of relevant parameters, including start time, end time, calculation scale (hours), number of desalination days, head stock (ton), blending schedule switch (yuan/time), tank field oil seed switch (yuan/time), crude oil relaxation cost (yuan/ton), tank storage relaxation cost coefficient, tanker storage cost coefficient, long transfer pipeline switch cost (yuan/time), pipeline rate fluctuation cost (yuan/(ton/hour)), business storage terminal mutual transfer cost (yuan/time), single time upper limit (seconds), terminal tank field maximum depreciation proportion, in-plant tank field maximum depreciation proportion, and head rate (ton/hour), as shown in table 2.
TABLE 2
Crude oil information: including crude oil code, crude oil name, type, physical property type, sulfur content, density, acid number, light heavy oil, marshalling, standing type, standing time, as shown in table 3.
TABLE 3
Crude oil sea transportation plan: refers to the plan for the tanker to deliver the crude oil, including the oil type, crude name, arrival dock, arrival time, estimated departure time, crude volume (ten thousand tons) and unit overdue cost (dollar/ton hour), as shown in table 4.
TABLE 4
Crude oil production plan: the throughput of each atmospheric and vacuum unit for each crude oil is shown in Table 5.
TABLE 5
Crude oil | Crude oil name | Production plan (ten thousand tons) | Atmospheric and vacuum equipment | Names of atmospheric and vacuum equipment |
RUS | Russian crude oil | 0.85 | CDU1 | Atmospheric and |
SUD | Nile crude oil | 0.00 | CDU1 | Atmospheric and |
RBB | Labe crude oil | 2.27 | CDU1 | Atmospheric and |
MEO | Merlo crude oil | 0.00 | CDU1 | Atmospheric and |
MIS | Milnas crude oil | 0.00 | CDU1 | Atmospheric and |
DAR | Darr crude oil | 0.00 | CDU2 | Atmospheric and vacuum 2 sets |
CBD | Crude carbene oil | 1.46 | CDU1 | Atmospheric and |
SAR | Salil crude oil | 0.00 | CDU1 | Atmospheric and |
DAQ | Daqing crude oil | 3.37 | CDU1 | Atmospheric and |
LIH | Liaohe crude oil | 3.98 | CDU2 | Atmospheric and vacuum 2 sets |
MNL | Maluna |
0 | CDU1 | Atmospheric and |
ESC | Eschlante crude oil | 0.60 | CDU2 | Atmospheric and vacuum 2 sets |
The oil mixing scheme comprises the following steps: the ratio of the atmospheric and vacuum equipment processing crude oil is shown in Table 6.
TABLE 6
Plan numbering | Atmospheric and vacuum equipment | Name of atmospheric and vacuum equipment | Crude oil | Crude oil name | Crude oil ratio | Beginning of the term plan | Remarks for |
1 | CDU1 | Atmospheric and |
RUS | Russian crude oil | 0.165653082 | Is that | |
1 | CDU1 | Atmospheric and |
CBD | Crude carbene oil | 0.174276119 | Is that | |
1 | CDU1 | Atmospheric and |
DAQ | Daqing crude oil | 0.6600708 | Is that | |
2 | CDU1 | Atmospheric and |
RUS | Russian crude oil | 0.094324264 | Whether or not | |
2 | CDU1 | Atmospheric and |
RBB | Labe crude oil | 0.344845164 | Whether or not | |
2 | CDU1 | Atmospheric and |
CBD | Crude carbene oil | 0.185038868 | Whether or not | |
2 | CDU1 | Atmospheric and |
DAQ | Daqing crude oil | 0.375791704 | Whether or not | |
6 | CDU2 | Atmospheric and vacuum 2 sets | LIH | Liaohe crude oil | 0.869346446 | Is that | |
6 | CDU2 | Atmospheric and vacuum 2 sets | ESC | Eschlante crude oil | 0.130653554 | Is that | |
7 | CDU2 | Atmospheric and vacuum 2 sets | DAR | Darr crude oil | 0.065 | Whether or not | |
7 | CDU2 | Atmospheric and vacuum 2 sets | LIH | Liaohe crude oil | 0.935 | Whether or not |
Tank area information: the information of crude oil stored in each storage tank comprises initial oil type, initial tank capacity, low limit alarm (safety lower limit), high limit alarm (safety upper limit), minimum tank quantity (physical lower limit), maximum tank quantity (physical upper limit), low limit relaxation, allowable switching, inventory cost (yuan/ten thousand ton hours), initial state, initial standing hours, payment while receiving, sulfur content (%), density (g/cm) 3 ). As shown in table 7.
TABLE 7
Pipeline crude oil plan: refers to the type of crude oil entering the plant via land lines, as shown in table 8.
TABLE 8
Starting point | Crude oil | Crude oil name |
Daqing | DAQ | Daqing crude oil |
Liaohe river | LIH | Liaohe crude oil |
Daqing Russian oil | RUS | Russian crude oil |
Pipeline information: refers to the topological structure and capacity limitations of the oil transport pipeline, as shown in table 9.
TABLE 9
Oil head information: refers to the section of crude oil that was initially present inside the long haul pipeline, as shown in table 10.
Watch 10
Example 1
In the embodiment 1, in the context of the embodiment, the parameter setting according to the actual demand scenario in the embodiment 1 includes the following steps:
TABLE 11
Module | Total number of variables | Continuous variable | Variable of integer | Variable 0-1 | Constraining |
Tanker-long transfer line entry | 142611 | 58466 | 84145 | 84143 | 588411 |
Long distance pipeline outlet-atmospheric pressure | 67307 | 29491 | 37816 | 37816 | 120150 |
And 2, setting manual intervention events and parameters according to actual scene requirements, and meeting the requirements of actual industrial scene application by setting the parameters and inputting the events. The parameter settings of example 1 share parameters with the examples, and the input event is set to null as shown in table 12.
TABLE 12
Type of event | Event description | Starting time | End time | Total amount of | Oil seed | Crude oil name | Starting point | Terminal point | Pot number | Batches of |
And 3, according to a divide-and-conquer solving strategy, carrying out division solving on the model in the step 1 under the parameter setting condition in the step 2, wherein the total process (input, modeling, divide-and-conquer solving and output) of the model shares the time of 227s.
And 4, checking the physical property values of the crude oil in the oil tank and the atmospheric and vacuum distillation unit in the step 3, wherein the physical properties mainly comprise density and sulfur content, and a trend tracking table (such as the oil tank) is shown in a table 13.
Watch 13
The density trend plots (e.g., CDU2 and H116) are shown in FIG. 6.
The trend plots for sulfur content (e.g., CDU1 and H111) are shown in FIG. 7.
And 5, checking the physical property in the step 4 to meet the requirement of a normal range, and performing report output, wherein the report output comprises the following steps:
scheduling Gantt chart: referring to the Gantt chart of operation of each resource (tank, pipeline, atmospheric and vacuum device), as shown in FIG. 8.
The material balance tables (sub-tank sections) for tank output and input are shown in tables 14 to 16.
TABLE 14
Watch 15
TABLE 16
The trend chart of the amount of stored oil in the tank is shown in table 17.
TABLE 17
Crude oil scheduling events summary table: refers to a schedule instruction table containing all the operational events, as shown in table 18.
Watch 18
A pipe delivery event table: the tubing events at the inlet and outlet of the long delivery line are indicated in table 19.
Watch 19
Watch 20
Crude oil processing report form: refer to the selected protocol and processed crude oil at each time of the atmospheric and vacuum unit as shown in table 21.
TABLE 21
Example 2:
in the embodiment 2, in the context of the embodiment, the parameter setting according to the actual demand scenario of the embodiment 2 includes the following steps:
TABLE 22
Module | Total number of variables | Continuous variable | Variable of integer | Variable 0-1 | Constraining |
Tanker-long transfer line entry | 142611 | 58466 | 84145 | 84143 | 588321 |
Long distance pipeline outlet-atmospheric pressure | 67307 | 29491 | 37816 | 37816 | 120150 |
And 2, setting manual intervention events and parameters according to actual scene requirements, and meeting the requirements of actual industrial scene application by setting the parameters and inputting the events. The parameter settings of example 2 are the same as the parameters of the example, the input events are the optimized output crude scheduling events of example 1, and 2 crude oil tank inlet events, 3 crude oil tank transfer events, 4 crude oil pipe inlet events, and 1 tank oil type switching event are specified, as shown in table 23.
TABLE 23
And 3, according to a divide-and-conquer solving strategy, carrying out division solving on the model in the step 1 under the parameter setting condition in the step 2, wherein the total process (input, modeling, divide-and-conquer solving and output) of the model shares the time of 159s.
And 4, checking the physical property values of the crude oil in the oil tank and the atmospheric and vacuum distillation unit in the step 3, wherein the physical properties mainly comprise density and sulfur content, and a trend tracking table (such as an oil tank) is shown in a table 24.
Watch 24
The density trend plots (e.g., CDU2 and H116) are shown in FIG. 10.
The trend plots for sulfur content (e.g., CDU1 and H111) are shown in FIG. 11.
And 5, checking the physical property in the step 4 to meet the requirement of a normal range, and performing report output, wherein the report output comprises the following steps:
scheduling Gantt chart: fig. 12 and 13 show operational gantt charts of various resources (tanks, pipelines, atmospheric and vacuum systems).
Tables for the balance of materials (sub-tank zones) for tank export and import are shown in tables 25 to 27.
TABLE 25
Watch 26
Watch 27
The trend chart of the amount of fuel stored in the tank is shown in table 28.
Watch 28
Crude oil scheduling events summary table: refers to a table of scheduling instructions that contains all the operational events, as shown in table 29.
Watch 29
A pipe delivery event table: the infusion events of the long infusion line inlet and outlet are shown in tables 30 and 31.
Watch 30
Watch 31
Crude oil processing report form: refer to the selected protocol and crude oil processed at each time of the atmospheric pressure unit as shown in table 32.
Watch 32
Example 3
In addition to the tank farm topologies shown in examples 1 and 2 above, the present invention is also applicable to other topologies where the wharf tank farm is carrying out crude oil transportation to the in-plant tank farm via a long transportation pipeline, such as the topology of the non-commercial tank farm shown in fig. 3, which consists of: wharf tank areas and in-plant tank areas, long-distance pipelines and oil storage tanks; the wharf tank field and the in-plant tank field carry out one-way crude oil transmission through long transmission pipelines, and the oil head also needs to be considered in the transmission of the long transmission pipelines.
A certain refinery comprises 1 oil tanker, 2 tank areas, 23 storage tanks in total, 2 sets of atmospheric and vacuum distillation units CDU1 and CDU2, the total scheduling time length is 7 days, the time scale is 2 hours, the density and the sulfur content are taken as main physical properties, the arrival time of the oil tanker is the scheduling starting time, a long-distance pipeline continuously transports crude oil in the whole period, the atmospheric and vacuum distillation units process the crude oil in the whole period according to an oil mixing scheme, an artificial intervention event can specify the crude oil to enter the tank, the occurrence time of the events of crude oil pipe transportation and tank oil type switching, the ending time, the oil type, the oil quantity, the tank number, the oil wheel batch and the like.
Based on the background, the method for setting the parameters of the scene meeting the actual requirements comprises the following steps:
Watch 33
Module | Total number of variables | Continuous variable | Variable of integer | Variable 0-1 | Constraining |
Tanker-long transfer line entry | 156253 | 63876 | 92377 | 92375 | 675659 |
Long distance pipeline outlet-atmospheric pressure | 67289 | 29482 | 37807 | 37807 | 120130 |
And 2, setting manual intervention events and parameters according to actual scene requirements, and meeting the requirements of actual industrial scene application by setting the parameters and inputting the events. Examples of input parameters are:
basic information: the codes and names associated with the refinery resource allocation are shown in table 34.
Watch 34
Setting parameters: the model calculates settings of relevant parameters including start time, end time, calculation scale (hours), number of desalination days, head stock (ton), blending schedule switch (yuan/time), tank field oil type switch (yuan/time), crude oil relaxation cost (yuan/ton), tank storage relaxation cost coefficient, tanker storage cost coefficient, long transfer pipeline switch cost (yuan/time), pipeline rate fluctuation cost (yuan/(ton/hour)), business storage terminal mutual transfer cost (yuan/time), single time upper limit (second), terminal tank field maximum stock-lowering proportion, in-plant tank field maximum stock-lowering proportion, and head rate (ton/hour), as shown in table 35.
Watch 35
Cost item | Numerical value | Remarks to note |
Starting time | 2020/7/2 0:00 | |
End time | 2020/7/9 0:00 | |
Calculating scale (hours) | 2 | |
Days of desalination | 5 | |
Oil head stock (ton) | 6500 | |
Oil mixing scheme switching (Yuan/Shi) | 1000 | |
Tank field oil seed switching (Yuan/Shi) | 10000 | |
Crude oil relaxation cost (Yuan/ton) | 10 | |
Tank slack cost factor | 5 | |
Tanker storage cost factor | 1.2 | |
Switching cost of long transmission pipeline (Yuan/Shi) | 1000 | |
Pipeline velocity fluctuation cost (Yuan/(ton/hr)) | 0.05 | |
Upper limit of single time (seconds) | 30 | |
Maximum warehouse descending proportion of wharf tank field | 0.1 | |
Maximum depreciation proportion of tank area in plant | 0.1 | |
Oil head rate (ton/hour) | 275 |
Crude oil information: including crude oil code, crude oil name, type, physical property type, sulfur content, density, acid number, light heavy oil, composition, standing type, standing time, as shown in table 36.
Watch 36
Marine crude oil transportation plan: refers to the plan for the tanker to deliver the crude oil, including oil type, crude name, arrival dock, arrival time, estimated departure time, crude volume (ten thousand tons) and unit cost per unit time out (dollar/ton hour), as shown in table 37.
Watch 37
Crude oil production plan: the production amounts of the respective atmospheric and vacuum apparatuses for the respective crude oils were as shown in Table 38.
Watch 38
The oil mixing scheme comprises the following steps: the ratio of the atmospheric and vacuum equipment to process crude oil is shown in Table 39.
Watch 39
Plan numbering | Atmospheric and vacuum equipment | Name of atmospheric and vacuum equipment | Crude oil | Crude oil name | Crude oil ratioExample (b) | Beginning of the term plan | Remarks to note |
1 | CDU1 | Atmospheric and |
RUS | Russian crude oil | 0.165653082 | Is that | |
1 | CDU1 | Atmospheric and |
CBD | Crude carbene | 0.174276119 | Is that | |
1 | CDU1 | Atmospheric and |
DAQ | Daqing crude oil | 0.6600708 | Is that | |
2 | CDU1 | Atmospheric and |
RUS | Russian crude oil | 0.094324264 | Whether or not | |
2 | CDU1 | Atmospheric and |
RBB | Labe crude oil | 0.344845164 | Whether or not | |
2 | CDU1 | Atmospheric and |
CBD | Crude carbene oil | 0.185038868 | Whether or not | |
2 | CDU1 | Atmospheric and |
DAQ | Daqing crude oil | 0.375791704 | Whether or not | |
6 | CDU2 | Atmospheric and vacuum 2 sets | LIH | Liaohe crude oil | 0.869346446 | Is that | |
6 | CDU2 | Atmospheric and vacuum 2 sets | ESC | Eschlanter crude | 0.130653554 | Is that | |
7 | CDU2 | Atmospheric and vacuum 2 sets | DAR | Darr crude oil | 0.065 | Whether or not | |
7 | CDU2 | Atmospheric and vacuum 2 sets | LIH | Liaohe crude oil | 0.935 | Whether or not |
Tank area information: the information of crude oil stored in each storage tank comprises initial oil type, initial tank capacity, low limit alarm (safety lower limit), high limit alarm (safety upper limit), minimum tank quantity (physical lower limit), maximum tank quantity (physical upper limit), low limit relaxation, allowable switching, inventory cost (yuan/ten thousand ton hours), initial state, initial standing hours, payment while receiving, sulfur content (%), density (g/cm) 3 ). As shown in table 40.
Watch 40
Pipeline crude oil plan: refers to the type of crude oil entering the plant via land transportation pipelines, as shown in table 41.
Table 41
Starting point | Crude oil | Crude oil name |
Daqing | DAQ | Daqing crude oil |
Liaohe river | LIH | Liaohe crude oil |
Daqing Russian oil | RUS | Russian crude oil |
Pipeline information: refers to the topology and capacity limitations of the oil transport lines, as shown in table 42.
Watch 42
Oil head information: refers to the section of crude oil that was initially present inside the long haul pipeline, as shown in table 43.
Watch 43
Inputting an event: the embodiment does not set human intervention and the input event is set to null as shown in table 44.
Watch 44
And 3, according to a divide-and-conquer solving strategy, carrying out division solving on the model in the step 1 under the parameter setting condition in the step 2, wherein the total process (input, modeling, divide-and-conquer solving and output) of the model shares 309s.
And 4, checking the physical property values of the crude oil in the oil tank and the atmospheric and vacuum distillation unit in the step 3, wherein the physical properties mainly comprise density and sulfur content, and a trend tracking table (such as an oil tank) is shown in a table 45.
TABLE 45
The density trend plots (e.g., CDU2 and H116) are shown in FIG. 14.
The sulfur content trend plots (e.g., CDU1 and H111) are shown in FIG. 15.
And 5, checking the physical property in the step 4 to meet the requirement of a normal range, and performing report output, wherein the report output comprises the following steps:
scheduling Gantt chart: referring to the Gantt diagram of the operation of each resource (tank, pipeline, atmospheric and vacuum), as shown in FIGS. 16 and 17.
Tank export and import material balance tables (sub-tank zones) as shown in tables 46 and 47.
TABLE 46
Watch 47
The trend chart of the amount of fuel stored in the tank is shown in table 48 below.
Watch 48
Crude oil scheduling events summary table: refers to the table of scheduling instructions that contains all the operational events, as shown in table 49.
Watch 49
A pipe delivery event table: the tubing events for long tubing line inlets and outlets are indicated in tables 50 and 51.
Watch 50
Watch 51
Crude oil processing report form: refer to the selected protocol and crude oil processed at each time of the atmospheric pressure unit, as shown in table 52.
Table 52
The foregoing is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, various modifications and decorations can be made without departing from the principle of the present invention, and these modifications and decorations should also be regarded as the protection scope of the present invention. Furthermore, although specific terms are employed herein, they are used in a generic and descriptive sense only and not for purposes of limitation.
Claims (11)
1. A crude oil scheduling optimization method applied to an oil refinery is characterized by comprising the following steps:
step (1), establishing a crude oil scheduling optimization model based on a discrete time strategy;
step (2), based on the mathematical model established in the step (1), flexibly defining model input parameters and manual intervention events;
step (3), according to a divide-and-conquer solving strategy, a crude oil scheduling model is optimized and solved to obtain a crude oil scheduling scheme and crude oil physical property values at all times in an event period;
and (4) checking the physical property values of the crude oil in the step (3): checking the physical properties of the crude oil stored in the tank at each moment and the physical properties of the atmospheric and vacuum processed crude oil, judging whether the checking is passed, if the checking is passed, executing the step (5), and otherwise, returning to the step (2) to check and adjust the input parameters and the manual intervention event;
and (5) forming a complete solution in the scheduling event period, and outputting a report, wherein the scheduling event summary table is used as the input of the next period scheduling to realize rolling scheduling.
2. The method of claim 1, wherein in step (1), the discrete-time-strategic crude scheduling optimization model comprises:
material balance model of crude oil transportation tanker: the oil unloading amount of the crude oil tanker in the period is equal to the carrying amount of the crude oil tanker in the period;
tank field storage tank material restraint model: the tank oil storage amount at the end of each time in the storage tank period is equal to the tank storage amount at the end of the last time plus the tank entering amount at the moment minus the tank exiting amount at the moment, the tank entering amount at the moment comprises the oil tanker oil unloading and tank entering amount and the other storage tank pipeline transmission and tank entering amount, and the tank exiting amount at the moment comprises the tank exiting amount transmitted to other storage tanks or an atmospheric and vacuum device by pipelines;
selecting a model for an atmospheric and vacuum device scheme: and selecting a processing scheme at each moment of the atmospheric and vacuum distillation unit, and processing the crude oil according to the processing scheme, wherein the total processing amount of the periodic oil seeds is more than or equal to the oil seed demand.
3. The method of claim 2, wherein the storage tank inventory is maintained within safe tank inventory upper and lower limits;
the pipeline transmission follows a topological structure, the transmission quantity is within the range of the upper limit and the lower limit of the pipeline transmission, and the transmission oil seeds are matched with the oil seeds of the starting tank and the target tank;
the oil head is considered in the long-distance pipeline transmission, and the oil type and the oil quantity at the inlet of the long-distance pipeline at the current moment are equal to the oil type and the oil quantity at the outlet of the long-distance pipeline at the corresponding moment after the current moment plus the transmission time of the long-distance pipeline;
the expression of the total cost obtained by adding the overdue cost of the oil tanker, the storage cost of the tank area, the switching cost of the atmospheric and vacuum oil mixing scheme, the switching cost of the oil type of the tank area and the fluctuation cost of the pipeline is as follows:
the total cost = tanker overdue cost + tank farm storage cost + atmospheric and vacuum oil mixing scheme switching cost + tank farm tank oil species switching cost + pipeline fluctuation cost.
4. The method according to claim 1, wherein in step (2), the parameters and the human intervention events are specifically:
the parameters include:
oil tanker input parameters: arrival time, estimated departure time, type and quantity of transported crude oil, arrival dock, unit overdue cost;
storage tank input parameters: the storage tank is located in a tank area, the upper limit and the lower limit of the capacity of a physical tank, the upper limit and the lower limit of the capacity of a safety tank, the current oil type and the stock, the initial state, the initial standing time, the stock cost, the physical property and whether the oil type is allowed to be switched or not;
pipeline input parameters: the pipeline transmission cost is the pipeline transmission cost, and the oil head information of the long-distance transmission pipeline is obtained;
the device inputs parameters: the crude oil proportioning is selected according to the initial scheme of the atmospheric and vacuum distillation unit;
the requirement input parameters are as follows: the periodic demand of each oil;
the human intervention event comprises: crude oil is fed into a tank, transferred to the tank, transported by a crude oil pipe and switched by the oil type of the tank; the human intervention event is derived from a scheduled event output scheduled in a last cycle.
5. The method of claim 1, wherein in step (3), the divide and conquer strategy comprises the following substeps:
step 31), module segmentation: the model is divided into three modules by taking a long transmission pipeline as a node: the long-distance pipeline comprises a long-distance pipeline front part, a long-distance pipeline and a long-distance pipeline rear part, wherein the long-distance pipeline front part comprises an oil tanker, a wharf tank area and a long-distance pipeline inlet, the long-distance pipeline comprises a long-distance pipeline inlet and a long-distance pipeline outlet, and the long-distance pipeline rear part comprises a long-distance pipeline outlet and an atmospheric and vacuum device;
step 32), solving before long-distance pipeline: adopting a strategy of integral modeling and solution emphasizing, neglecting the oil head transmission of the long-distance pipeline, assuming that the transmission time of the long-distance pipeline is fixed time, establishing an integral crude oil scheduling model, wherein an objective function comprises overall important cost, but the cost more emphasizes the minimum before the long-distance pipeline, and obtaining overall feasible solution and optimal solution/near-optimal solution before the long-distance pipeline;
step 33), solving the long-distance pipeline: simulating the inlet oil type and the oil quantity of the long-distance pipeline obtained in the step 32) according to the transmission rate obtained in the step 32) by adopting a simulation strategy and considering oil head transmission to obtain the outlet oil type and the oil quantity of the long-distance pipeline at each moment;
step 34), solving after long-distance pipeline transportation: modeling and solving are carried out only for the long-distance pipeline outlet-atmospheric and vacuum device by adopting a module modeling and solving strategy, and the cost of the target function is the minimum of the long-distance pipeline outlet-atmospheric and vacuum device, so that the optimal solution/near-optimal solution behind the long-distance pipeline is obtained;
step 35), the solution of the front part of the long transmission pipeline in the global variable in the step 32) and all the solutions after the long transmission pipeline in the step 34) are combined to form the solution of the whole crude oil scheduling.
6. The method of claim 5, wherein in step (3), the divide and conquer solution strategy further introduces slack, comprising: tank storage relaxation, crude oil production plan relaxation, and depalletizing relaxation.
7. The method of claim 1, wherein in step (5), the report comprises:
tank farm oil seed balance sheet: the inlet and outlet balance table of each oil type of each tank area in the whole period;
tank field storage tank balance table: the inlet and outlet balance tables of all storage tanks of all tank areas in the whole period;
tank trend table: storing oil seeds and oil mass in each storage tank at each moment;
scheduling event summary table: the scheduling event summary table comprises a crude oil tank inlet table, a crude oil pipe transmission table, a crude oil tank transfer table, an overland transportation table, a crude oil processing table, a tank oil type switching table and a physical property output table.
8. The method of claim 7, wherein the crude oil enters the tank: the method comprises the steps that the tanker transports crude oil to a wharf tank, and the oil tanker comprises a starting time, an ending time, an oil type, an oil quantity, a starting point and an end point;
the crude oil pipeline transportation: the wharf tank is used for conveying crude oil to an in-plant tank through a long conveying pipeline, and the in-plant tank comprises starting time, ending time, oil type, oil quantity, a starting point and a finishing point;
the crude oil transferring: the method refers to the mutual transmission between other crude oil storage tanks except for crude oil pipeline transmission, and comprises a starting time, an ending time, an oil type, an oil quantity, a starting point and an end point;
the land transportation comprises the following steps: the method comprises the steps that the ground transportation starting point conveys crude oil to an inner tank of a factory through a ground transportation pipeline, and the ground transportation starting point comprises starting time, ending time, oil type, oil quantity, a starting point and an end point;
crude oil processing: the method comprises the steps that crude oil is conveyed to an atmospheric and vacuum distillation unit by an in-plant tank for processing, and the method comprises an atmospheric and vacuum distillation unit processing scheme, processing amount, crude oil type and crude oil amount conveyed by the in-plant tank, starting time, ending time, starting point and ending point;
switching of the oil types of the tanks: the switching of the oil storage tank in the oil storage tank is realized, and the switching comprises switching time, oil seeds after switching and tank stock after switching;
physical property output table: comprises an atmospheric and vacuum distillation unit processing crude oil physical tracking list and a storage physical tracking list.
9. A crude oil scheduling optimization device applied to an oil refinery is characterized by comprising:
the model establishing unit is used for establishing a crude oil scheduling optimization model based on a discrete time strategy;
the input unit is used for supporting flexible definition of model input parameters and manual intervention events based on the mathematical model established by the model establishing unit;
the model optimization unit is used for optimizing and solving the crude oil scheduling model according to the divide-and-conquer solving strategy to obtain a crude oil scheduling scheme and crude oil physical property values at all times in an event period;
the physical property checking unit is used for checking the physical property value of the crude oil obtained by the model optimizing unit;
and the output unit is used for forming a complete solution in the scheduling event period and outputting a report.
10. A crude oil scheduling optimization system applied to an oil refinery is characterized by comprising:
the method comprises the following steps: one or more processors; a memory for storing one or more programs; the processor is configured to execute program instructions stored in the memory that when executed perform the method of crude oil schedule optimization as set forth in any one of claims 1 to 8 as applied to a refinery.
11. A computer readable storage medium, having a computer program stored thereon, wherein the computer program, when executed by one or more processors, implements the method of crude scheduling optimization as claimed in any one of claims 1 to 8 for use in a refinery.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111102934.2A CN115841165A (en) | 2021-09-18 | 2021-09-18 | Crude oil scheduling optimization method and device applied to oil refinery |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111102934.2A CN115841165A (en) | 2021-09-18 | 2021-09-18 | Crude oil scheduling optimization method and device applied to oil refinery |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115841165A true CN115841165A (en) | 2023-03-24 |
Family
ID=85574350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111102934.2A Pending CN115841165A (en) | 2021-09-18 | 2021-09-18 | Crude oil scheduling optimization method and device applied to oil refinery |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115841165A (en) |
-
2021
- 2021-09-18 CN CN202111102934.2A patent/CN115841165A/en active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Integrated production, inventory and delivery problems: Complexity and algorithms | |
Agra et al. | Combined ship routing and inventory management in the salmon farming industry | |
Prokudin et al. | Application of information technologies for the optimization of itinerary when delivering cargo by automobile transport | |
Irina et al. | Management of the Logistics Component of the Grain HarvestingProcess with Consideration of the Choice of Automobile Transport Technology Based on the Energetic Criterion | |
Cafaro et al. | Rigorous formulation for the scheduling of reversible-flow multiproduct pipelines | |
Kozachenko et al. | Creation of export-oriented network of grain elevators in Ukraine | |
CN114742502B (en) | Finished oil pipeline conveying coordination optimization system | |
Takeyasu et al. | Optimization technique by genetic algorithms for international logistics | |
CN115841165A (en) | Crude oil scheduling optimization method and device applied to oil refinery | |
Tsyganov | Mechanisms for managing large-scale transport and economic systems | |
EP4012871B1 (en) | Method for the directed transmission of energy in the form of energy packets | |
Zaynagabdinov et al. | Optimum planning use of equipment in agriculture | |
CN114091994A (en) | Data processing method and device for goods | |
WO2021151137A1 (en) | Method for optimizing transport parameters | |
Bamoumen et al. | A MILP Model for the Scheduling of a Multiproduct Tree-Structure Pipeline Network in Mining Industry | |
DE102012023579B4 (en) | Method for the needs-based supply of assembly and/or work stations in a synchronized assembly line production, especially in the vehicle industry | |
CN110490457A (en) | One kind being based on financial data acquisition process and monitoring system platform | |
Borodin et al. | An interval programming approach for an operational transportation planning problem | |
Panin et al. | Modeling the Optimal Mode of Operation of the Harvesting and Transport Links if Grain Is Stored in Flexible Polyethylene Containers | |
Kirillova et al. | Comparative characteristics of export transactions of a chemical industry enterprise | |
Shilovsky et al. | Algorithm for the Development and Delivery of a Multi-Item Set of Spare Parts and Maintenance Supplies for Geographically Dispersed Consumers | |
US11636415B2 (en) | Systems and methods for supply management | |
CN110825045B (en) | Method for improving switching efficiency of multi-device multi-product processing schemes | |
Wang et al. | Liner Ship Fleet Planning: Models and Algorithms | |
Mostafaei et al. | Hierarchical Decomposition Approach for Detailed Scheduling of Pipeline Systems with Branches |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication |