CN114236101B - Method for calculating sulfur content of ship fuel without measuring carbon concentration - Google Patents
Method for calculating sulfur content of ship fuel without measuring carbon concentration Download PDFInfo
- Publication number
- CN114236101B CN114236101B CN202111439933.7A CN202111439933A CN114236101B CN 114236101 B CN114236101 B CN 114236101B CN 202111439933 A CN202111439933 A CN 202111439933A CN 114236101 B CN114236101 B CN 114236101B
- Authority
- CN
- China
- Prior art keywords
- ship
- solar
- diffusion
- fuel consumption
- determining
- 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.)
- Active
Links
- 239000000446 fuel Substances 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 39
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 title claims abstract description 23
- 229910052717 sulfur Inorganic materials 0.000 title claims abstract description 23
- 239000011593 sulfur Substances 0.000 title claims abstract description 23
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 title claims abstract description 15
- 229910052799 carbon Inorganic materials 0.000 title claims abstract description 15
- 238000009792 diffusion process Methods 0.000 claims abstract description 60
- 239000003344 environmental pollutant Substances 0.000 claims abstract description 40
- 231100000719 pollutant Toxicity 0.000 claims abstract description 40
- 238000012544 monitoring process Methods 0.000 claims abstract description 34
- 239000000295 fuel oil Substances 0.000 claims description 9
- 239000003921 oil Substances 0.000 claims description 3
- 229910052698 phosphorus Inorganic materials 0.000 claims description 2
- 238000007689 inspection Methods 0.000 abstract description 2
- 238000009434 installation Methods 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 18
- MWUXSHHQAYIFBG-UHFFFAOYSA-N nitrogen oxide Inorganic materials O=[N] MWUXSHHQAYIFBG-UHFFFAOYSA-N 0.000 description 9
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000001514 detection method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- GJEAMHAFPYZYDE-UHFFFAOYSA-N [C].[S] Chemical compound [C].[S] GJEAMHAFPYZYDE-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 239000000809 air pollutant Substances 0.000 description 1
- 231100001243 air pollutant Toxicity 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 150000003568 thioethers Chemical class 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/26—Oils; Viscous liquids; Paints; Inks
- G01N33/28—Oils, i.e. hydrocarbon liquids
- G01N33/2835—Specific substances contained in the oils or fuels
- G01N33/287—Sulfur content
Landscapes
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medicinal Chemistry (AREA)
- Oil, Petroleum & Natural Gas (AREA)
- General Chemical & Material Sciences (AREA)
- Food Science & Technology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for calculating sulfur content of ship fuel without measuring carbon concentration, which aims at remote sensing monitoring of ship tail gas, and determines atmospheric stability (the more unstable the atmospheric stability is, the stronger the diffusion capability is) according to conventional meteorological conditions such as cloud cover, wind speed, sunlight and the like, so as to further determine atmospheric diffusion parameters, and inversion is carried out on emission rate of emission source pollutants based on Gaussian line source model through the diffusion parameters and the measured tail gas concentration; and determining the fuel consumption of the ship according to the related parameters of the engine of the ship, and calculating the sulfur content of the ship fuel according to the emission rate of the strong source pollutants. The method mainly solves the problems of high difficulty, low efficiency and the like in the traditional boarding inspection, and simultaneously solves the problems of narrow coverage, high economic cost, large influence on the monitoring result by the installation position of equipment and the like in the main stream telemetry method, and has the characteristics of economy, high efficiency and wide coverage.
Description
Technical Field
The invention relates to a method for calculating sulfur content of ship fuel oil without measuring carbon concentration, and belongs to the technical field of remote sensing and monitoring of ship tail gas.
Background
In recent years, with the acceleration of urban process and the high-speed development of economy in China, the atmospheric pollution is gradually increased. The exhaust emission pollutant of the ship is one of the atmospheric pollution sources, and the main types of the exhaust emission pollutant comprise sulfur dioxide, nitrogen oxides, atmospheric particulates and the like. According to the statistics of the related data, sulfides and nitrogen oxides discharged by the ship account for nearly 10% of the total discharge, so that the sustainable development of socioeconomic performance is restricted, and the physical health of people is threatened. According to the notification of the transportation department about the implementation of the emission control area of the air pollutants of the stamp-pad ship, the emission control area is expanded to the inland river on the basis of the established emission control area of the Bohai sea (Jinjin Ji), long triangle and bead triangle water area ships; from 1 month 1 day 2020, the marine vessel enters a inland river control area, and marine fuel with sulfur content not more than 0.1% m/m is used; the emission of nitrogen oxides meets the corresponding emission limit requirement of International convention for preventing pollution of ships.
How to effectively monitor the emission of pollutants from a ship in a ship and to maximize the effect of an emission control zone is an important and urgent challenge in the current or even future period. Due to the particularity of monitoring the tail gas emission pollutants of the ship (the sailing water area is wide, the flowing pollution source and the meteorological environment are complex), and the problems of great difficulty, low efficiency and the like in the traditional boarding inspection are solved, the current prior art cannot meet the high-standard monitoring requirement, and the maritime department also lacks effective monitoring and management means for the tail gas emission of the ship. Therefore, the research on the remote measurement technology of the pollutant in the tail gas emission of the ship has important practical significance.
Developed countries and regions such as europe and america start relatively early in the aspect of monitoring the exhaust emission of ships, and research on monitoring the exhaust emission pollutants of ships is carried out in the respectively established ship emission control regions. Sniffing is still a commonly adopted technical method at the present stage, and in addition, research on the direction of optical telemetry is also gradually carried out. The sniffing technology is used for calculating the sulfur content of the fuel oil based on the sulfur-carbon content ratio in the ship tail gas, and the optical telemetry technology is used for inverting the gas types and the gas concentrations by utilizing the characteristic absorption property of the gas components in a spectrum band. According to the layout mode of the equipment, the equipment is divided into a fixed type and a movable type. Fixed equipment is typically loaded in fixed locations on land, ports, bridges, etc.; the mobile device is mounted on a mobile platform such as a ship or an unmanned aerial vehicle. From the device type, mainly comprises a laser radar, a differential spectrometer, a sniffer, a portable multi-gas meter and the like. In the monitoring equipment, the laser radar has wide monitoring coverage, but has higher manufacturing cost, maintenance cost and updating cost, and is suitable for the dock areas with centralized dock berths; the sniffing method is mature, but belongs to passive detection of gas concentration, the position relationship between equipment and pollution sources has larger influence on the result, and the screening limit value is difficult to determine, so that the sniffing method is mainly applicable to harbors with the leading wind direction of sea and land for a long time; the portable multi-gas meter is rapid in detection and relatively reliable in result, but has strong operation specialization and lower supervision coverage rate due to the need of boarding; the differential spectrometer has the advantages of quick detection, high automation degree, wide coverage and strong stability, but the related research in the current stage is less, and the practical application is still in the starting stage. According to the practical experience of the aspect of monitoring the ship tail gas outside the country, the sniffing method is a mainstream ship tail gas monitoring technical means at present, but the sniffing method has the defects of narrow coverage, high economic cost, large influence of equipment installation positions on monitoring results and the like due to the limitation of a monitoring principle.
Disclosure of Invention
The technical problems to be solved by the invention are as follows: the method for calculating the sulfur content of the ship fuel without measuring the carbon concentration has the advantages of wide coverage, low economic cost, high monitoring efficiency and the like.
The invention adopts the following technical scheme for solving the technical problems:
a method for calculating sulfur content of ship fuel without measuring carbon concentration, comprising the following steps:
step 1, determining a solar time angle according to longitude and observation time of a monitoring point, determining a solar altitude angle by combining a solar inclination angle and latitude of the monitoring point, determining a solar incidence level by combining cloud amount data, and determining atmospheric stability by combining ground wind speed data;
step 2, inquiring a power function value table of diffusion parameters according to the atmospheric stability and the leeward distance to determine diffusion parameters of the emission of the ship pollutants;
step 3, determining pollutant discharge rate according to a Gaussian line source diffusion equation, diffusion parameters and measured pollutant concentration;
step 4, determining the fuel consumption of the ship according to the basic fuel consumption coefficient of the ship main engine, the fuel consumption of the auxiliary engine, the power of the main engine, the fuel consumption rate of the main engine, the power of the auxiliary engine and the fuel consumption rate of the auxiliary engine;
and 5, determining the sulfur content of the ship fuel according to the pollutant discharge rate and the ship fuel consumption.
As a preferable scheme of the invention, the specific process of the step 1 is as follows:
step 11, determining a solar time angle according to the longitude of the monitoring point and the observation time:
ω＝15t+λ-300
wherein ω represents a solar time angle; t represents an hour part of Beijing time at the time of observation; λ represents the geographical longitude of the location of the monitoring point;
step 12, determining a solar altitude angle according to the solar hour angle, the solar inclination angle and the latitude of the monitoring point:
h＝arcsin(sinψsinσ+cosψcosσcosω)
wherein h represents the solar altitude; psi represents the geographic latitude of the position of the monitoring point; sigma represents the solar tilt angle;
wherein, the value of theta isd represents the date on which the observation time is located in one year;
step 13, inquiring a solar incidence grade table according to solar altitude and cloud cover data, and determining a solar incidence grade;
and 14, inquiring an atmospheric stability level table according to the solar incidence level and the ground wind speed data, and determining the atmospheric stability.
As a preferable scheme of the invention, the specific process of the step 2 is as follows:
the diffusion parameter is a function of the atmospheric stability and the lee distance in the atmospheric diffusion mode, expressed exponentially:
in sigma y Sum sigma z Respectively representing transverse diffusion parameters and longitudinal diffusion parameters, X represents the distance between the downwind direction and the pollution source, and gamma 1 、α 1 Respectively representing the index and coefficient of the calculated lateral diffusion parameter, gamma 2 、α 2 Respectively representing an index and a coefficient for calculating the longitudinal diffusion parameter;
γ 1 、α 1 、γ 2 and alpha 2 Is determined from a power function value table of the diffusion parameter.
As a preferable scheme of the invention, the specific process of the step 3 is as follows:
when the pollutant is continuously discharged along the horizontal direction, the pollutant is regarded as a line source, the line source with limited length is divided into n equal parts according to a point source diffusion equation, and when the included angle between the wind direction and the line source is alpha, a wind system coordinate is introduced, and y is the same as that of the wind system coordinate w The axis is wind direction, x w The axis is the direction perpendicular to the wind direction, deducing a limited-length line source diffusion model to calculate a downwind monitoring point P (x p ,y p The concentration formula at z) is:
wherein C (x) p ,y p Z, α) represents the concentration of the downwind P point, Q represents the pollutant discharge rate, u represents the average wind speed, L represents the length of the pollutant source, σ y Sum sigma z Respectively represents a transverse diffusion parameter and a longitudinal diffusion parameter, x p ,y p Z represents the coordinate of a downwind monitoring point P, h' represents the height of a discharge source from the ground, alpha represents the included angle between the source and the wind direction, and x 1 The downwind distance between a certain point P '(0, y', z) and the point P on the line source is given; the pollutant discharge rate is determined according to the concentration formula.
As a preferred scheme of the present invention, the specific process of the step 4 is as follows:
S＝p*A 1 +A 2
A 1 ＝a 1 *b 1
A 2 ＝a 2 *b 2
wherein S is the fuel consumption of the ship; p is the basic oil consumption of the main engineCoefficients; a is that 1 The fuel consumption of the host; a is that 2 The fuel consumption of auxiliary machinery; a, a 1 Is the host power; b 1 Is the fuel consumption rate of the main engine; a, a 2 Is auxiliary power; b 2 Is the fuel consumption rate of auxiliary machinery.
As a preferable scheme of the invention, the specific process of the step 5 is as follows:
wherein Q represents the pollutant discharge rate, and S represents the fuel consumption of the ship.
Compared with the prior art, the technical scheme provided by the invention has the following technical effects:
according to the method, the ship tail gas monitoring points are set, after tail gas is measured, the components and the calculation method of the atmospheric stability are defined by combining key parameters such as weather, date and coordinates. The continuous leakage source characteristic of the tail gas emission of the ship is combined, a Gaussian line source diffusion model is built, and the pollutant emission source is reversely pushed out; and the instrument is used for selecting a sample ship to collect the tail gas emission source intensity data in a short distance, the established tail gas pollutant diffusion model and the emission source intensity calculation method are verified and calibrated, a reverse thrust algorithm of the emission value is determined, and the sulfur content in the fuel oil is reversely deduced by using the measured sulfur dioxide concentration. Avoiding high error of calculation of carbon concentration when sulfur content is measured by sniffing.
Drawings
FIG. 1 is a flow chart of a method for calculating sulfur content of ship fuel without measuring carbon concentration;
FIG. 2 is a graph of the atmospheric stability calculation process of the method of the present invention;
FIG. 3 is a graph of a model of a Gaussian line source according to the method of the invention.
Detailed Description
Embodiments of the present invention are described in detail below, examples of which are illustrated in the accompanying drawings. The embodiments described below by referring to the drawings are exemplary only for explaining the present invention and are not to be construed as limiting the present invention.
Referring to fig. 1, a flowchart of a method for calculating sulfur content of ship fuel without measuring carbon concentration according to the present invention comprises the following specific steps:
step 1: determining the atmospheric stability: the solar time angle is determined through the local longitude and the observation time, the solar altitude angle is determined through the solar dip angle and the local latitude, the solar incidence level is determined through cloud amount data, and the atmospheric stability is determined through the wind speed.
As shown in fig. 2, the atmospheric stability determination method includes the steps of:
step 11: the solar time angle is determined from the local longitude and the observation time.
ω＝15t+λ-300
Wherein ω represents a solar time angle; t represents an hour part of Beijing time at the time of observation; λ represents the geographical longitude of the location.
Step 12: and determining the solar altitude according to the solar time angle, the solar dip angle and the local latitude.
Wherein σ represents the solar tilt angle; the value of theta isd represents the day of the year on which the current date is.
h＝arcsin(sinψsinσ+cosψcosσcosω)
Wherein h represents the solar altitude; psi represents the geographical latitude of the location.
Step 13: the solar incidence level was determined from a solar altitude and cloud cover look-up table (table 1).
Table 1 solar incidence level table
Step 13: atmospheric stability level was determined from sun incidence level and ground wind speed table look-up (table 2).
Table 2 atmospheric stability level table
Step 2: determining diffusion parameters: the diffusion parameter is a function of the atmospheric stability and the leeward distance in the atmospheric diffusion mode. And determining the diffusion parameters of the pollutant emission according to the power function value table of the diffusion parameters.
The diffusion parameter is a function of the atmospheric stability and the leeward distance in the atmospheric diffusion mode. Reference to national standard GB/T13201-91, expressed in exponential form:
in sigma y Sum sigma z Respectively representing transverse diffusion parameters and longitudinal diffusion parameters, X represents the distance between the downwind direction and the pollution source, and gamma 1 、α 1 Respectively representing the index and coefficient of the calculated lateral diffusion parameter, gamma 2 、α 2 Respectively representing the index and the coefficient of the calculated longitudinal diffusion parameter. Gamma ray 1 、α 1 、γ 2 、α 2 Is determined from a power function value table (table 3) of the diffusion parameters.
TABLE 3 diffusion parameter power function table
Step 3: the pollutant discharge rate is determined from the gaussian line source diffusion equation, the diffusion parameters and the measured pollutant concentration.
And determining the emission source pollutant emission rate according to the Gaussian line source diffusion model. When the pollutant is continuously discharged along the horizontal direction, the pollutant can be regarded as a line source, and the pollutant concentration of the line source discharged in the transverse wind direction is generally equal, so that the variable y can be integrated by a Gaussian model of point source diffusion (point source summation method), and the Gaussian diffusion model of the line source can be obtained.
According to the point source diffusion equation, the finite length line source is divided into n equal parts, when the included angle between the wind direction and the line source is alpha, a wind system coordinate is introduced (figure 3), and the integral mode of the finite length line source diffusion when the wind direction is arbitrary is deduced.
The deduced finite length line source diffusion mode calculates the downwind P-point (x p ,y p The concentration formula at z) is:
wherein C (x) p ,y p Z, α) represents the concentration of the downwind P point, Q represents the pollutant discharge rate, L represents the length of the pollutant source, u represents the average wind speed, σ y Sum sigma z Respectively a transverse diffusion parameter and a longitudinal diffusion parameter, h' represents the height of a discharge source from the ground, alpha represents the included angle between the source and the wind direction, and x 1 Y are a point P '(0, y', z) and a calculated point P (x) p ,y p Z) downwind and lateral distance between.
Step 4: and determining the fuel consumption of the ship according to the basic fuel consumption coefficient of the main engine, the fuel consumption of the auxiliary engine, the main engine power, the main engine fuel consumption rate, the auxiliary engine power and the auxiliary engine fuel consumption rate.
And determining the fuel consumption of the ship by combining the navigational speed information of the ship AIS data and the information of the tonnage, the host, the auxiliary machine, the rated power and the like of the ship in the ship type database. The implementation process is as follows:
S＝p*A 1 +A 2
A 1 ＝a 1 *b 1
A 2 ＝a 2 *b 2
wherein S is the fuel consumption of the ship; p is the basic oil consumption coefficient of the host; a is that 1 The fuel consumption of the host; a is that 2 The fuel consumption of auxiliary machinery; a, a 1 Is the host power; b 1 Is the fuel consumption rate of the main engine; a, a 2 Is auxiliary power; b 2 Is the fuel consumption rate of auxiliary machinery.
Step 5: and determining the sulfur content of the ship fuel according to the emission rate of the emission source pollutant and the ship fuel consumption.
Because of S and SO 2 The molecular mass ratio of (2) is 32:64=1/2, so
Examples
2019, 12, 25 am 10:06 min-10: 07 th division, jiangsu province Nanjing city in Yangtze river three bridge monitoring point monitoring a ship SO 2 The concentration was 51.5. Mu.g/m 3 And the system calculates to identify the ship as out of standard and sends out alarm. After reporting to Jiangsu maritime bureau, the ship is detected by maritime related departments, and the ship is verified to be a ship with exceeding fuel sulfur content.
Local longitude 118.62614, latitude 31.961408. Then t=10, d=359, λ= 118.6261, ψ= 31.9614.
ω＝15t+λ-300＝-31.3739
h＝arcsin(sinψsinσ+sosψcosσcosω)＝21.1266
And obtaining the atmospheric stability between the D-E stages by checking and combining cloud quantity data and wind speed conditions.
The diffusion parameter is a function of the atmospheric stability and the leeward distance in the atmospheric diffusion mode. Reference to national standard GB/T13201-91, sigma y 、σ z Expressed in exponential form:
wherein X represents the distance from the downwind to the pollution source, gamma 1 、γ 2 、α 1 、α 2 By looking up Table 3, in case of wind, the atmospheric stability is D-E class, gamma 1 ＝0.078563，α 1 ＝0.925118，γ 2 ＝0.092753，α 2 ＝0.78837。
Combining the parameters, and calculating by using a Gaussian line source model to obtain: x is x 1 Emission source emission rate Q was 156642.7491 μg/s (0.156643 g/s) =10 meters.
The power of the main engine of the ship is 200kw, the fuel consumption rate of the main engine is 0.192 kg/kw.h, the power of the auxiliary engine is 21kw, the fuel consumption rate of the auxiliary engine is 0.185 kg/kw.h, the ship is a bulk carrier, and the basic fuel consumption coefficient of the main engine can take a value of 0.72. Then
S＝p*a 1 *b 1 +a 2 *b 2 ＝31.53kg/h
Namely, the fuel consumption of the ship is 31.53kg/h (8.76 g/s).
The above embodiments are only for illustrating the technical idea of the present invention, and the protection scope of the present invention is not limited thereto, and any modification made on the basis of the technical scheme according to the technical idea of the present invention falls within the protection scope of the present invention.
Claims (6)
1. A method for calculating sulfur content of ship fuel without measuring carbon concentration, which is characterized by comprising the following steps:
step 1, determining a solar time angle according to longitude and observation time of a monitoring point, determining a solar altitude angle by combining a solar inclination angle and latitude of the monitoring point, determining a solar incidence level by combining cloud amount data, and determining atmospheric stability by combining ground wind speed data;
step 2, inquiring a power function value table of diffusion parameters according to the atmospheric stability and the leeward distance to determine diffusion parameters of the emission of the ship pollutants;
step 3, determining pollutant discharge rate according to a Gaussian line source diffusion equation, diffusion parameters and measured pollutant concentration;
step 4, determining the fuel consumption of the ship according to the basic fuel consumption coefficient of the ship main engine, the fuel consumption of the auxiliary engine, the power of the main engine, the fuel consumption rate of the main engine, the power of the auxiliary engine and the fuel consumption rate of the auxiliary engine;
and 5, determining the sulfur content of the ship fuel according to the pollutant discharge rate and the ship fuel consumption.
2. The method for calculating the sulfur content of the ship fuel oil without measuring the carbon concentration according to claim 1, wherein the specific process of the step 1 is as follows:
step 11, determining a solar time angle according to the longitude of the monitoring point and the observation time:
ω＝15t+λ-300
wherein ω represents a solar time angle; t represents an hour part of Beijing time at the time of observation; λ represents the geographical longitude of the location of the monitoring point;
step 12, determining a solar altitude angle according to the solar hour angle, the solar inclination angle and the latitude of the monitoring point:
h＝arcsin(sinψsinσ+cosψcosσcosω)
wherein h represents the solar altitude; psi represents the geographic latitude of the position of the monitoring point; sigma represents the solar tilt angle;
wherein, the value of theta isd represents the date on which the observation time is located in one year;
step 13, inquiring a solar incidence grade table according to solar altitude and cloud cover data, and determining a solar incidence grade;
and 14, inquiring an atmospheric stability level table according to the solar incidence level and the ground wind speed data, and determining the atmospheric stability.
3. The method for calculating the sulfur content of the ship fuel oil without measuring the carbon concentration according to claim 1, wherein the specific process of the step 2 is as follows:
the diffusion parameter is a function of the atmospheric stability and the lee distance in the atmospheric diffusion mode, expressed exponentially:
in sigma y Sum sigma z Respectively representing transverse diffusion parameters and longitudinal diffusion parameters, X represents the distance between the downwind direction and the pollution source, and gamma 1 、α 1 Respectively representing the index and coefficient of the calculated lateral diffusion parameter, gamma 2 、α 2 Respectively representing an index and a coefficient for calculating the longitudinal diffusion parameter;
γ 1 、α 1 、γ 2 and alpha 2 Is determined from a power function value table of the diffusion parameter.
4. The method for calculating the sulfur content of the ship fuel oil without measuring the carbon concentration according to claim 1, wherein the specific process of the step 3 is as follows:
when the pollutant is continuously discharged along the horizontal direction, the pollutant is regarded as a line source, the line source with limited length is divided into n equal parts according to a point source diffusion equation, and when the included angle between the wind direction and the line source is alpha, a wind system coordinate is introduced, and y is the same as that of the wind system coordinate w The axis is wind direction, x w The axis is the direction perpendicular to the wind direction, deducing the limited length line source diffusion mode to calculate the wind direction monitoring point P (x p ，y p The concentration formula at z) is:
wherein C (x) p ，y p Z, α) represents the concentration of the downwind P point, Q represents the pollutant discharge rate, u represents the average wind speed, L represents the length of the pollutant source, σ y Sum sigma z Respectively represents a transverse diffusion parameter and a longitudinal diffusion parameter, x p ，y p Z represents the coordinate of a downwind monitoring point P, h' represents the height of a discharge source from the ground, alpha represents the included angle between the source and the wind direction, and x 1 The downwind distance between a certain point P '(0, y', z) and the point P on the line source is given; the pollutant discharge rate is determined according to the concentration formula.
5. The method for calculating sulfur content of ship fuel oil without measuring carbon concentration according to claim 1, wherein the specific process of step 4 is as follows:
S＝p*A 1 +A 2
A 1 ＝a 1 *b 1
A 2 ＝a 2 *b 2
wherein S is the fuel consumption of the ship; p is the basic oil consumption coefficient of the host; a is that 1 The fuel consumption of the host; a is that 2 The fuel consumption of auxiliary machinery; a, a 1 Is the host power; b 1 Is the fuel consumption rate of the main engine; a, a 2 Is auxiliary power; b 2 Is the fuel consumption rate of auxiliary machinery.
6. The method for calculating sulfur content of ship fuel oil without measuring carbon concentration according to claim 1, wherein the specific process of step 5 is as follows:
wherein Q represents the pollutant discharge rate, and S represents the fuel consumption of the ship.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111439933.7A CN114236101B (en) | 2021-11-30 | 2021-11-30 | Method for calculating sulfur content of ship fuel without measuring carbon concentration |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111439933.7A CN114236101B (en) | 2021-11-30 | 2021-11-30 | Method for calculating sulfur content of ship fuel without measuring carbon concentration |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114236101A CN114236101A (en) | 2022-03-25 |
| CN114236101B true CN114236101B (en) | 2023-11-28 |
Family
ID=80752026
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111439933.7A Active CN114236101B (en) | 2021-11-30 | 2021-11-30 | Method for calculating sulfur content of ship fuel without measuring carbon concentration |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114236101B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107589100A (en) * | 2017-09-08 | 2018-01-16 | 交通运输部天津水运工程科学研究所 | A kind of bunker oil sulfur content sniff estimation algorithm |
| CN109060014A (en) * | 2018-07-16 | 2018-12-21 | 交通运输部水运科学研究所 | A kind of marine exhaust discharge compliance determines and violation ship identification method and system |
| KR20200086538A (en) * | 2019-01-09 | 2020-07-17 | 한국해양과학기술원 | Emission quantity monitoring system and the estimating method for find dust of ship |
| CN111582580A (en) * | 2020-05-08 | 2020-08-25 | 北京交通大学 | Travel path planning method considering population pollutant exposure |
| CN113125683A (en) * | 2021-04-23 | 2021-07-16 | 武汉理工大学 | Mobile monitoring device and method for sulfur content of marine fuel oil |
| CN113125682A (en) * | 2021-04-23 | 2021-07-16 | 武汉理工大学 | Buoy type device and method for monitoring sulfur content of fuel oil for regional multi-ship |
| CN113360850A (en) * | 2021-06-17 | 2021-09-07 | 东莞市环境监测中心站 | Atmospheric pollutant source tracing analysis method |
-
2021
- 2021-11-30 CN CN202111439933.7A patent/CN114236101B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107589100A (en) * | 2017-09-08 | 2018-01-16 | 交通运输部天津水运工程科学研究所 | A kind of bunker oil sulfur content sniff estimation algorithm |
| CN109060014A (en) * | 2018-07-16 | 2018-12-21 | 交通运输部水运科学研究所 | A kind of marine exhaust discharge compliance determines and violation ship identification method and system |
| KR20200086538A (en) * | 2019-01-09 | 2020-07-17 | 한국해양과학기술원 | Emission quantity monitoring system and the estimating method for find dust of ship |
| CN111582580A (en) * | 2020-05-08 | 2020-08-25 | 北京交通大学 | Travel path planning method considering population pollutant exposure |
| CN113125683A (en) * | 2021-04-23 | 2021-07-16 | 武汉理工大学 | Mobile monitoring device and method for sulfur content of marine fuel oil |
| CN113125682A (en) * | 2021-04-23 | 2021-07-16 | 武汉理工大学 | Buoy type device and method for monitoring sulfur content of fuel oil for regional multi-ship |
| CN113360850A (en) * | 2021-06-17 | 2021-09-07 | 东莞市环境监测中心站 | Atmospheric pollutant source tracing analysis method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114236101A (en) | 2022-03-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Yang et al. | An AIS-based emission inventory and the impact on air quality in Tianjin port based on localized emission factors | |
| CN112749478A (en) | Atmospheric pollution source-tracing diffusion analysis system and method based on Gaussian diffusion model | |
| CN109978335B (en) | Ship fuel compliance judgment and violation identification method | |
| CN109060014A (en) | A kind of marine exhaust discharge compliance determines and violation ship identification method and system | |
| CN113624921B (en) | Multi-mode comprehensive pollution tracing method | |
| CN111814367A (en) | Port ship emission monitoring and supervision cloud service system | |
| CN111811572A (en) | Ship exhaust emission real-time monitoring method based on big data | |
| Yang et al. | Evaluation of the effect of air pollution control during the Beijing 2008 Olympic Games using Lidar data | |
| CN105539797A (en) | Navigation method and system of wind power assisting navigation ship based on ECDIS | |
| Wang et al. | Analysis of ship emission characteristics under real-world conditions in China | |
| CN112214721B (en) | AIS data-based dynamic ship emission list establishing method | |
| Zou et al. | Effects of cleaner ship fuels on air quality and implications for future policy: A case study of Chongming Ecological Island in China | |
| Sun et al. | Ship emission of nitrous acid (HONO) and its impacts on the marine atmospheric oxidation chemistry | |
| CN109975492B (en) | Coastal atmosphere combined pollution sky-space-ground integrated monitoring and early warning system | |
| CN114819704A (en) | Tracing method and system for mobile source pollutant emission | |
| CN114236101B (en) | Method for calculating sulfur content of ship fuel without measuring carbon concentration | |
| He et al. | Selection of onshore sites based on monitoring possibility evaluation of exhausts from individual ships for Yantian Port, China | |
| CN114002381A (en) | Atmospheric pollution source-tracing diffusion analysis method and device | |
| CN114239426A (en) | Yangtze river trunk ship emission list generation method based on water flow data assimilation | |
| Wu et al. | Development of a spectrum-based ship fuel sulfur content real-time evaluation method | |
| Hu et al. | Using the OSPM model on pollutant dispersion in an urban street canyon | |
| Peng et al. | System for identifying ships using high-sulfur-content fuel oil in port waters | |
| Deng et al. | SO2 compliance monitoring and emission characteristics analysis of navigating ships: A case study of Shanghai waters in emission control areas, China | |
| CN116465824A (en) | Ship fuel oil sulfur content calculation method based on optical technology and sulfur-carbon ratio method | |
| CN112364300A (en) | Near-ground wind speed statistics downscaling correction method based on relative slope length |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |