CN112766862A - Method for dividing and locating distribution area of petrochemical product supply chain - Google Patents
Method for dividing and locating distribution area of petrochemical product supply chain Download PDFInfo
- Publication number
- CN112766862A CN112766862A CN202110157188.0A CN202110157188A CN112766862A CN 112766862 A CN112766862 A CN 112766862A CN 202110157188 A CN202110157188 A CN 202110157188A CN 112766862 A CN112766862 A CN 112766862A
- Authority
- CN
- China
- Prior art keywords
- distribution
- oil depot
- towns
- oil
- radius
- 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
- 238000009826 distribution Methods 0.000 title claims abstract description 185
- 238000000034 method Methods 0.000 title claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 7
- 241001061225 Arcos Species 0.000 claims description 6
- 239000003348 petrochemical agent Substances 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 5
- 238000012790 confirmation Methods 0.000 claims description 3
- 239000003921 oil Substances 0.000 description 115
- 238000005457 optimization Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000007774 longterm Effects 0.000 description 3
- 239000003208 petroleum Substances 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 2
- 208000001848 dysentery Diseases 0.000 description 2
- 239000000446 fuel Substances 0.000 description 2
- 239000007789 gas Substances 0.000 description 2
- 239000003502 gasoline Substances 0.000 description 2
- 239000000314 lubricant Substances 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 206010039203 Road traffic accident Diseases 0.000 description 1
- 239000010426 asphalt Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 239000003225 biodiesel Substances 0.000 description 1
- 229910002092 carbon dioxide Inorganic materials 0.000 description 1
- 239000001569 carbon dioxide Substances 0.000 description 1
- 239000013064 chemical raw material Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- -1 diesel Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000295 fuel oil Substances 0.000 description 1
- 239000003350 kerosene Substances 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000012188 paraffin wax Substances 0.000 description 1
- 239000002006 petroleum coke Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 239000002904 solvent Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/08—Logistics, e.g. warehousing, loading or distribution; Inventory or stock management
- G06Q10/083—Shipping
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q30/00—Commerce
- G06Q30/02—Marketing; Price estimation or determination; Fundraising
- G06Q30/0201—Market modelling; Market analysis; Collecting market data
- G06Q30/0204—Market segmentation
- G06Q30/0205—Location or geographical consideration
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- Economics (AREA)
- Finance (AREA)
- Accounting & Taxation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Data Mining & Analysis (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for dividing and selecting addresses of distribution areas of a petrochemical product supply chain, which comprises the following steps: listing sequence codes and coordinate information of all towns in the selected area, determining distribution radius ranges of distribution areas, comparing and selecting the most suitable distribution radius, dividing each distribution area, carrying out first adjustment on the position of the oil depot, carrying out second adjustment on the position of the oil depot, carrying out third adjustment on the position of the oil depot, and carrying out fourth adjustment on the position of the oil depot; the method comprises the steps of firstly listing sequence codes and coordinate information of all towns in a selected area, determining distribution radius ranges of distribution areas through calculation, then comparing and selecting the most suitable distribution radius, dividing each distribution area, firstly adjusting the position of an oil depot for the first time, then adjusting the position of the oil depot for the second time, adjusting the position of the oil depot for the third time, and finally adjusting the position of the oil depot for the fourth time, so that the transportation cost can be greatly reduced.
Description
Technical Field
The invention relates to the technical field of petrochemical product distribution, in particular to a method for dividing and locating distribution areas of a petrochemical product supply chain.
Background
Petrochemical products are petrochemical products for short, that is, products which can be produced by petrochemical industry, and the petrochemical products can be divided into: petroleum fuel, petroleum solvent and chemical raw material, lubricant, paraffin, petroleum asphalt, petroleum coke and other 6 types. Among them, the production of each fuel is the largest, accounting for about 90% of the total production; the variety of various lubricants is the most, and the yield accounts for about 5 percent. The product standards are established by all countries to meet the requirements of production and use.
Oil refineries, oil depots and towns are considered as three major entities in the petrochemical supply chain. Typically, petrochemicals need to be transported from refineries, through oil depots, to towns. Generally, a large area is divided into a plurality of distribution areas for storing and distributing petrochemical products. While a large area may be divided into "area with refinery" and "area without refinery", the two areas differ in dividing distribution areas: the first distribution sector division of "a sector with a refinery" needs to be centered around the refinery, while the "sector without a refinery" needs to consider only the maximum averaging of the distributions in the sectors at the time of the sector division.
The problems of unreasonable site selection of an oil depot, high occupied transportation road, long delivery time, low oil depot occupancy rate and the like mainly exist in the current gasoline transportation. The oil depot site selection is a node for logistics transportation optimization, and the current discussion of China on petrochemical product transportation and distribution optimization is basically limited to theoretical mode discussion, so that problems are rarely solved by establishing a mathematical model and are verified by simulation. At present, the supply chain optimization work of petrochemical products in China is only optimized from the aspects of purchase, configuration, production and processing, and the scheme of optimizing from the aspect of oil depot site selection is in short supply.
Based on the above, the present invention provides a method for dividing and locating distribution areas of a petrochemical supply chain, so as to solve the above-mentioned problems.
Disclosure of Invention
The invention aims to provide a method for dividing and selecting sites of petrochemical product supply chain distribution areas, which can greatly reduce the transportation cost by firstly listing sequence codes and coordinate information of all towns in a selected area, determining the distribution radius range of the distribution areas through calculation, then comparing and selecting the most suitable distribution radius, dividing each distribution area, firstly adjusting the position of an oil depot for the first time, then adjusting the position of the oil depot for the second time, adjusting the position of the oil depot for the third time, and finally adjusting the position of the oil depot for the fourth time, so as to solve the problems in the background technology.
In order to achieve the purpose, the invention provides the following technical scheme: a method for dividing and locating distribution areas of a petrochemical product supply chain comprises the following steps:
s1, listing sequence codes and coordinate information of all towns in the selected region;
s2, determining the distribution radius range of the distribution area: the method comprises the steps of adding a calculation formula of a spherical area to a research area, calculating the distance between all towns including an oil depot and an oil refinery by a latitude and longitude matrix calculation method and using MATLAB, and determining an effective distribution range by the method, wherein the method comprises the following specific steps:
minimum delivery range value: means that the smallest distribution distance from any town to any other town in the area under study, that is, the distances from all towns in the area to other towns are calculated and compared to obtain the shortest transportation distance, which is the smallest distribution radius (the smallest distribution range value) because any town can at least reach another town to consider the radius (range) as meaningful.
Maximum delivery range value: the maximum value of the distribution range in the research area is that any town can include all other towns, the distribution range can be infinite (because the larger the distribution range is, all towns can be included more easily, but the distribution range is not the optimal upper limit range), but the program only selects the minimum distribution range which can include all other towns as the maximum distribution range value, in other words, the maximum distribution range value is the calculation upper limit of the program for calculating the optimal distribution area dividing radius;
s3, comparing and selecting the most suitable dispensing radius: selecting a distribution radius within the range of the minimum effective distribution radius and the maximum effective distribution radius, wherein the distance from the distribution radius to any town is shortest, and the radius is used as the optimal radius to carry out the first division on the selected region;
s4, dividing each distribution area into: through the confirmation of the optimal delivery radius, the various delivery areas will be divided by the MATLAB program. Each distribution area has at least one oil depot, and theoretically, the distribution distance from a central point (a temporary central point of the oil depot position) to each town in each distribution area is less than or equal to 15 kilometers. The oil refinery has its own oil depot which will be used as the oil depot of the first divided region to distribute petrochemical products to towns within 15 km, after the first distribution range is confirmed, the distribution region will not be used as the candidate point any more and will disappear from the figure, and the rest towns will be divided step by 15 km of distribution radius. The division conditions are processed from more to less in sequence. After the division is finished, a few remote towns are not contained in the divided distribution areas;
s5, adjusting the position of the oil depot for the first time: after distribution areas are divided, the oil depot position is firstly pushed from the distribution radius (15 kilometers) central point position of each distribution area by taking the demand of different towns on petrochemical products as weight. At this time, the oil depot positions in the distribution areas are relatively optimal, but relatively distant towns which exceed the optimal distribution radius by 15 kilometers are not distributed, and all the relatively distant towns which exceed the optimal distribution radius are considered in the next step;
s6, adjusting the position of the oil depot for the second time: on the basis of the positions of the oil depots in S5, all towns exceeding the optimal distribution radius are calculated and compared by embedding MATLAB into a hemiversine formula matrix, the towns are classified into the oil depots which are close to the towns, the demands and the weights of the distribution areas for petrochemical products change due to the addition of the towns, the positions of the oil depots which are newly added into the town distribution areas can be secondarily shifted according to the weights (the demands of the petrochemical products) of the towns in the areas through recalculation, at the moment, the optimal distribution distances of the towns at the boundaries of the distribution areas also change due to the change of the positions of the oil depots, and only the nearest distribution distances are optimal from the oil depots to the towns through comparison, namely, the distribution distances from the oil depots to the towns are confirmed through matrix calculation and comparison again, and the distribution distances from the oil depots to the towns are compared, relatively short distribution distances are combined into respective matrix sets, and the latest distribution areas are confirmed;
and S7, adjusting the position of the oil depot for the third time: since the weight of the primary stream from the oil refinery to the oil depot is much larger than that of the secondary stream from the oil depot to the town, and the transportation modes mainly used in the transportation of the primary stream are ship transportation and railway transportation, the nearest port (or along the river) and the railway platform (or along the railway) are the final positions of the oil depots in each distribution area based on the optimized positions of the oil depots in S6;
and S8, adjusting the position of the oil depot for the fourth time: and finally, dividing the distribution areas into different distribution areas, namely calculating and comparing points of all town positions in the matrix and points of the final oil depot positions in the matrix by using a hemiversive formula based on the final oil depot positions along river routes or railway lines of the distribution areas in S7 through matrix calculation and comparison in a program, finally dividing the distribution areas with short distribution distances from all towns to all oil depots as the final distribution areas, and performing cluster division in the form of the matrix, wherein finally, the different distribution areas have different distribution radiuses, and the oil depot positions are optimized according to the characteristics of the different distribution areas.
Preferably, at least two digits after the decimal point are reserved for the coordinates of all towns, and in the operation, the more digits reserved after the decimal point are, the larger the calculation amount is, and the more accurate the calculation amount is.
Preferably, the dividing rule in S3 is: the shortest possible delivery radius is used to account for as many towns as possible, and the number of towns in each divided delivery area is as average as possible.
Preferably, the calculation formula in S2 is:
the limiting conditions are as follows:
xij∈{0,1},i＝1,2,3,...,n；j＝1,2,3,...,m
parameters are as follows:
i: town index (i ═ 1,2,3, …, n);
j: refinery index (j ═ 1,2,3, …, m);
j: the number of distribution areas to be tested;
RA: refinery coordinates;
TL: oil depot coordinates;
TLj: the coordinates of the oil depot to be tested;
SG: coordinates of each town;
n: total number of towns;
m: total number of oil depots;
and Vi: total consumption of petrochemical products in town i;
nj: demand for petrochemicals in town j;
TKp: the transportation cost of primary logistics from an oil refinery to an oil depot;
TKs: the transportation cost of secondary logistics from an oil depot to a town;
RA ═ e, f: latitude and longitude of the refinery;
TLj ═ (cj, dj): latitude and longitude of the oil depot;
SGi ═ (ai, bi): latitude and longitude of the town;
f (RA, TLj): primary stream transport distance from refinery to depot j, from (e, f) to (cj, dj);
f(RA,TLj)＝R*arcos[sin(e)*sin(cj)+cos(e)*cos(cj)*cos(f-dj)]*Pi/180
r is the radius of the earth;
g (TLj, SGi): the secondary stream transport distance from depot j to town i, from (cj, dj) to (ai, bi);
g(TLj,SGi)＝R*arcos[sin(cj)*sin(ai)+cos(cj)*cos(ai)*cos(dj-bi)]*Pi/180
r is the radius of the earth.
Preferably, xijWhich may be a binary decision variable, which may be "1" or "0", which varies depending on whether town i is classified into a distribution area, by which judgment the shortest transportation distance can be determined.
Preferably, the hemiversine formula in S6 is as follows:
d＝R·c
Compared with the prior art, the invention has the beneficial effects that:
1. the invention can greatly improve the transport efficiency of the petrochemical supply chain and reduce the transport cost by comparing different influence factors and adjusting the distribution area of the petrochemical products and the position of the oil depot.
2. In the invention, the sequence codes and the coordinate information of all towns in the selected area are listed firstly, the distribution radius range of the distribution area is determined by calculation, then the most suitable distribution radius is compared and selected, each distribution area is divided, the position of the oil depot is firstly adjusted for the first time, then the position of the oil depot is adjusted for the second time and the position of the oil depot is adjusted for the third time, and finally the position of the oil depot is adjusted for the fourth time, so that the transportation cost can be greatly reduced.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic diagram illustrating the division of petrochemical distribution areas and the adjustment of the location of a tank farm according to the present invention;
FIG. 2 is a schematic diagram showing the shortest time for loading the first logistics goods before and after the optimization according to the present invention;
FIG. 3 is a schematic diagram showing the comparison between the number of tankers required for the first logistics and the amount of petrochemical cargo required for loading according to the present invention;
FIG. 4 is a schematic diagram showing the comparison of the occupancy rates of primary logistics relative to the river before and after the optimization of the present invention;
FIG. 5 is a schematic diagram showing the comparative ratio of the occupation ratio of the secondary logistics roads before and after the structural optimization of the present invention;
FIG. 6 is a statistical comparison of transportation costs and other costs of the structure of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example 1
Referring to fig. 1, the present invention provides a technical solution for dividing and locating distribution areas of a petrochemical supply chain: the method comprises the following steps:
s1, listing sequence codes and coordinate information of all towns in the selected region;
s2, (please refer to FIG. 1: a, b) determining the distribution radius range of the distribution area: the method comprises the steps of adding a calculation formula of a spherical area to a research area, calculating the distance between all towns including an oil depot and an oil refinery by a latitude and longitude matrix calculation method and using MATLAB, and determining an effective distribution range by the method, wherein the method comprises the following specific steps:
minimum delivery range value: means that the smallest distribution distance from any town to any other town in the area under study, that is, the distances from all towns in the area to other towns are calculated and compared to obtain the shortest transportation distance, which is the smallest distribution radius (the smallest distribution range value) because any town can at least reach another town to consider the radius (range) as meaningful.
Maximum delivery range value: the maximum value of the distribution range in the research area is that any town can include all other towns, the distribution range can be infinite (because the larger the distribution range is, all towns can be included more easily, but the distribution range is not the optimal upper limit range), but the program only selects the minimum distribution range which can include all other towns as the maximum distribution range value, in other words, the maximum distribution range value is the calculation upper limit of the program for calculating the optimal distribution area dividing radius;
s3, (please refer to FIG. 1: a, b) compares and selects the most suitable dispensing radius: selecting a distribution radius within the range of the minimum effective distribution radius and the maximum effective distribution radius, wherein the distance from the distribution radius to any town is shortest, and the radius is used as the optimal radius to carry out the first division on the selected region;
s4, (please refer to FIG. 1: c-g) dividing each distribution area: through the confirmation of the optimal delivery radius, the various delivery areas will be divided by the MATLAB program. Each distribution area has at least one oil depot, and theoretically, the distribution distance from a central point (a temporary central point of the oil depot position) to each town in each distribution area is less than or equal to 15 kilometers. The oil refinery has its own oil depot which will be used as the first partitioned oil depot to distribute petrochemical products to towns within 15 km, as can be seen from fig. 1d, after the first distribution range is confirmed, the distribution area will no longer be the candidate point and will disappear from the drawing, and the rest towns will be gradually partitioned by 15 km of distribution radius. The division conditions are processed from more to less in sequence. After the division is completed, a few more remote towns are not included in the divided distribution area, such as town "m" in fig. 1;
s5, the position of the oil depot is adjusted for the first time (please refer to the figure 1: g-h): after the distribution areas are divided, the location of the oil depot is shifted from the location of the center of the distribution radius (15 km) of each distribution area by using the demand for petrochemicals in different towns as a weight (see h in fig. 1). At this time, the oil depot positions in the distribution areas are relatively optimal, but relatively distant towns which exceed the optimal distribution radius by 15 kilometers are not distributed, and all the relatively distant towns which exceed the optimal distribution radius are considered in the next step;
s6, adjusting the position of the oil depot for the second time (please refer to the figure 1: h-i): on the basis of the positions of the oil depots in S5, all towns exceeding the optimal distribution radius are calculated and compared by embedding MATLAB into a hemiversine formula matrix, the towns are classified into the oil depots which are close to the towns, the demands and the weights of the distribution areas for petrochemical products change due to the addition of the towns, the positions of the oil depots which are newly added into the town distribution areas can be secondarily shifted according to the weights (the demands of the petrochemical products) of the towns in the areas through recalculation, at the moment, the optimal distribution distances of the towns at the boundaries of the distribution areas also change due to the change of the positions of the oil depots, and only the nearest distribution distances are optimal from the oil depots to the towns through comparison, namely, the distribution distances from the oil depots to the towns are confirmed through matrix calculation and comparison again, and the distribution distances from the oil depots to the towns are compared, relatively short distribution distances are combined into respective matrix sets, and the latest distribution areas are confirmed;
s7, the position of the oil depot is regulated for the third time (please refer to the figure 1: i-j): since the weight of the primary stream from the oil refinery to the oil depot is much larger than that of the secondary stream from the oil depot to the town, and the transportation modes mainly used in the transportation of the primary stream are ship transportation and railway transportation, the nearest port (or along the river) and the railway platform (or along the railway) are the final positions of the oil depots in each distribution area based on the optimized positions of the oil depots in S6;
s8, (please refer to FIG. 1: j) adjusting the position of the oil depot for the fourth time: and finally, dividing the distribution areas into different distribution areas, namely calculating and comparing points of all town positions in the matrix and points of the final oil depot positions in the matrix by using a hemiversive formula based on the final oil depot positions along river routes or railway lines of the distribution areas in S7 through matrix calculation and comparison in a program, finally dividing the distribution areas with short distribution distances from all towns to all oil depots as the final distribution areas, and performing cluster division in the form of the matrix, wherein finally, the different distribution areas have different distribution radiuses, and the oil depot positions are optimized according to the characteristics of the different distribution areas.
At least two digits after the decimal point are required to be reserved for the coordinates of all towns, and in the operation, the more digits reserved after the decimal point are, the larger the calculation amount is, and the more accurate the calculation amount is.
Wherein the dividing principle in S3 is as follows: the shortest possible delivery radius is used to account for as many towns as possible, and the number of towns in each divided delivery area is as average as possible.
Wherein the calculation formula in S2 is:
the limiting conditions are as follows:
xij∈{0,1},i＝1,2,3,...,n；j＝1,2,3,...,m
parameters are as follows:
i: town index (i ═ 1,2,3, …, n);
j: refinery index (j ═ 1,2,3, …, m);
j: the number of distribution areas to be tested;
RA: refinery coordinates;
TL: oil depot coordinates;
TLj: the coordinates of the oil depot to be tested;
SG: coordinates of each town;
n: total number of towns;
m: total number of oil depots;
and Vi: total consumption of petrochemical products in town i;
nj: demand for petrochemicals in town j;
TKp: the transportation cost of primary logistics from an oil refinery to an oil depot;
TKs: the transportation cost of secondary logistics from an oil depot to a town;
RA ═ e, f: latitude and longitude of the refinery;
TLj ═ (cj, dj): latitude and longitude of the oil depot;
SGi ═ (ai, bi): latitude and longitude of the town;
f (RA, TLj): primary stream transport distance from refinery to depot j, from (e, f) to (cj, dj);
f(RA,TLj)＝R*arcos[sin(e)*sin(cj)+cos(e)*cos(cj)*cos(f-dj)]*Pi/180
r is the radius of the earth;
g (TLj, SGi): the secondary stream transport distance from depot j to town i, from (cj, dj) to (ai, bi);
g(TLj,SGi)＝R*arcos[sin(cj)*sin(ai)+cos(cj)*cos(ai)*cos(dj-bi)]*Pi/180
r is the radius of the earth.
Wherein x isijWhich may be a binary decision variable, which may be "1" or "0", which varies depending on whether town i is classified into a distribution area, by which judgment the shortest transportation distance can be determined.
Wherein, the hemiversine formula in S6 is:
d＝R.c
Example 2
Taking optimization of petrochemical product supply chain in Germany Ruhr region as an example, optimizing results through simulation contrast
a. Loading time of cargo in primary stream (from refinery to oil depot)
Due to the large volume (including volume) of petrochemical products, the location of the loading time in the transport is very important, especially when the primary distribution of the petrochemical products has a great impact on the overall supply chain. In addition, the idling (return trip) of the tanker is also a very important factor, and when the tanker is delivered, it needs to return to the refinery in time to prepare for the next delivery, which means that the number of tankers is very important in the whole logistics.
It can be seen from the results of fig. 2 that the port requires 29 tankers (the shortest loading time including the return empty) for distribution before and after optimization, but the optimized loading time is saved by about 4.28 hours (138,477 seconds-123,079 seconds) compared with that before optimization.
b. The amount of goods to be transported in the primary logistics
The transportation of petrochemicals is a bulk cargo transportation. In general, water transportation and railway transportation are the main transportation means thereof, by which a large amount of goods can be efficiently transported within a short transportation time. In order to minimize the pressure in the reservoir, it is important to have a quantitative hold during the dispensing process.
As can be seen from fig. 3, when the primary distribution of petrochemical products is carried out by the tanker, the cargo transportation amount in the luer region is obviously reduced after the distribution region range and the oil depot position are adjusted. The main reason for the reduction of the freight volume in the luer district is that the optimized distribution area is divided according to the oil refinery, and all towns near the oil refinery are directly distributed by the oil depot of the oil refinery, that is, in the distribution area of the oil refinery, the primary stream and the secondary stream are integrated, so that the freight distribution volume of towns near the oil refinery can be reduced. The overall reduction in cargo transport is about 12.9%.
c. Relative occupancy rate of primary logistics water path
A successful logistics system must be achieved with the shortest latency and most efficient return time. Excessive transportation costs and transportation costs need to be avoided as much as possible by accurate delivery times, and only so much as to increase the competitive advantage of petrochemical transportation.
Through comparison of results in fig. 4, it can be found that the river occupancy rate of primary logistics transportation is obviously reduced after optimization. Furthermore, it can be seen from the comparison of the graphs that all the transportation time ends in advance after optimization, that is, all the oil ships complete the distribution in advance (including the idle return), which means that the tasks of the total transportation amount can be completed with less river occupancy and shorter distribution time through the repartition of the distribution area and the adjustment of the oil depot position. The whole efficiency is greatly improved.
d. Relative road occupancy rate of secondary logistics (from oil depot to gas station)
Transportation from the oil depot to the gas stations in various towns via highways, provincial roads and common streets is the main transportation mode for secondary logistics of petrochemical products. The relative occupancy of the haul road is a very important assessment index. In the model, the occupation of all roads is simulated through the roads in the model, and specific data in the transportation process are obtained. The following data results were obtained by model operation.
e. Capacity load of oil depot
The reserve of the oil depot is a very important index for judging whether the oil depot can normally operate or not. Long term high reserves can cause a number of cost and cost problems, and in addition, long term high reserves can cause parts of the tank to become inoperative, for example, certain specialty oil products are not suitable for long term storage, such as summer diesel and winter diesel, and seasonal demand is not suitable for use and storage in the opposite season. Therefore, in the secondary stream, the reserve of the reservoir is very important.
The Duisburg oil depot and the Hu nxe oil depot are two important transport oil depots, and the oil storage capacity of the two important transport oil depots is very large. The storable oil product is: gasoline, diesel, kerosene, heavy oil, natural gas, biodiesel, and the like (total reservoir reserves in the luer region of about 20 hundred million liters). It can be found by simulation that after the global optimization, the luer region decreased by 1% in total (about 200 kiloliters) from the first week to the end of the fourth week of simulation compared to the existing reservoir storage state. Through the decline of oil depot reserves, can make the whole operation of oil depot more high-efficient, also can make the use of each oil storage tank of oil depot more nimble simultaneously.
f. Statistical data of transportation cost
Finally, in order to further compare statistical data after the optimization of the oil depot position, the article respectively counts the primary logistics through a waterway, a railway and a highway and calculates according to the highway by matching with the secondary logistics, and adds other transportation costs (including noise treatment cost, carbon dioxide treatment cost, other harmful substance treatment cost and traffic accident treatment cost) on the basis of the basic transportation cost, and detailed comparison statistics is carried out on the transportation of the petrochemical products in the whole Ruhr region. By comparison, it can be seen that the area has a total transportation cost reduction of about 10.16% after optimization according to the optimal transportation method
In the description herein, references to the description of "one embodiment," "an example," "a specific example" or the like are intended to mean that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above do not necessarily refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
The preferred embodiments of the invention disclosed above are intended to be illustrative only. The preferred embodiments are not intended to be exhaustive or to limit the invention to the precise embodiments disclosed. Obviously, many modifications and variations are possible in light of the above teaching. The embodiments were chosen and described in order to best explain the principles of the invention and the practical application, to thereby enable others skilled in the art to best utilize the invention. The invention is limited only by the claims and their full scope and equivalents.
Claims (6)
1. A method for dividing and locating distribution areas of a petrochemical product supply chain is characterized by comprising the following steps:
s1, listing sequence codes and coordinate information of all towns in the selected region;
s2, determining the distribution radius range of the distribution area: the method comprises the steps of adding a calculation formula of a spherical area to a research area, calculating the distance between all towns including an oil depot and an oil refinery by a calculation method of a latitude and longitude matrix and using MATLAB;
s3, comparing and selecting the most suitable dispensing radius: selecting a distribution radius within the range of the minimum effective distribution radius and the maximum effective distribution radius, wherein the distance from the distribution radius to any town is shortest, and the radius is used as the optimal radius to carry out the first division on the selected region;
s4, dividing each distribution area into: through the confirmation of the optimal delivery radius, the various delivery areas will be divided by the MATLAB program. Each distribution area is provided with at least one oil depot;
s5, adjusting the position of the oil depot for the first time: after distribution areas are divided, the oil depot position is firstly pushed from the distribution radius central point position of each distribution area according to the requirements of different towns on petrochemical products as weights;
s6, adjusting the position of the oil depot for the second time: on the basis of the oil depot position in S5, by the calculation of embedding MATLAB into a hemiversine formula matrix, all towns exceeding the optimal distribution radius are calculated and compared, and are classified into oil depot distribution areas close to the towns, because the towns are newly classified into the towns, the demands and weights of the distribution areas for petrochemical products are changed, and by the calculation again, the oil depot position newly added into the town distribution areas can be secondarily shifted according to the weights of the towns in the areas, and at the moment, due to the change of the oil depot position, the optimal distribution distance of the boundary towns of the distribution areas can also be changed;
and S7, adjusting the position of the oil depot for the third time: since the weight of the primary stream from the oil refinery to the oil depot is much greater than that of the secondary stream from the oil depot to the town, and the transportation modes mainly used in the transportation of the primary stream are ship transportation and railway transportation, the nearest port and railway platform are the final positions of the oil depots in each distribution area based on the optimized position of the oil depot in S6;
and S8, adjusting the position of the oil depot for the fourth time: and finally, dividing the distribution areas into different distribution areas, namely calculating and comparing points of all town positions in the matrix and points of the final oil depot positions in the matrix by using a hemiversive formula based on the final oil depot positions along river routes or railway lines of the distribution areas in S7 through matrix calculation and comparison in a program, finally dividing the distribution areas with short distribution distances from all towns to all oil depots as the final distribution areas, and performing cluster division in the form of the matrix, wherein finally, the different distribution areas have different distribution radiuses, and the oil depot positions are optimized according to the characteristics of the different distribution areas.
2. The method of claim 1, wherein the petrochemical supply chain distribution area comprises: all town coordinates need to retain at least two digits after the decimal point.
3. The method of claim 1, wherein the petrochemical supply chain distribution area comprises: the dividing principle in S3 is as follows: the shortest possible delivery radius is used to account for as many towns as possible, and the number of towns in each divided delivery area is as average as possible.
4. The method of claim 1, wherein the petrochemical supply chain distribution area comprises: the calculation formula in S2 is:
the limiting conditions are as follows:
xij∈{0,1},i＝1,2,3,...,n；j＝1,2,3,...,m
parameters are as follows:
i: town index (i ═ 1,2,3, …, n);
j: refinery index (j ═ 1,2,3, …, m);
j: the number of distribution areas to be tested;
RA: refinery coordinates;
TL: oil depot coordinates;
TLj: the coordinates of the oil depot to be tested;
SG: coordinates of each town;
n: total number of towns;
m: total number of oil depots;
and Vi: total consumption of petrochemical products in town i;
nj: demand for petrochemicals in town j;
TKp: the transportation cost of primary logistics from an oil refinery to an oil depot;
TKs: the transportation cost of secondary logistics from an oil depot to a town;
RA ═ e, f: latitude and longitude of the refinery;
TLj ═ (cj, dj): latitude and longitude of the oil depot;
SGi ═ (ai, bi): latitude and longitude of the town;
f (RA, TLj): primary stream transport distance from refinery to depot j, from (e, f) to (cj, dj);
f(RA,TLj)＝R*arcos[sin(e)*sin(cj)+cos(e)*cos(cj)*cos(f-dj)]*Pi/180
r is the radius of the earth;
g (TLj, SGi): the secondary stream transport distance from depot j to town i, from (cj, dj) to (ai, bi);
g(TLj,SGi)＝R*arcos[sin(cj)*sin(ai)+cos(cj)*cos(ai)*cos(dj-bi)]*Pi/180
r is the radius of the earth.
5. The method of claim 4, wherein the petrochemical supply chain distribution area comprises: x is the number ofijWhich may be a binary decision variable, which may be "1" or "0", which varies depending on whether town i is classified into a distribution area, by which judgment the shortest transportation distance can be determined.
6. The method of claim 1, wherein the petrochemical supply chain distribution area comprises: the hemiversine formula in S6 is:
d＝R·c
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110157188.0A CN112766862A (en) | 2021-02-04 | 2021-02-04 | Method for dividing and locating distribution area of petrochemical product supply chain |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110157188.0A CN112766862A (en) | 2021-02-04 | 2021-02-04 | Method for dividing and locating distribution area of petrochemical product supply chain |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112766862A true CN112766862A (en) | 2021-05-07 |
Family
ID=75705014
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110157188.0A Pending CN112766862A (en) | 2021-02-04 | 2021-02-04 | Method for dividing and locating distribution area of petrochemical product supply chain |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112766862A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109636002A (en) * | 2018-11-13 | 2019-04-16 | 深圳大学 | A kind of site selecting method and system of the Cold Chain Logistics home-delivery center based on goods damage |
| CN110298474A (en) * | 2019-04-28 | 2019-10-01 | 武汉理工大学 | A kind of logistics centers location method of the ship stores spare part based on AIS and GIS |
| CN110533354A (en) * | 2018-05-25 | 2019-12-03 | 步步高商业连锁股份有限公司 | A kind of chain shop allied logistic distribution partitioning method based on equal goods metric density line |
-
2021
- 2021-02-04 CN CN202110157188.0A patent/CN112766862A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110533354A (en) * | 2018-05-25 | 2019-12-03 | 步步高商业连锁股份有限公司 | A kind of chain shop allied logistic distribution partitioning method based on equal goods metric density line |
| CN109636002A (en) * | 2018-11-13 | 2019-04-16 | 深圳大学 | A kind of site selecting method and system of the Cold Chain Logistics home-delivery center based on goods damage |
| CN110298474A (en) * | 2019-04-28 | 2019-10-01 | 武汉理工大学 | A kind of logistics centers location method of the ship stores spare part based on AIS and GIS |
Non-Patent Citations (6)
| Title |
|---|
| 史婉莹: ""我国成品油分销油库配送半径优化研究"", 《中国优秀硕士学位论文全文数据库 经济与管理科学辑》 * |
| 王亮: ""成品油物流系统的构建与优化"", 《中国博士学位论文全文数据库 经济与管理科学辑》 * |
| 王柳娜: "基于供应链的成品油二次物流优化研究", 《发现》 * |
| 田源等: "铁路运营物资区域联合储备中心选址", 《物流技术》 * |
| 白冰等: "石化产品配送区域划分与油库选址调整的优化解决方案研究", 清华大学出版社等 * |
| 陈卫: "论成品油二次配送优化与企业核心竞争力", 《辽宁经济职业技术学院(辽宁经济管理干部学院学报)》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Messaoudi et al. | Site selection methodology for the wind-powered hydrogen refueling station based on AHP-GIS in Adrar, Algeria | |
| CN103278833A (en) | Line recommendation system and method based on Beidou satellite/GPS (global positioning system) data | |
| Trench | How pipelines make the oil market work–Their networks, operation and regulation | |
| Bagočius et al. | Selecting a location for a liquefied natural gas terminal in the Eastern Baltic Sea | |
| CN108647821B (en) | Parameter self-learning-based differential evolution logistics distribution path optimization method | |
| CN112766862A (en) | Method for dividing and locating distribution area of petrochemical product supply chain | |
| CN116542591A (en) | Multi-mode intermodal route optimization method considering uncertain conditions | |
| CN114492904A (en) | Transportation path optimization method of logistics management system | |
| Myint et al. | Myanmar country report | |
| CN115841032A (en) | Method for calculating optimal bus motorized sharing rate of medium and small cities | |
| CN109543917B (en) | Underground logistics system network planning method based on 0-1 planning and prim algorithm | |
| CN111652782B (en) | Rural garbage collection operation system based on multiple metadata | |
| Ramalho et al. | External costs in short sea shipping based intermodal transport chains | |
| Svenson | Optimized route selection for logging trucks | |
| Timilsina et al. | Why have CO2 emissions increased in the transport sector in Asia? Underlying factors and policy options | |
| Teng et al. | Container multimodal transport routing selection modeling considering carbon emission cost | |
| Lomotko et al. | Modeling the railway and automobile supply chain on the basis of «green» logistics | |
| Bos | Direction indirect | |
| Myers Jaffe | Exploring the role of natural gas in US Trucking (Revised Version) | |
| Uemura et al. | Formulating Policy Options for Promoting Natural Gas Utilization in the East Asia Summit | |
| Boske et al. | Analysis of Texas biofuel supply chains Originating in the United States and Brazil | |
| CN114757570A (en) | LNG intelligent allocation management method based on big data and Internet of things system | |
| MIROSHNICHENKO | Modernization of transport and Customs infrastructure for the supply of liquefied natural gas market in Ukraine | |
| Svetlana | TRANSPORT IN INDONESIA: RECENT DEVELOPMENT TRENDS | |
| Lee et al. | Impacts of a Snake River drawdown on energy consumption and environmental emissions in transporting eastern Washington wheat and barley |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210507 |
|
| RJ01 | Rejection of invention patent application after publication |