CN113419037B - Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas - Google Patents

Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas Download PDF

Info

Publication number
CN113419037B
CN113419037B CN202110982164.9A CN202110982164A CN113419037B CN 113419037 B CN113419037 B CN 113419037B CN 202110982164 A CN202110982164 A CN 202110982164A CN 113419037 B CN113419037 B CN 113419037B
Authority
CN
China
Prior art keywords
ship
tail gas
concentration
target
information
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
Application number
CN202110982164.9A
Other languages
Chinese (zh)
Other versions
CN113419037A (en
Inventor
邱梦春
陆生忠
温作乐
梁大荣
汪鲁见
吴垒
王超
屈颖
周城
陶波
于志伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Chunlai Technology Co Ltd
Original Assignee
Hangzhou Chunlai Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Chunlai Technology Co Ltd filed Critical Hangzhou Chunlai Technology Co Ltd
Priority to CN202110982164.9A priority Critical patent/CN113419037B/en
Publication of CN113419037A publication Critical patent/CN113419037A/en
Application granted granted Critical
Publication of CN113419037B publication Critical patent/CN113419037B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/0004—Gaseous mixtures, e.g. polluted air
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization

Abstract

The invention relates to a method, a terminal and a system for monitoring concentration of pollutants discharged by multi-ship tail gas, wherein the monitoring method comprises the following steps: acquiring the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter, and establishing a tail gas information matrix; acquiring static information and dynamic information of a ship and establishing a ship information matrix; identifying a target time interval in which the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter has a peak value, and intercepting a corresponding tail gas information sub-matrix from a tail gas information matrix; polling target time intervals corresponding to the ship tail gas telemeters, and screening target ships meeting target conditions in a ship information matrix; correlating the telemetering concentration and meteorological information of the atmospheric pollutants monitored by the target ship and at least two ship tail gas telemeters; and inverting the concentration of the atmospheric pollutants discharged by the tail gas of each ship based on a ship smoke plume Gaussian diffusion model and a nonlinear multi-component analysis algorithm. The invention realizes the concentration monitoring of the pollutants discharged by the tail gas of a plurality of ships.

Description

Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas
Technical Field
The invention belongs to the technical field of ship tail gas emission monitoring, and particularly relates to a method, a terminal and a system for monitoring the concentration of pollutants discharged by multiple ship tail gases.
Background
Along with the departure of the implementation scheme of the ship atmospheric pollutant emission control area, the ship emission control area is subjected to capacity expansion and expanded to inland rivers; and the emission of atmospheric pollutants such as sulfur oxides, nitrogen oxides, particulate matters, volatile organic compounds and the like is required to be reduced by measures such as improving the ship fuel standard, the ship engine standard and the like.
The implementation of the ship emission control zone embodiment requires effective monitoring means to provide security. The current monitoring means is single and low-efficiency, and the fuel oil sulfur content is only checked on board by a random way, so that the effective monitoring of a large number of ships entering and leaving a port cannot be realized. In addition, the existing ship tail gas monitoring method mainly comprises the following two steps:
the first is a method for remotely measuring ship exhaust gas by an sniffing method, and for example, a device and a method for remotely measuring sniffed exhaust gas are disclosed in patent document No. CN 113155921A. However, the remote measurement result is susceptible to environmental factors, such as wind direction, wind speed, ship speed and course, and has the defects of missing measurement, error measurement and the like.
The second is a method for remotely measuring the exhaust gas of a ship based on a spectrum remote measuring technique, for example, a inland ship monitoring system disclosed in patent document with publication number CN112132385A, which includes a spectrum exhaust gas analyzer. However, the tail gas analysis result is easy to install, and the measurement distance is too long (more than 200 m), so that CO is generated2Saturated absorption and the like;
moreover, the two ship tail gas remote measuring methods only can monitor the fuel oil sulfur content of a single ship, and cannot meet the monitoring requirement that multiple ships sail at the same time in a measuring area, which often exists in actual situations; in addition, there is no effective means for monitoring atmospheric pollutants such as sulfur oxides, nitrogen oxides, and carbon oxides emitted from ship exhaust.
Disclosure of Invention
Based on the above-mentioned shortcomings and drawbacks of the prior art, an object of the present invention is to solve at least one or more of the above-mentioned problems of the prior art, in other words, to provide a method, a terminal and a system for monitoring the concentration of pollutants in exhaust emission of multiple ships, which satisfy one or more of the above-mentioned needs.
In order to achieve the purpose, the invention adopts the following technical scheme:
a monitoring method for the concentration of pollutants discharged by multi-ship tail gas is characterized in that a plurality of ship tail gas telemeters are distributed in a ship discharge control area, and the telemetering areas of the ship tail gas telemeters are mutually independent;
the monitoring method comprises the following steps:
s1, acquiring the telemetered concentration of the atmospheric pollutants monitored by each ship tail gas telemeter, and establishing a tail gas information matrix formed by the time point and the telemetered concentration of the atmospheric pollutants;
acquiring static information and dynamic information of a ship, wherein the static information comprises a ship MMSI and inherent parameters of the ship, and the dynamic information comprises the ship MMSI, a current navigation time point, a navigation speed, a course and a ship coordinate position; according to the MMSI of the ship and the static information and the dynamic information of the ship integrated at the current sailing time point, establishing a ship information matrix;
s2, identifying a target time interval of the peak value of the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter, and intercepting in a tail gas information matrix according to the target time interval to obtain a corresponding tail gas information sub-matrix;
s3, polling target time intervals corresponding to the ship tail gas telemeters, and screening target ships meeting target conditions in a ship information matrix;
s4, correlating the target ship and the telemetering concentration of the atmospheric pollutants monitored by at least two ship tail gas telemeters corresponding to the target ship based on the tail gas information submatrix, and correlating the telemetering concentration of the atmospheric pollutants with meteorological information corresponding to a target time interval;
s5, inverting the concentration of the atmospheric pollutants discharged by the tail gas of each ship based on the ship smoke plume Gaussian diffusion model and the nonlinear multi-component analysis algorithm.
Preferably, in step S3, the target conditions include condition (1) and condition (2):
condition (1): condition (1):
Figure 227307DEST_PATH_IMAGE001
；tis the current point in time of the voyage of the ship,δfor maximum time of reception of dynamic information, TmaxIs the maximum value of the target time interval, TminIs the minimum value of the target time interval;
condition (2): the distance between the ship coordinate position and the corresponding correlation midpoint position of the ship tail gas telemeter is smaller than the effective detection distance of the corresponding ship tail gas telemeter.
Preferably, in step S3, the target conditions are:
the distance between the ship coordinate position and the corresponding correlation midpoint position of the ship tail gas telemeter is smaller than the effective detection distance of the corresponding ship tail gas telemeter.
Preferably, in step S4, the meteorological information includes a current environmental wind direction, a wind speed, a temperature, a humidity, and an air pressure.
Preferably, the step S5 includes:
based on the received vessel coordinate position of the target vessel: (x 0，y 0) Course of the vehiclewSpeed of flightvExpressed in a geodetic coordinate system, thentPosition coordinates of the voyage point of the target ship at the time (c:)X，Y）：
Figure 82130DEST_PATH_IMAGE002
Converting the position coordinates of the navigation point in the geodetic coordinate system into the ship coordinate position (x 0，y 0) Is at the origin and has a wind direction angle ofxAxis in the cross-wind directionyPosition coordinates of a navigation point in a wind direction coordinate system established by axes: (X′，Y′）：
Figure 65130DEST_PATH_IMAGE003
Wherein the content of the first and second substances,Lis the distance between the correlation midpoint position of the ship tail gas telemeter and a target ship,Pis the downwind angle;
suppose that the target vessel navigates tot n At the moment, based on a Gaussian plume diffusion model, the concentration of the atmospheric pollutants is calculated as follows:
Figure 398022DEST_PATH_IMAGE004
wherein the content of the first and second substances,
Figure 353340DEST_PATH_IMAGE005
for ships att i The atmospheric pollutants discharged at any moment diffuse tot n Any space point of time
Figure 379064DEST_PATH_IMAGE006
The concentration of the (C) in the (C),Q i for ships att i The total amount of atmospheric pollutants emitted at a time,
Figure 789973DEST_PATH_IMAGE007
are respectively as
Figure 926556DEST_PATH_IMAGE008
The diffusion coefficient in the direction of the light,uwhich is the wind speed,His a point in space
Figure 798697DEST_PATH_IMAGE006
Height from water surface;
the method comprises the steps of obtaining and accumulating the concentrations of the atmospheric pollutants at a plurality of space points in a target time interval, establishing a nonlinear equation with the telemetering concentration of the atmospheric pollutants monitored by a ship tail gas telemeter corresponding to a target ship, and inversing the concentrations of the atmospheric pollutants emitted by the target ships at different moments by utilizing a nonlinear multi-component analysis algorithm.
Preferably, the nonlinear multi-component analysis algorithm is a nonlinear least squares method or a principal component analysis method.
Preferably, the atmospheric pollutants are sulfur oxides, nitrogen oxides, carbon oxides or particulate matter.
Preferably, after the step S5, the method further includes the following steps:
s6, calculating the sulfur content of the fuel oil of the ship
Figure 933006DEST_PATH_IMAGE009
；
Wherein the content of the first and second substances,C S is the concentration of sulfur oxides in the vessel,C C is the concentration of carbon oxides of the vessel.
The invention also provides a terminal for monitoring the concentration of pollutants discharged by the tail gas of multiple ships, which comprises:
the acquisition module is used for acquiring the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter and establishing a tail gas information matrix formed by time points and the telemetering concentration of the atmospheric pollutants; the system is also used for acquiring static information and dynamic information of the ship, wherein the static information comprises the MMSI of the ship and inherent parameters of the ship, and the dynamic information comprises the MMSI of the ship, the current time point of navigation, the speed, the course and the coordinate position of the ship; according to the MMSI of the ship and the static information and the dynamic information of the ship integrated at the current sailing time point, establishing a ship information matrix;
the identification module is used for identifying a target time interval of the peak value of the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter, and intercepting the target time interval in a tail gas information matrix to obtain a corresponding tail gas information sub-matrix;
the screening module is used for polling the target time intervals corresponding to the ship tail gas telemeters and screening target ships meeting target conditions in the ship information matrix;
the correlation module is used for correlating the target ship with the telemetering concentration of the atmospheric pollutants monitored by at least two ship tail gas telemeters corresponding to the target ship based on the tail gas information submatrix and is also used for correlating the telemetering concentration of the atmospheric pollutants with meteorological information corresponding to a target time interval;
and the inversion module is used for inverting the concentration of the atmospheric pollutants discharged by the tail gas of each ship based on the ship smoke plume Gaussian diffusion model and the nonlinear multi-component analysis algorithm.
The invention also provides a monitoring system for the concentration of pollutants discharged by the tail gas of multiple ships, which comprises:
the ship tail gas telemeters are distributed in the ship emission control area; the remote measuring areas of the ship tail gas telemeters are mutually independent;
the monitoring terminal according to the above scheme is in communication connection with each ship tail gas telemeter;
the AIS receiver is in communication connection with an AIS transmitter and a monitoring terminal of the ship respectively;
and the meteorological analyzer is in communication connection with the monitoring terminal.
Compared with the prior art, the invention has the beneficial effects that:
(1) the invention realizes the concentration monitoring of the tail gas emission pollutants in the navigation process of a plurality of ships in the ship emission control area; compared with the existing single-ship monitoring technology, the monitoring efficiency is effectively improved;
(2) according to the invention, the corresponding tail gas information sub-matrix is obtained by intercepting in the tail gas information matrix according to the target time interval, so that the subsequent data processing amount is effectively reduced, the data processing efficiency is improved, and the requirement on the data processing capacity of the monitoring terminal is reduced.
(3) The invention can realize the monitoring of the concentration of the atmospheric pollutants such as sulfur oxide, nitrogen oxide, carbon oxide and the like discharged by the tail gas of each ship, and can also monitor the sulfur content of fuel oil of the ship.
Drawings
Fig. 1 is a flowchart of a method for monitoring the concentration of pollutants emitted from multiple ship exhaust according to embodiment 1 of the present invention;
FIG. 2 is a schematic diagram of converting a geodetic coordinate system into a wind direction coordinate system according to embodiment 1 of the present invention;
FIG. 3 is a schematic diagram of a multi-vessel exhaust emission pollutant concentration monitoring terminal according to embodiment 1 of the present invention;
FIG. 4 is a schematic diagram of the construction of a system for monitoring the concentration of pollutants discharged from multiple ships in accordance with embodiment 1 of the present invention;
FIG. 5 is a schematic plan view of a multi-vessel exhaust emission pollutant concentration monitoring system applied to a port area according to embodiment 1 of the present invention;
fig. 6 is a schematic plan layout view of a monitoring system for concentration of pollutants discharged from multiple ships according to embodiment 3 of the present invention, applied to an inland river region.
Detailed Description
In order to more clearly illustrate the embodiments of the present invention, the following description will explain the embodiments of the present invention with reference to the accompanying drawings. It is obvious that the drawings in the following description are only some examples of the invention, and that for a person skilled in the art, other drawings and embodiments can be derived from them without inventive effort.
Example 1:
the embodiment is exemplified and explained in a ship emission control area by taking a port area as an example, and is arrangedRThe remote measuring areas of the ship tail gas telemeters are mutually independent. Wherein the content of the first and second substances,Ris an integer greater than 1.
As shown in fig. 1, the method for monitoring the concentration of pollutants emitted from multiple ships in this embodiment includes the following steps:
s1, obtainingRThe telemetering concentration of the atmospheric pollutants monitored by the ship tail gas telemeter establishes a digital tail gas information matrix formed by time points and the telemetering concentration of the atmospheric pollutantsD 1、D 2、…、D R ；
Wherein the atmospheric pollutant includes SO2、CO2、NO x PM, etc., the time for the ship tail gas telemeter to acquire the telemetered concentration is generally 500ms, and then a tail gas information matrixD j The expression of (a) is as follows:
Figure 625019DEST_PATH_IMAGE010
wherein the content of the first and second substances,j=1，2，…，R，t m，j is as followsjShip tail gasFirst of telegauge monitoringmAt a point in time, the time of the measurement,mis an integer greater than 1;
Figure 565293DEST_PATH_IMAGE011
are respectively the firstjA ship tail gas telemeter is arranged on the secondmCO monitored at each time point2、NO x 、SO2And the remote measured concentration of PM.
The method comprises the steps that the AIS receiver is used for obtaining static information and dynamic information of discrete ship navigation, wherein the static information comprises ship MMSI and ship intrinsic parameters (such as ship name, ship length, ship type, draft, ship host power information and the like), and the dynamic information comprises ship MMSI, current navigation time point, navigation speed, course, ship coordinate position and the like; according to the correlation between the MMSI of the ship with the dynamic information and the MMSI of the ship with the static information, integrating the static information and the dynamic information of the ship at the current sailing time point to establish a ship information matrix;
wherein, the static information is received within 6min generally; the dynamic information is about 2-12 s according to different navigational speeds; considering that the static information of the same ship passing through the ship emission control area is constant, the ship MMSI of the static information is only required to be matched with the ship MMSI of the dynamic information in a correlation manner so as to integrate the static information and the dynamic information of the ship at the current sailing time point and establish a character-type ship information matrixK ship1、K ship2、…、K Mship，MThe number of ships in the ship emission control area.
S2, identifying the target time interval of the appearance peak value of the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeterT 1、T 2、…、T R Intercepting the corresponding exhaust information sub-matrix in the exhaust information matrix according to the target time intervalD * 1、D * 2、…、D * R ；
Specifically, a man-kenidel trend test method is adopted to identify a target time interval in which the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter has a peak value; the determination of the target time interval, for example, takes the time points corresponding to the half-high peak points on both sides of the peak point where the telemetric concentration appears as the two end points of the target time interval.
S3, polling the target time interval corresponding to each ship tail gas telemeterT j Screening target ships meeting target conditions in the ship information matrix; in addition, information about ships that do not satisfy the target condition is discarded.
Specifically, the target conditions include condition (1) and condition (2):
condition (1):
Figure 229624DEST_PATH_IMAGE012
；tis the current point in time of the voyage of the ship,δthe maximum time for receiving the dynamic information is defaulted to 12 s;j=1，2，…，R；T j，max is a target time intervalT j The maximum value of (a) is,T j，min is a target time intervalT j Is measured.
Condition (2): the distance between the ship coordinate position and the corresponding correlation midpoint position of the ship tail gas telemeter is smaller than the effective detection distance of the corresponding ship tail gas telemeter.
The effective detection distance of the ship tail gas telemeter is related to the wind speed and the wind direction in meteorological information and is generally 500-2000 m;
s4, correlating the target ship and the telemetering concentration of the atmospheric pollutants monitored by at least two ship tail gas telemeters corresponding to the target ship based on the tail gas information submatrix, and correlating the telemetering concentration of the atmospheric pollutants with meteorological information corresponding to a target time interval;
the meteorological information comprises the current environmental wind direction, the wind speed, the temperature, the humidity, the air pressure and the like, and can be acquired through a meteorological analyzer.
In addition, the embodiment correlates the relevant information according to the exhaust information submatrix, thereby reducing the data processing amount and improving the efficiency of correlation matching.
S5, inverting the concentration of the atmospheric pollutants discharged by the tail gas of each ship based on the ship smoke plume Gaussian diffusion model and the nonlinear multi-component analysis algorithm.
Specifically, from the vessel coordinate position of the target vessel received by the AIS receiver: (x 0，y 0) Course of the vehiclewSpeed of flightvExpressed in a geodetic coordinate system, thentPosition coordinates of the voyage point of the target ship at the time (c:)X，Y）：
Figure 393889DEST_PATH_IMAGE013
；
As shown in fig. 2, the position coordinates of the waypoint in the geodetic coordinate system are converted into ship coordinate positions: (x 0，y 0) Is at the origin and has a wind direction angle ofxAxis in the cross-wind directionyPosition coordinates of a navigation point in a wind direction coordinate system established by axes: (X′，Y′）：
Figure 570268DEST_PATH_IMAGE014
Wherein the content of the first and second substances,Lis the distance between the correlation midpoint position of the ship tail gas telemeter and a target ship,Pis the downwind angle;
suppose that the target vessel navigates tot n At the moment, based on a Gaussian plume diffusion model, the concentration of the atmospheric pollutants is calculated as follows:
Figure 783074DEST_PATH_IMAGE015
wherein the content of the first and second substances,
Figure 364228DEST_PATH_IMAGE005
for ships att i The atmospheric pollutants discharged at any moment diffuse tot n Any space point of time
Figure 637078DEST_PATH_IMAGE006
The concentration of the (C) in the (C),Q i for ships att i The total amount of atmospheric pollutants emitted at a time,
Figure 38103DEST_PATH_IMAGE007
are respectively as
Figure 54601DEST_PATH_IMAGE008
The diffusion coefficient in the direction of the light,uwhich is the wind speed,His a point in space
Figure 286999DEST_PATH_IMAGE006
Height from water surface;
the method comprises the steps of obtaining and accumulating the concentrations of the atmospheric pollutants at a plurality of space points in a target time interval, establishing a nonlinear equation with the telemetering concentration of the atmospheric pollutants monitored by a ship tail gas telemeter corresponding to a target ship, and inversing the concentrations of the atmospheric pollutants emitted by the target ships at different moments by utilizing a nonlinear multi-component analysis algorithm.
Wherein the target time interval isfArea monitored by ship exhaust gas telemeter associated with ship exhaust gas telemeter by ship-discharged atmospheric pollutants
Figure 996329DEST_PATH_IMAGE016
The cumulative formula of the concentration of (A) is as follows:
Figure 875862DEST_PATH_IMAGE017
pin the area monitored by ship tail gas telemeterXThe number of the intervals is equally divided in the axial direction, the distance between the main engine and the auxiliary engine of the marine tail gas telemeter is generally 50-200 m and 20cm is used as the interval,pgenerally 250 to 1000;qin the area monitored by ship tail gas telemeterYThe number of the intervals is equally divided in the axial direction, the Y-axis coverage of the ship tail gas telemeter is about 20-60 m, 2cm is taken as an interval,qgenerally 1000 to 3000.
With atmospheric pollutants SO2For example, assume that it goes through step S3After screening, there arek 1A target vessel,k 2The ship tail gas telemeter meets target conditions, and the actually monitored gas concentrations are respectively as follows:
Figure 430471DEST_PATH_IMAGE018
the above-mentionedk 2SO monitored by ship tail gas telemeter2At a concentration ofk 1Discharge of a target vessel to be diffusedk 2And (4) comprehensive monitoring results in a remote area of the ship tail gas remote-measuring instrument.
Suppose thatk 1SO discharged from a target vessel2Respectively in a concentration of
Figure 986217DEST_PATH_IMAGE019
Then the SO is obtained according to the above step S22And (3) a plurality of point data of the interval of the peak value full width at half maximum, and considering that a certain nonlinear relation exists, the nonlinear multi-component analysis model has the following formula:
Figure 866449DEST_PATH_IMAGE020
according to a nonlinear multi-component analysis algorithm, such as a nonlinear least square method, a principal component analysis method and the like, and according to a maximum principle and a typical correlation maximum principle, the SO emitted by each ship monitored in the ship emission control area at different moments is inverted2Concentration of
Figure 242066DEST_PATH_IMAGE021
；
By analogy, the CO emitted by each ship monitored in the ship emission control area can be obtained2Concentration, NO x Concentration, PM concentration.
S6, calculating the sulfur content of the fuel oil of each ship;
specifically, the firstfFuel oil sulfur content of ship
Figure 131525DEST_PATH_IMAGE009
。
Wherein the content of the first and second substances,
Figure 276198DEST_PATH_IMAGE022
is as followsfSO discharged from a ship2The concentration of the active ingredients in the mixture is,
Figure 124069DEST_PATH_IMAGE023
is as followsfCO discharged from ship2And (4) concentration.
Correspondingly, as shown in fig. 3, the terminal for monitoring the concentration of pollutants discharged from multiple ship exhaust according to the embodiment includes:
an acquisition module to acquireRThe telemetering concentration of the atmospheric pollutants monitored by the ship tail gas telemeter establishes a digital tail gas information matrix formed by time points and the telemetering concentration of the atmospheric pollutantsD 1、D 2、…、D R ；
Wherein the atmospheric pollutant includes SO2、CO2、NO x PM, etc., the time for the ship tail gas telemeter to acquire the telemetered concentration is generally 500ms, and then a tail gas information matrixD j The expression of (a) is as follows:
Figure 721403DEST_PATH_IMAGE024
wherein the content of the first and second substances,j=1，2，…，R，t m，j is as followsjThe first of the ship tail gas telemeter monitoringmAt a point in time, the time of the measurement,mis an integer greater than 1;
Figure 680132DEST_PATH_IMAGE025
are respectively the firstjA ship tail gas telemeter is arranged on the secondmCO monitored at each time point2、NO x 、SO2And the remote measured concentration of PM.
The method comprises the steps that the AIS receiver is used for obtaining static information and dynamic information of discrete ship navigation, wherein the static information comprises ship MMSI and ship intrinsic parameters (such as ship name, ship length, ship type, draft, ship host power information and the like), and the dynamic information comprises ship MMSI, current navigation time point, navigation speed, course, ship coordinate position and the like; according to the correlation between the MMSI of the ship with the dynamic information and the MMSI of the ship with the static information, integrating the static information and the dynamic information of the ship at the current sailing time point to establish a ship information matrix;
wherein, the static information is received within 6min generally; the dynamic information is about 2-12 s according to different navigational speeds; considering that the static information of the same ship passing through the ship emission control area is constant, only the static information of the ship MMSI and the dynamic information of the ship MMSI are matched in a correlation mode to integrate the static information and the dynamic information of the ship at the current sailing time point and establish a character-type ship information matrixK ship1、K ship2、…、K Mship，MThe number of ships in the ship emission control area;
the identification module is used for identifying a target time interval of the peak value of the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeterT 1、T 2、…、T R Intercepting the corresponding exhaust information sub-matrix in the exhaust information matrix according to the target time intervalD * 1、D * 2、…、D * R ；
Specifically, a man-kenidel trend test method is adopted to identify a target time interval in which the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter has a peak value; the determination of the target time interval, for example, takes the time points corresponding to the half-high peak points on both sides of the peak point where the telemetric concentration appears as the two end points of the target time interval.
A screening module for polling the corresponding target time interval of each ship tail gas telemeterT j Screening target ships meeting target conditions in the ship information matrix; in addition, information about ships that do not satisfy the target condition is discarded.
Specifically, the target conditions include condition (1) and condition (2):
condition (1):
Figure 410803DEST_PATH_IMAGE012
；tis the current point in time of the voyage of the ship,δthe maximum time for receiving the dynamic information is defaulted to 12 s;j=1，2，…，R；T j，max is a target time intervalT j The maximum value of (a) is,T j，min is a target time intervalT j Is measured.
Condition (2): the distance between the ship coordinate position and the corresponding correlation midpoint position of the ship tail gas telemeter is smaller than the effective detection distance of the corresponding ship tail gas telemeter.
The effective detection distance of the ship tail gas telemeter is related to the wind speed and the wind direction in meteorological information and is generally 500-2000 m;
the correlation module is used for correlating the target ship with the telemetering concentration of the atmospheric pollutants monitored by at least two ship tail gas telemeters corresponding to the target ship based on the tail gas information submatrix and is also used for correlating the telemetering concentration of the atmospheric pollutants with meteorological information corresponding to a target time interval;
the meteorological information comprises the current environmental wind direction, the wind speed, the temperature, the humidity and the air pressure, and can be acquired through a meteorological analyzer.
And the inversion module is used for inverting the concentration of the atmospheric pollutants discharged by the tail gas of each ship based on the ship smoke plume Gaussian diffusion model and the nonlinear multi-component analysis algorithm.
Specifically, from the vessel coordinate position of the target vessel received by the AIS receiver: (x 0，y 0) Course of the vehiclewSpeed of flightvExpressed in a geodetic coordinate system, thentPosition coordinates of the voyage point of the target ship at the time (c:)X，Y）：
Figure 429575DEST_PATH_IMAGE013
；
Converting position coordinates of a navigation point in a geodetic coordinate system into position coordinates of a shipCoordinate position (x 0，y 0) Is at the origin and has a wind direction angle ofxAxis in the cross-wind directionyPosition coordinates of a navigation point in a wind direction coordinate system established by axes: (X′，Y′）：
Figure 310943DEST_PATH_IMAGE026
Wherein the content of the first and second substances,Lis the distance between the correlation midpoint position of the ship tail gas telemeter and a target ship,Pis the downwind angle;
suppose that the target vessel navigates tot n At the moment, based on a Gaussian plume diffusion model, the concentration of the atmospheric pollutants is calculated as follows:
Figure 11046DEST_PATH_IMAGE027
wherein the content of the first and second substances,
Figure 927049DEST_PATH_IMAGE005
for ships att i The atmospheric pollutants discharged at any moment diffuse tot n Any space point of time
Figure 788826DEST_PATH_IMAGE006
The concentration of the (C) in the (C),Q i for ships att i The total amount of atmospheric pollutants emitted at a time,
Figure 626332DEST_PATH_IMAGE007
are respectively as
Figure 926863DEST_PATH_IMAGE008
The diffusion coefficient in the direction of the light,uwhich is the wind speed,His a point in space
Figure 900635DEST_PATH_IMAGE006
Height from water surface;
the method comprises the steps of obtaining and accumulating the concentrations of the atmospheric pollutants at a plurality of space points in a target time interval, establishing a nonlinear equation with the telemetering concentration of the atmospheric pollutants monitored by a ship tail gas telemeter corresponding to a target ship, and inversing the concentrations of the atmospheric pollutants emitted by the target ships at different moments by utilizing a nonlinear multi-component analysis algorithm.
Wherein the target time interval isfArea monitored by ship exhaust gas telemeter associated with ship exhaust gas telemeter by ship-discharged atmospheric pollutants
Figure 198893DEST_PATH_IMAGE028
The cumulative formula of the concentration of (A) is as follows:
Figure 526625DEST_PATH_IMAGE029
pin the area monitored by ship tail gas telemeterXThe number of the intervals is equally divided in the axial direction, the distance between the main engine and the auxiliary engine of the marine tail gas telemeter is generally 50-200 m and 20cm is used as the interval,pgenerally 250 to 1000;qin the area monitored by ship tail gas telemeterYThe number of the intervals is equally divided in the axial direction, the Y-axis coverage of the ship tail gas telemeter is about 20-60 m, 2cm is taken as an interval,qgenerally 1000 to 3000.
With atmospheric pollutants SO2For example, assume that after the filtering of step S3, there isk 1A target vessel,k 2The ship tail gas telemeter meets target conditions, and the actually monitored gas concentrations are respectively as follows:
Figure 834109DEST_PATH_IMAGE030
the above-mentionedk 2SO monitored by ship tail gas telemeter2At a concentration ofk 1Discharge of a target vessel to be diffusedk 2And (4) comprehensive monitoring results in a remote area of the ship tail gas remote-measuring instrument.
Suppose thatk 1SO discharged from a target vessel2Respectively in a concentration of
Figure 459125DEST_PATH_IMAGE031
Then the SO is obtained according to the above step S22And (3) a plurality of point data of the interval of the peak value full width at half maximum, and considering that a certain nonlinear relation exists, the nonlinear multi-component analysis model has the following formula:
Figure 193863DEST_PATH_IMAGE020
according to a nonlinear multi-component analysis algorithm, such as a nonlinear least square method, a principal component analysis method and the like, and according to a maximum principle and a typical correlation maximum principle, the SO emitted by each ship monitored in the ship emission control area at different moments is inverted2Concentration;
by analogy, the CO emitted by each ship monitored in the ship emission control area can be obtained2Concentration, NO x Concentration, PM concentration.
The inversion module is also used for calculating the fuel oil sulfur content of each ship;
specifically, the firstfFuel oil sulfur content of ship
Figure 474803DEST_PATH_IMAGE009
。
Wherein the content of the first and second substances,
Figure 382716DEST_PATH_IMAGE022
is as followsfSO discharged from a ship2The concentration of the active ingredients in the mixture is,
Figure 331081DEST_PATH_IMAGE023
is as followsfCO discharged from ship2And (4) concentration.
As shown in fig. 4 and 5, the monitoring system for the concentration of pollutants discharged from multiple ships according to the present embodiment is disposed in a port area, and includes:
Rthe ship tail gas telemeter 1 is arranged in a ship emission control area; the remote measuring areas of the ship tail gas telemeters are mutually independent;
the monitoring terminal 2 of the embodiment is in communication connection with each ship tail gas telemeter;
the AIS receiver 3 is in communication connection with the AIS transmitter and the monitoring terminal of each ship 0 respectively; the AIS receiver is used for receiving static information and dynamic information of navigation of each ship in the ship emission control area;
the meteorological analyzer 4 is in communication connection with the monitoring terminal; the meteorological analyzer is used for acquiring meteorological information including current environmental wind direction, wind speed, temperature, humidity, air pressure and the like.
In addition, the monitoring terminal can be in communication connection with the cloud server, and remote monitoring is achieved.
Example 2:
the method for monitoring the concentration of pollutants discharged by multiple ship tail gases in the embodiment is different from the method in embodiment 1 in that:
the target condition only retains the condition (2) and meets the requirements of different applications; namely, the target conditions of the present embodiment are:
the distance between the ship coordinate position and the corresponding correlation midpoint position of the ship tail gas telemeter is smaller than the effective detection distance of the corresponding ship tail gas telemeter;
the effective detection distance of the ship tail gas telemeter is related to the wind speed and the wind direction in meteorological information, and is generally 500-2000 m.
Other steps can be referred to example 1.
Accordingly, the monitoring terminal and the monitoring system are adjusted accordingly.
Example 3:
the method for monitoring the concentration of pollutants discharged by multiple ship tail gases in the embodiment is different from the method in embodiment 1 in that:
the monitoring method is applied to inland river areas, and the layout of each ship tail gas telemeter is shown in figure 6;
other steps can be referred to example 1.
The foregoing has outlined rather broadly the preferred embodiments and principles of the present invention and it will be appreciated that those skilled in the art may devise variations of the present invention that are within the spirit and scope of the appended claims.

Claims (6)

1. A monitoring method for the concentration of pollutants discharged by multi-ship tail gas is characterized in that a plurality of ship tail gas telemeters are distributed in a ship discharge control area, and the telemetering areas of the ship tail gas telemeters are mutually independent;
the monitoring method comprises the following steps:
s1, acquiring the telemetered concentration of the atmospheric pollutants monitored by each ship tail gas telemeter, and establishing a tail gas information matrix formed by the time point and the telemetered concentration of the atmospheric pollutants;
acquiring static information and dynamic information of a ship, wherein the static information comprises a ship MMSI and inherent parameters of the ship, and the dynamic information comprises the ship MMSI, a current navigation time point, a navigation speed, a course and a ship coordinate position; according to the MMSI of the ship and the static information and the dynamic information of the ship integrated at the current sailing time point, establishing a ship information matrix;
s2, identifying a target time interval of the peak value of the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter, and intercepting in a tail gas information matrix according to the target time interval to obtain a corresponding tail gas information sub-matrix;
s3, polling target time intervals corresponding to the ship tail gas telemeters, and screening target ships meeting target conditions in a ship information matrix;
wherein the target conditions include condition (1) and condition (2):
condition (1):
Figure 505108DEST_PATH_IMAGE001
；tis the current point in time of the voyage of the ship,δfor maximum time of reception of dynamic information, TmaxIs the maximum value of the target time interval, TminIs the minimum value of the target time interval;
condition (2): the distance between the ship coordinate position and the corresponding correlation midpoint position of the ship tail gas telemeter is smaller than the effective detection distance of the corresponding ship tail gas telemeter;
alternatively, the target conditions are:
the distance between the ship coordinate position and the corresponding correlation midpoint position of the ship tail gas telemeter is smaller than the effective detection distance of the corresponding ship tail gas telemeter;
s4, correlating the target ship and the telemetering concentration of the atmospheric pollutants monitored by at least two ship tail gas telemeters corresponding to the target ship based on the tail gas information submatrix, and correlating the telemetering concentration of the atmospheric pollutants with meteorological information corresponding to a target time interval; the meteorological information comprises the current environmental wind direction, wind speed, temperature, humidity and air pressure;
s5, inverting the concentration of the atmospheric pollutants discharged by the tail gas of each ship based on the ship smoke plume Gaussian diffusion model and the nonlinear multi-component analysis algorithm, wherein the inversion comprises the following steps:
based on the received vessel coordinate position of the target vessel: (x 0，y 0) Course of the vehiclewSpeed of flightvExpressed in a geodetic coordinate system, thentPosition coordinates of the voyage point of the target ship at the time (c:)X，Y）：
Figure 848060DEST_PATH_IMAGE002
Converting the position coordinates of the navigation point in the geodetic coordinate system into the ship coordinate position (x 0，y 0) Is at the origin and has a wind direction angle ofxAxis in the cross-wind directionyPosition coordinates of a navigation point in a wind direction coordinate system established by axes: (X′，Y′）：
Figure 126594DEST_PATH_IMAGE003
Wherein the content of the first and second substances,Lis the distance between the correlation midpoint position of the ship tail gas telemeter and a target ship,Pis the downwind angle;
suppose that the target vessel navigates tot n At the moment, based on a Gaussian plume diffusion model, the concentration of the atmospheric pollutants is calculated as follows:
Figure 822149DEST_PATH_IMAGE004
wherein the content of the first and second substances,
Figure 377895DEST_PATH_IMAGE005
for ships att i The atmospheric pollutants discharged at any moment diffuse tot n Any space point of time
Figure 913919DEST_PATH_IMAGE006
The concentration of the (C) in the (C),Q i for ships att i The total amount of atmospheric pollutants emitted at a time,
Figure 820695DEST_PATH_IMAGE007
are respectively as
Figure 84055DEST_PATH_IMAGE008
The diffusion coefficient in the direction of the light,uwhich is the wind speed,His a point in space
Figure 25466DEST_PATH_IMAGE009
Height from water surface;
the method comprises the steps of obtaining and accumulating the concentrations of the atmospheric pollutants at a plurality of space points in a target time interval, establishing a nonlinear equation with the telemetering concentration of the atmospheric pollutants monitored by a ship tail gas telemeter corresponding to a target ship, and inversing the concentrations of the atmospheric pollutants emitted by the target ships at different moments by utilizing a nonlinear multi-component analysis algorithm.
2. The method of monitoring of claim 1, wherein the nonlinear multicomponent analysis algorithm is a nonlinear least squares method or a principal component analysis method.
3. The method of monitoring of claim 1, wherein the atmospheric pollutant is a sulfur oxide, a nitrogen oxide, a carbon oxide, or a particulate matter.
4. The monitoring method according to claim 3, wherein after the step S5, the method further comprises the steps of:
s6, calculating the sulfur content of the fuel oil of the ship
Figure 466812DEST_PATH_IMAGE010
；
Wherein the content of the first and second substances,C S is the concentration of sulfur oxides in the vessel,C C is the concentration of carbon oxides of the vessel.
5. The monitoring method of claim 1 is applied to a monitoring terminal for the concentration of pollutants discharged from multiple ship exhaust gases, and the monitoring terminal comprises:
the acquisition module is used for acquiring the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter and establishing a tail gas information matrix formed by time points and the telemetering concentration of the atmospheric pollutants; the system is also used for acquiring static information and dynamic information of the ship, wherein the static information comprises the MMSI of the ship and inherent parameters of the ship, and the dynamic information comprises the MMSI of the ship, the current time point of navigation, the speed, the course and the coordinate position of the ship; according to the MMSI of the ship and the static information and the dynamic information of the ship integrated at the current sailing time point, establishing a ship information matrix;
the identification module is used for identifying a target time interval of the peak value of the telemetering concentration of the atmospheric pollutants monitored by each ship tail gas telemeter, and intercepting the target time interval in a tail gas information matrix to obtain a corresponding tail gas information sub-matrix;
the screening module is used for polling the target time intervals corresponding to the ship tail gas telemeters and screening target ships meeting target conditions in the ship information matrix;
the correlation module is used for correlating the target ship with the telemetering concentration of the atmospheric pollutants monitored by at least two ship tail gas telemeters corresponding to the target ship based on the tail gas information submatrix and is also used for correlating the telemetering concentration of the atmospheric pollutants with meteorological information corresponding to a target time interval;
and the inversion module is used for inverting the concentration of the atmospheric pollutants discharged by the tail gas of each ship based on the ship smoke plume Gaussian diffusion model and the nonlinear multi-component analysis algorithm.
6. A monitoring system for concentration of pollutants discharged from tail gas of multiple ships, comprising:
the ship tail gas telemeters are distributed in the ship emission control area; the remote measuring areas of the ship tail gas telemeters are mutually independent;
the monitoring terminal of claim 5, in communication with each marine exhaust telemeter;
the AIS receiver is in communication connection with an AIS transmitter and a monitoring terminal of the ship respectively;
and the meteorological analyzer is in communication connection with the monitoring terminal.
CN202110982164.9A 2021-08-25 2021-08-25 Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas Active CN113419037B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110982164.9A CN113419037B (en) 2021-08-25 2021-08-25 Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110982164.9A CN113419037B (en) 2021-08-25 2021-08-25 Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas

Publications (2)

Publication Number Publication Date
CN113419037A CN113419037A (en) 2021-09-21
CN113419037B true CN113419037B (en) 2021-11-02

Family

ID=77719408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110982164.9A Active CN113419037B (en) 2021-08-25 2021-08-25 Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas

Country Status (1)

Country Link
CN (1) CN113419037B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114354553B (en) * 2021-12-14 2023-08-01 杭州春来科技有限公司 Sniffing monitoring method and system for sulfur content of ship fuel oil
CN113933258B (en) * 2021-12-21 2022-04-12 杭州春来科技有限公司 VOCs pollutant tracing method, terminal and system based on navigation monitoring
CN114218231B (en) * 2022-02-21 2022-05-17 杭州春来科技有限公司 Ship tail gas monitoring data processing method and system and computer readable storage medium
CN115032338B (en) * 2022-05-26 2023-03-21 武汉理工大学 Port ship atmospheric pollutant emission monitoring site location method and system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109060014A (en) * 2018-07-16 2018-12-21 交通运输部水运科学研究所 A kind of marine exhaust discharge compliance determines and violation ship identification method and system
CN111256755A (en) * 2020-01-20 2020-06-09 武汉理工大学 Mobile ship exhaust emission tracing device and method
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

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9719972B2 (en) * 2015-03-31 2017-08-01 International Business Machines Corporation System and method for air-pollutant source-localization using parked motor vehicles
US11731747B2 (en) * 2019-06-19 2023-08-22 Schlumberger Technology Corporation Marine emissions detection

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109060014A (en) * 2018-07-16 2018-12-21 交通运输部水运科学研究所 A kind of marine exhaust discharge compliance determines and violation ship identification method and system
CN111256755A (en) * 2020-01-20 2020-06-09 武汉理工大学 Mobile ship exhaust emission tracing device and method
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

Also Published As

Publication number Publication date
CN113419037A (en) 2021-09-21

Similar Documents

Publication Publication Date Title
CN113419037B (en) Method, terminal and system for monitoring concentration of pollutants discharged by multi-ship tail gas
CN109060014B (en) Ship exhaust emission compliance judgment and illegal ship identification method
CN107389880B (en) Ship tail gas mobile monitoring method
CN107589100B (en) A kind of bunker oil sulfur content sniff estimation algorithm
Ridley et al. Distributions of NO, NOx, NOy, and O3 to 12 km altitude during the summer monsoon season over New Mexico
US20210131964A1 (en) Remote sensing-based detection system and method for gaseous pollutant from diesel vehicle exhaust
CN109978335B (en) Ship fuel compliance judgment and violation identification method
CN109063910A (en) A kind of Pollution From Ships object discharge method of real-time based on big data
CN111811572A (en) Ship exhaust emission real-time monitoring method based on big data
CN111238570B (en) Ship exhaust pollutant detection system and method
CN111289040A (en) Navigation monitoring system
CN113125683B (en) Marine fuel oil sulfur content mobile monitoring device and method
Zhou et al. Monitoring the compliance of sailing ships with fuel sulfur content regulations using unmanned aerial vehicle (UAV) measurements of ship emissions in open water
CN212621763U (en) Boats and ships plume SO based on unmanned aerial vehicle2Discharge collecting and detecting device
CN215263401U (en) Portable emission test system
CN113960252A (en) Method and system for monitoring tail gas of navigation channel ship and data processing equipment
Walden et al. Measurement report: Characterization of uncertainties in fluxes and fuel sulfur content from ship emissions in the Baltic Sea
CN113052042A (en) Pollutant emission source monitoring device and method
Wu et al. Development of a spectrum-based ship fuel sulfur content real-time evaluation method
Menaga et al. Air quality monitoring system using vehicles based on the IoT
Pszenny et al. Overview of the 1988 GCE/CASE/WATOX studies of biogeochemical cycles in the North Atlantic region
CN111855901A (en) Ship emission testing method based on power emission factor
Peng et al. System for identifying ships using high-sulfur-content fuel oil in port waters
Zhou et al. Monitoring compliance with fuel sulfur content regulations of sailing ships by unmanned aerial vehicle (UAV) measurements of ship emissions in open water
CN114236101B (en) Method for calculating sulfur content of ship fuel without measuring carbon concentration

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