CN117112969B - Vehicle pollutant emission list accounting method - Google Patents
Vehicle pollutant emission list accounting method Download PDFInfo
- Publication number
- CN117112969B CN117112969B CN202311386470.1A CN202311386470A CN117112969B CN 117112969 B CN117112969 B CN 117112969B CN 202311386470 A CN202311386470 A CN 202311386470A CN 117112969 B CN117112969 B CN 117112969B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- emission
- target
- pollutant
- load
- 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
- 239000003344 environmental pollutant Substances 0.000 title claims abstract description 84
- 231100000719 pollutant Toxicity 0.000 title claims abstract description 84
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000012937 correction Methods 0.000 claims abstract description 17
- 230000001133 acceleration Effects 0.000 claims description 22
- 239000002826 coolant Substances 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 101150040844 Bin1 gene Proteins 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000005284 excitation Effects 0.000 claims description 2
- 101100540662 Glycine max VSPA gene Proteins 0.000 claims 3
- 238000011156 evaluation Methods 0.000 abstract description 3
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 36
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 description 30
- 229910002091 carbon monoxide Inorganic materials 0.000 description 30
- 229910052757 nitrogen Inorganic materials 0.000 description 18
- CURLTUGMZLYLDI-UHFFFAOYSA-N Carbon dioxide Chemical compound O=C=O CURLTUGMZLYLDI-UHFFFAOYSA-N 0.000 description 14
- 229910002092 carbon dioxide Inorganic materials 0.000 description 7
- 239000001569 carbon dioxide Substances 0.000 description 7
- 238000010276 construction Methods 0.000 description 5
- 239000000356 contaminant Substances 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Data Mining & Analysis (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Physics (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Quality & Reliability (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
The invention discloses a vehicle pollutant emission list accounting method, which is based on one of two preset mutually independent Bin interval dividing methods, and establishes a multi-scale Bin interval of which the target vehicle covers cold start, low load and medium and high load operation stages; based on the reference emission factor, the Bin interval mileage ratio and the correction coefficient, obtaining an emission factor capable of representing emission characteristics of a target area, and constructing a multi-scale emission factor set; and constructing and obtaining a multi-scale vehicle pollutant emission list in the target area according to the multi-scale emission factor set and all vehicle driving mileage of the target area. The invention solves the problems that the existing emission factors and the list evaluation results are macroscopic and cannot characterize the multi-scale traffic emission characteristics in the area, and provides a new method for quantifying and evaluating the emission of the vehicle pollutants in the traffic field.
Description
Technical Field
The invention relates to the technical field of vehicle pollutant emission accounting, in particular to a vehicle pollutant emission list accounting method.
Background
The microscopic operating mode (Bin) interval of the vehicle (hereinafter referred to as the vehicle Bin interval or Bin interval) is the basis for accounting for emission factors and pollutant emission schedules. In the existing research, the Bin intervals of the vehicle are mainly divided based on specific power (VSP), and the emission factors and lists obtained based on the accounting can characterize the emission characteristics of the vehicle under different powers. However, there are certain limitations to the corresponding emissions factors and inventory due to the manner in which specific power (VSP) is calculated. Such as: specific Power (VSP) is calculated through vehicle speed, acceleration and gradient, and engine torque, rotational speed and efficiency are absent from the calculation, so that the emission factors and lists cannot directly consider the influence of the running characteristics of the engine on the emission of vehicle pollutants. The specific power (VSP) only characterizes the power variation in the vehicle operating state and cannot intuitively distinguish the vehicle operating phases of cold start, low load, etc., so that the emission factors and the list cannot fully reflect the emission intensity of the vehicle pollutants in the high emission intervals of cold start, low load, etc. The emission characteristics of different vehicle types are different, for example, the light vehicle needs to consider the influence of the engine on emission, the engine of the heavy vehicle is in a low-load running state for a long time, and the important quantification of the emission characteristics of the engine of the heavy vehicle in low-load running is particularly important, so that the construction of a targeted emission factor and a list is important. However, since specific power (VSP) is a relatively macroscopic technical indicator, it is difficult to specifically characterize emissions factors and lists for specific vehicle types. Therefore, the evaluation results of emission factors and lists obtained based on VSP-Bin interval accounting are macroscopic and cannot characterize the multi-scale traffic emission characteristics in the area.
Disclosure of Invention
The invention aims at solving the technical defects existing in the prior art, provides a vehicle pollutant emission list accounting method, solves the problems that the existing emission factors and list evaluation results are macroscopic and cannot characterize the multi-scale traffic emission characteristics in the area, and provides a new method for quantifying and evaluating the vehicle pollutant emission in the traffic field.
A vehicle pollutant emission inventory accounting method comprising the steps of:
s1, dividing the vehicle operation stage into cold start, low load and medium and high load;
the cold start operation phase is defined as a vehicle operation phase in a preset time period after the engine is started when the temperature of the engine coolant does not exceed a preset temperature threshold value, and is set as Bin1; the low-load operation stage is a vehicle operation stage in which the engine load is lower than a full-load preset percentage value in a non-cold start operation stage; the medium-high load operation phase is a vehicle operation phase except for cold start and low load;
s2, when a pollutant emission list of a target vehicle is calculated, selecting a vehicle Bin interval from two mutually independent Bin interval dividing methods according to the model specification of the target vehicle in the low-load and medium-high-load operation stages, and establishing a multi-scale Bin interval of the target vehicle;
for a first target vehicle, constructing a multi-scale Bin interval by adopting a Bin interval dividing method I; the Bin interval dividing method comprises the following steps:
for a low-load operation stage, carrying out Bin interval division based on two variables of engine load and vehicle speed;
for the medium-high load operation stage, dividing a Bin interval is carried out based on a vehicle speed per acceleration VSPA and two variables of the vehicle speed, and the vehicle speed per acceleration VSPA per second reflects the driving excitation degree of the vehicle;
for a second target vehicle, constructing a multi-scale Bin interval by adopting a Bin interval dividing method II; the second method for dividing the Bin section adopts a mutual information method, fuses vehicle characteristic parameters including vehicle speed, acceleration, gradient, rotating speed and torque, constructs a characteristic fitting parameter FFP, and simultaneously, develops the division of the vehicle Bin section based on the characteristic fitting parameter FFP and two variables of engine load, and comprises the following steps:
according to each parameter change range [ min, max ] set in advance, carrying out normalization processing on a vehicle characteristic parameter F and a target pollutant p to be calculated, and mapping data to a [0,1] interval, wherein the vehicle characteristic parameter F comprises vehicle speed, acceleration, gradient, rotating speed and torque;
calculating the condition mutual information h between the vehicle characteristic parameter F and the target pollutant p one by one F,p As a coefficient of the vehicle characteristic parameter F, wherein the condition mutual information h F,p Representing the correlation coefficient of the vehicle characteristic parameter F and the target pollutant p;
for all the vehicle characteristic parameters F, calculating mutual information I between every two parameters a,b Wherein, the method comprises the steps of, wherein,calculate the sum of all mutual informationAs a constant term coefficient, where mutual information I a,b Representing a redundant relationship between the vehicle characteristic parameters F;
obtaining characteristic fitting parameters FFP for target pollutant p p The formula of (2) is shown below;
；
s3, acquiring actual emission data of a target vehicle, and calculating a vehicle reference emission factor of the target vehicle according to the constructed multi-scale Bin interval;
，
in the method, in the process of the invention,the standard emission factor of the target pollutant p in g/km under the jth Bin section of the stage operation stage is represented; the stage represents three vehicle operation stages of cold start, low load, medium load and high load;the emission quantity of the target pollutant p in the jth Bin section of the stage operation stage is expressed in g, and the emission quantity is obtained by summing the emission rates of the pollutants p in the Bin sections;the vehicle driving mileage under the jth Bin section in the stage operation stage is expressed, and the unit is km;
s4, calculating a cold start emission factor, a low load emission factor, a medium and high load emission factor and a comprehensive emission factor based on a vehicle reference emission factor and a Bin interval driving mileage ratio, and constructing a multi-scale emission factor set in a target area;
，
，
in the method, in the process of the invention,the emission factor of the target pollutant p in the stage operation stage is expressed in g/km;the comprehensive emission factor of the target pollutant p is expressed in g/km;representing the driving mileage duty ratio of the jth Bin interval;the correction coefficient representing the reference emission factor is equal to the product of the average running speed correction factor, the weather information correction factor and the oil standard correction factor in the target area;
s5, acquiring the driving mileage of the bicycle in the target area, calculating the product of the multi-scale emission factor set and the driving mileage of the bicycle in the target area, and acquiring the pollutant emission intensity and the comprehensive emission intensity of the bicycle in the cold start, low load and medium and high load operation stages in the target area;
，
，
in the method, in the process of the invention,the emission intensity of the target pollutant p in the stage operation stage is expressed as g;the comprehensive emission intensity of the target pollutant p is expressed in g;indicating the mileage of a bicycle.
S6, obtaining the driving mileage of all vehicles in the target area, repeating the step S5, calculating and obtaining the pollutant emission intensity and the comprehensive emission intensity of all vehicles in the target area in the cold start, low load, medium and high load operation stages, and constructing a multi-scale vehicle pollutant emission list in the target area.
Wherein the preset temperature threshold is 70 ℃, and the preset time period is 300s.
Wherein the preset percentage value is 40%.
Wherein, the vehicle per acceleration VSPA is obtained through the following calculation;
wherein v represents the vehicle speed; a represents acceleration.
The first target vehicle comprises a large-sized passenger car, a bus, a medium-sized cargo car and a heavy-duty cargo car; the second target vehicle comprises a mini-type passenger car, a taxi, a medium-sized passenger car, a mini-type cargo car, a light-duty cargo car and a motorcycle.
According to the invention, two mutually independent Bin interval dividing methods are adopted, a multi-scale Bin interval covering cold start, low load and medium and high load operation stages is established, a multi-scale emission factor set capable of reflecting high emission characteristics of cold start, low load and the like and capable of directly representing engine operation characteristics is constructed based on the multi-scale Bin interval, a multi-scale vehicle pollutant emission list in a target area is further calculated and obtained, and a novel method is provided for quantifying and evaluating vehicle pollutant emission in the traffic field.
Drawings
Fig. 1 is a flow chart of a vehicle pollutant emission inventory accounting method according to an embodiment of the present invention.
Fig. 2 is a table of multi-scale Bin section division results of a heavy-duty cargo vehicle as a target vehicle according to a Bin section division method.
FIGS. 3, 4 and 5 are diagrams of CO and CO at different Bin intervals for a heavy truck contaminant as a target vehicle according to an embodiment of the present invention 2 、No x A histogram of the baseline emission factor.
Fig. 6 is a bar chart of mileage duty ratio for each operation stage of a heavy truck as a target vehicle according to an embodiment of the present invention.
Fig. 7 is a schematic view showing a list of taillights of a heavy-duty truck as a target vehicle according to an embodiment of the present invention.
Fig. 8 is a schematic illustration of a multi-scale exhaust emissions inventory for a plurality of heavy trucks in a target area in accordance with an embodiment of the present invention.
FIG. 9 is a diagram of a Bin interval constructed based on a Bin interval partitioning method II according to an embodiment of the invention.
Detailed Description
The invention is described in further detail below with reference to the drawings and the specific examples. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
According to the multi-scale emission factor set construction and vehicle pollutant emission list accounting method, a multi-scale Bin interval which is covered by a target vehicle and has cold start, low load and medium and high load operation stages is established according to one of two mutually independent Bin interval dividing methods; based on the reference emission factor, the Bin interval mileage ratio and the correction coefficient, obtaining an emission factor capable of representing emission characteristics of a target area, and constructing a multi-scale emission factor set; and constructing and obtaining a multi-scale vehicle pollutant emission list in the target area according to the multi-scale emission factor set and all vehicle driving mileage of the target area.
The two mutually independent Bin interval dividing methods are respectively as follows:
the method comprises the following steps: carrying out Bin interval division from engine load, VSPA and vehicle speed variables;
the second method is as follows: and constructing FFP by adopting a mutual information method, carrying out Bin interval division based on the FFP and the engine load, and then calculating a vehicle reference emission factor based on the multi-scale Bin interval and the acquired actual emission data of the vehicle.
As shown in fig. 1, the method for checking a vehicle pollutant emission list according to the embodiment of the invention comprises the following steps:
s1, dividing the vehicle operation stage into three types of cold start, low load and medium and high load;
wherein, the cold start operation phase is defined as a vehicle operation phase in 300s after the engine is started when the temperature of the engine coolant is not more than 70 ℃, and is set as Bin1;
the low load operating phase is defined as a vehicle operating phase in which the engine load is less than 40% of full load during the non-cold start operating phase;
the medium-high load operation phase is defined as a vehicle operation phase other than cold start and low load;
s2, when a pollutant emission list of a target vehicle is calculated, selecting a vehicle Bin interval from two mutually independent Bin interval dividing methods according to the model specification of the target vehicle in the low-load and medium-high-load operation stages, and establishing a multi-scale Bin interval of the target vehicle;
according to the embodiment of the invention, two mutually independent Bin interval dividing methods are provided for dividing the Bin interval of the vehicle according to different vehicle types and specifications aiming at low-load and medium-high-load operation stages; selecting a specific Bin interval dividing method according to the model specification of a target vehicle when a vehicle pollutant emission list is calculated subsequently;
wherein, according to the technical guideline (trial run) of the road motor vehicle emission list, the specific model specification includes: mini-size passenger car, small-size passenger car, taxi, medium-size passenger car, large-size passenger car, bus, mini-size cargo car, light cargo car, medium-size cargo car, heavy cargo car, motorcycle, etc.
In the following, a first target vehicle, such as a heavy truck, is taken as an example to describe the construction of a vehicle Bin section:
the invention is trueIn the application case, the actual emission data of hundreds of heavy trucks meeting the national six-emission standard from the remote monitoring platform and the hub test are collected and arranged, wherein the data comprise vehicle operation parameters such as vehicle speed, acceleration, gradient, torque, rotating speed, engine coolant temperature, engine load and the like, and carbon monoxide (CO) and carbon dioxide (CO) 2 ) Oxides of Nitrogen (NO) x ) And the pollutant emission rate parameters, the data frequency is 1Hz, are used as the data basis for accounting the pollutant emission list of the vehicle.
For heavy trucks, a Bin interval dividing method I is adopted to construct corresponding vehicle multi-scale Bin intervals. The specific data processing mode of the Bin interval dividing method I is as follows:
for the low-load operation phase, the Bin interval division is performed based on two variables of the engine load and the vehicle speed. The engine load is divided into 4 intervals according to [0%, 10%), [10%, 20%), [20%, 30%), [30%, 40%; the speed of the vehicle is divided into 3 sections according to a low speed [0,60 ] km/h, a medium speed [60,90 ] km/h and a high speed [90, ++ infinity) km/h. Thus 12 intervals are obtained for the low load operation phase in total.
For the middle-high load operation stage, the Bin interval is divided based on two variables such as the vehicle speed per acceleration VSPA (Vehicle speed per acceleration) and the vehicle speed.
Wherein VSPA is obtained by calculation of the following formula.
，
Wherein v represents the vehicle speed; a represents acceleration. The second-by-second vehicle speed per acceleration VSPA reflects the degree of driving excitement of the vehicle.
Wherein VSPA is divided into 9 intervals according to [ - ≡3), [ -3, -1), [ -1, 0), [0, 1), [1, 3), [3, 5), [5, 10), [15, + fact +; the speed of the vehicle is divided into 3 sections according to a low speed [0,60 ] km/h, a medium speed [60,90 ] km/h and a high speed [90, ++ infinity) km/h. Thus, 27 intervals are obtained for the middle-high load operation stage in total.
The multi-scale Bin section division result of the heavy cargo vehicle obtained by adopting the Bin section division method I is shown in fig. 2, and the multi-scale Bin section constructed and obtained to cover cold start, low load and medium and high load running states should comprise 40 Bin sections. Fig. 2 shows only 35 Bin sections, because in the actual driving process, there is a phenomenon that the frequency of occurrence of the vehicle running situation under certain Bin sections is extremely low, so that the Bin sections with the mileage ratio of less than 0.02% are additionally processed, and are combined with the Bin sections with the similar vehicle running situation.
When a second target vehicle (such as a light truck) is subjected to multi-scale Bin interval division, a corresponding vehicle multi-scale Bin interval can be constructed by adopting a Bin interval division method II, and the specific implementation flow is as follows:
firstly, acquiring actual road emission data of a light truck, wherein the data comprise vehicle operation parameters such as vehicle speed, acceleration, gradient, torque, rotating speed, engine coolant temperature, engine load and the like, and carbon monoxide (CO) and carbon dioxide (CO) 2 ) Oxides of Nitrogen (NO) x ) And the pollutant emission rate parameters, the data frequency is 1Hz, are used as the data basis for accounting the pollutant emission list of the vehicle.
Secondly, according to a calculation theory of a Bin interval dividing method II, a mutual information method is adopted, vehicle characteristic parameters of vehicle speed, acceleration, gradient, rotation speed and torque are fused, a characteristic fitting parameter FFP is constructed, meanwhile, based on the characteristic fitting parameter FFP and two variables of engine load, the division of the Bin interval of the vehicle is carried out, and the method comprises the following steps:
according to each parameter change range [ min, max ] set in advance, carrying out normalization processing on a vehicle characteristic parameter F and a target pollutant p to be calculated, and mapping data to a [0,1] interval, wherein the vehicle characteristic parameter F comprises vehicle speed, acceleration, gradient, rotating speed and torque;
based on the normalized vehicle characteristic parameters and target pollutant data, calculating the condition mutual information h between the vehicle characteristic parameters F and the target pollutant p one by one F,p As a coefficient of the vehicle characteristic parameter F, wherein the condition mutual information h F,p Representing the characteristic parameter F of the vehicle and the target pollutant pCorrelation coefficients;
for all the vehicle characteristic parameters F, calculating mutual information I between every two parameters a,b Wherein, the method comprises the steps of, wherein,calculate the sum of all mutual informationAs a constant term coefficient, where mutual information I a,b Representing a redundant relationship between the vehicle characteristic parameters F;
obtaining characteristic fitting parameters FFP for target pollutant p p The formula of (2) is shown below;
。
for target pollutants carbon monoxide (CO), carbon dioxide (CO) 2 ) Oxides of Nitrogen (NO) x ) Respectively obtaining characteristic fitting parameters FFP p The formula of (2) is as follows;
CO：，
CO 2 ：，
NO x ：，
wherein F is 1 Representing the vehicle speed; f (F) 2 Representing acceleration; f (F) 3 Representing torque; f (F) 4 Indicating the rotation speed; f (F) 5 Indicating grade.
For target pollutants carbon monoxide (CO), carbon dioxide (CO) 2 ) Oxides of Nitrogen (NO) x ) Fitting parameters FFP based on characteristics of corresponding target contaminants p And two variables of engine load, and dividing vehicle Bin intervals to obtain 19 Bin intervals in totalAs shown in fig. 9, fig. 9 shows a Bin section schematic constructed for a light truck using a Bin section division method two.
In fig. 9, for carbon monoxide (CO), carbon dioxide (CO 2 ) Oxides of Nitrogen (NO) x ) Constant target pollutant, characteristic fitting parameter FFP p The range of values of (2) is different because the target pollutant is different, and the characteristic fitting parameter FFP is calculated p When the condition mutual information value of the vehicle characteristic parameter and the target pollutant P is different, the characteristic fitting parameter FFP is caused p The range of values of (2) is different.
Therefore, when the emission factors of different target pollutants P are calculated later, the parameters FFP are fitted according to the characteristics of the target pollutants P according to the different target pollutants P p And dividing the actual road emission data.
S3, calculating a heavy truck reference emission factor serving as the target vehicle according to the actual emission data of the heavy truck and the multi-scale Bin interval division result, wherein the actual emission data and the multi-scale Bin interval division result are obtained in the previous steps.
，
In the method, in the process of the invention,the standard emission factor of the target pollutant p in g/km under the jth Bin section of the stage operation stage is represented, wherein the stage represents three vehicle operation stages of cold start, low load, medium and high load and the like;representing the discharge amount of the target pollutant p in g in the jth Bin interval of the stage operation stage, wherein the discharge amount is obtained by summing the discharge rates of the pollutants p in the Bin intervals;the vehicle mileage in km at the jth Bin section in the stage of operation is represented.
The calculation result of the reference emission factor of the heavy truck as the target vehicle is shown in fig. 3, 4 and 5, and CO are covered 2 、NO x And the reference emission factors of three pollutants.
S4, based on a standard emission factor of the heavy-duty truck serving as a target vehicle and a Bin interval driving mileage ratio, respectively calculating a cold start emission factor, a low load emission factor, a medium and high load emission factor and a comprehensive emission factor by the following two formulas, and further constructing a multi-scale emission factor set in the target area, wherein the multi-scale emission factor set is shown in a table 1.
，
，
In the method, in the process of the invention,the emission factor of the target pollutant p in the stage operation stage is expressed in g/km;the comprehensive emission factor of the target pollutant p is expressed in g/km;representing the driving mileage duty ratio of the jth Bin interval;and the correction coefficient representing the reference emission factor is equal to the product of the average running speed correction factor, the weather information correction factor and the oil standard correction factor in the target area.
The Bin range is shown in fig. 6. Further, for CO, the correction coefficient is 0.68; for CO 2 The correction coefficient is 0.9; for NO x The correction coefficient was 0.71. Based on this, a set of constructed multi-scale emission factors capable of characterizing the emission characteristics of the target region is shown in Table 1Shown.
TABLE 1
| Emission factor | CO/(g/km) | CO 2 /(g/km) | NO x /(g/km) |
| Cold start emission factor | 0.008 | 2.96 | 0.012 |
| Low load emission factor | 0.095 | 138.55 | 0.026 |
| Medium and high load emission factor | 0.099 | 149.13 | 0.002 |
| Comprehensive emission factor | 0.202 | 290.64 | 0.040 |
S5, acquiring the driving mileage of the bicycle in the target area, and calculating the product of the multi-scale emission factor set and the driving mileage of the bicycle by the following two formulas to obtain the pollutant emission intensity and the comprehensive emission intensity of the bicycle in the target area in the cold start, low load, medium and high load operation stages.
，
，
In the method, in the process of the invention,the emission intensity of the target pollutant p in the stage operation stage is expressed as g;the comprehensive emission intensity of the target pollutant p is expressed in g;indicating the mileage of the bicycle. The list of single vehicle emissions for a heavy utility vehicle in the target area is shown in fig. 7 and includes carbon monoxide (CO), carbon dioxide (CO) 2 ) Oxides of Nitrogen (NO) x ) Respective cold start emissions, low load emissions, medium and high load emissions, and integrated emissions.
S6, acquiring the single vehicle driving mileage of all vehicles in the target area, repeating the step S5, and calculating to obtain a single vehicle emission list of all vehicles. And adding the single vehicle emission intensities of all vehicles, so as to calculate and obtain the pollutant emission intensity and the comprehensive emission intensity of all vehicles in the target area in the cold start, low load, medium and high load operation stages, and further construct a multi-scale vehicle pollutant emission list in the target area. As shown in fig. 8, a multi-scale vehicle exhaust emissions list is provided for all heavy-duty trucks in the target area, including carbon monoxide (CO), carbon dioxide (CO 2 ) Oxides of Nitrogen (NO) x ) Respective cold start emissions, low load emissions, and medium and high load emissionsDischarge amount and comprehensive discharge amount.
The above description of construction of the multi-scale vehicle exhaust emission list by using the heavy truck as the target vehicle is the same as the construction of the multi-scale vehicle exhaust emission list by other vehicles except for the Bin interval division method, and will not be described again.
While the fundamental and principal features of the invention and advantages of the invention have been shown and described, it will be apparent to those skilled in the art that the invention is not limited to the details of the foregoing exemplary embodiments, but may be embodied in other specific forms without departing from the spirit or essential characteristics thereof;
the present embodiments are, therefore, to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein.
Furthermore, it should be understood that although the present disclosure describes embodiments, not every embodiment is provided with a separate embodiment, and that this description is provided for clarity only, and that the disclosure is not limited to the embodiments described in detail below, and that the embodiments described in the examples may be combined as appropriate to form other embodiments that will be apparent to those skilled in the art.
Claims (5)
1. A vehicle pollutant emission inventory accounting method, comprising the steps of:
s1, dividing the vehicle operation stage into cold start, low load and medium and high load;
the cold start operation phase is defined as a vehicle operation phase in a preset time period after the engine is started when the temperature of the engine coolant does not exceed a preset temperature threshold value, and is set as Bin1; the low-load operation stage is a vehicle operation stage in which the engine load is lower than a full-load preset percentage value in a non-cold start operation stage; the medium-high load operation phase is a vehicle operation phase except for cold start and low load;
s2, when a pollutant emission list of a target vehicle is calculated, selecting a vehicle Bin interval from two mutually independent Bin interval dividing methods according to the model specification of the target vehicle in the low-load and medium-high-load operation stages, and establishing a multi-scale Bin interval of the target vehicle;
for a first target vehicle, constructing a multi-scale Bin interval by adopting a Bin interval dividing method I; the Bin interval dividing method comprises the following steps:
for a low-load operation stage, carrying out Bin interval division based on two variables of engine load and vehicle speed;
for the medium-high load operation stage, dividing a Bin interval is carried out based on a vehicle speed per acceleration VSPA and two variables of the vehicle speed, and the vehicle speed per acceleration VSPA per second reflects the driving excitation degree of the vehicle;
for a second target vehicle, constructing a multi-scale Bin interval by adopting a Bin interval dividing method II; the second method for dividing the Bin section adopts a mutual information method, fuses vehicle characteristic parameters including vehicle speed, acceleration, gradient, rotating speed and torque, constructs a characteristic fitting parameter FFP, and simultaneously, develops the division of the vehicle Bin section based on the characteristic fitting parameter FFP and two variables of engine load, and comprises the following steps:
according to each parameter change range [ min, max ] set in advance, carrying out normalization processing on a vehicle characteristic parameter F and a target pollutant p to be calculated, and mapping data to a [0,1] interval, wherein the vehicle characteristic parameter F comprises vehicle speed, acceleration, gradient, rotating speed and torque;
calculating the condition mutual information h between the vehicle characteristic parameter F and the target pollutant p one by one F,p As a coefficient of the vehicle characteristic parameter F, wherein the condition mutual information h F,p Representing the correlation coefficient of the vehicle characteristic parameter F and the target pollutant p;
for all the vehicle characteristic parameters F, calculating mutual information I between every two parameters a,b Wherein, the method comprises the steps of, wherein,calculate the sum of all mutual information +.>As a constant term coefficient, where mutual information I a,b Representing a redundant relationship between the vehicle characteristic parameters F;
obtaining characteristic fitting parameters FFP for target pollutant p p The formula of (2) is shown below;
；
s3, acquiring actual emission data of a target vehicle, and calculating a vehicle reference emission factor of the target vehicle according to the constructed multi-scale Bin interval;
，
in the method, in the process of the invention,the standard emission factor of the target pollutant p in g/km under the jth Bin section of the stage operation stage is represented; the stage represents three vehicle operation stages of cold start, low load, medium load and high load; />The emission quantity of the target pollutant p in the jth Bin section of the stage operation stage is expressed in g, and the emission quantity is obtained by summing the emission rates of the target pollutant p in the corresponding Bin section every second; />The vehicle driving mileage under the jth Bin section in the stage operation stage is expressed, and the unit is km;
s4, calculating a cold start emission factor, a low load emission factor, a medium and high load emission factor and a comprehensive emission factor based on a vehicle reference emission factor and a Bin interval driving mileage ratio, and constructing a multi-scale emission factor set in a target area;
，
，
in the method, in the process of the invention,the emission factor of the target pollutant p in the stage operation stage is expressed in g/km; />The comprehensive emission factor of the target pollutant p is expressed in g/km; />Representing the driving mileage duty ratio of the jth Bin interval; />The correction coefficient representing the reference emission factor is equal to the product of the average running speed correction factor, the weather information correction factor and the oil standard correction factor in the target area;
s5, acquiring the driving mileage of the bicycle in the target area, calculating the product of the multi-scale emission factor set and the driving mileage of the bicycle in the target area, and acquiring the pollutant emission intensity and the comprehensive emission intensity of the bicycle in the cold start, low load and medium and high load operation stages in the target area;
，
，
in the method, in the process of the invention,the emission intensity of the target pollutant p in the stage operation stage is expressed as g; />The comprehensive emission intensity of the target pollutant p is expressed in g; />Representing the driving mileage of a bicycle;
s6, obtaining the driving mileage of all vehicles in the target area, repeating the step S5, calculating and obtaining the pollutant emission intensity and the comprehensive emission intensity of all vehicles in the target area in the cold start, low load, medium and high load operation stages, and constructing a multi-scale vehicle pollutant emission list in the target area.
2. The vehicle pollutant emission inventory accounting method of claim 1, wherein said preset temperature threshold is 70 ℃ and said preset time period is 300s.
3. The vehicle pollutant emission inventory accounting method of claim 1, wherein said predetermined percentage value is 40%.
4. The vehicle pollutant emission inventory accounting method of claim 1, wherein said vehicle per acceleration VSPA is obtained by calculation of the following formula;
wherein v represents the vehicle speed; a represents acceleration.
5. The vehicle pollutant emission inventory accounting method of claim 1, wherein said first target vehicle comprises a large passenger vehicle, a bus, a medium cargo vehicle, and a heavy cargo vehicle; the second target vehicle comprises a mini-type passenger car, a taxi, a medium-sized passenger car, a mini-type cargo car, a light-duty cargo car and a motorcycle.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311386470.1A CN117112969B (en) | 2023-10-25 | 2023-10-25 | Vehicle pollutant emission list accounting method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311386470.1A CN117112969B (en) | 2023-10-25 | 2023-10-25 | Vehicle pollutant emission list accounting method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117112969A CN117112969A (en) | 2023-11-24 |
| CN117112969B true CN117112969B (en) | 2024-01-02 |
Family
ID=88813239
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311386470.1A Active CN117112969B (en) | 2023-10-25 | 2023-10-25 | Vehicle pollutant emission list accounting method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117112969B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19753006A1 (en) * | 1997-11-30 | 1999-06-17 | Wwu Wissenschaftliche Werkstat | Measuring exhaust gases of any vehicle or craft to determine efficiency of anti-pollution components e.g. catalysts |
| CN105806627A (en) * | 2016-03-15 | 2016-07-27 | 南京汽车集团有限公司 | Vehicle-mounted emission testing device and testing method thereof |
| CN110209990A (en) * | 2019-05-22 | 2019-09-06 | 中山大学 | A kind of single vehicle trajectory of discharge calculation method based on vehicle identification detection data |
| KR102238001B1 (en) * | 2020-10-30 | 2021-04-08 | 재단법인 서울특별시 서울기술연구원 | Device and method for calculating air pollutant emissions from road mobile source |
| CN116070839A (en) * | 2022-12-23 | 2023-05-05 | 成都市机动车排气污染防治技术保障中心 | Multi-source data road micro-environment air quality contribution tracing method |
| CN116307837A (en) * | 2023-02-10 | 2023-06-23 | 中汽研汽车检验中心(昆明)有限公司 | Regional motor vehicle emission evaluation method and system based on multiple monitoring technologies |
-
2023
- 2023-10-25 CN CN202311386470.1A patent/CN117112969B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE19753006A1 (en) * | 1997-11-30 | 1999-06-17 | Wwu Wissenschaftliche Werkstat | Measuring exhaust gases of any vehicle or craft to determine efficiency of anti-pollution components e.g. catalysts |
| CN105806627A (en) * | 2016-03-15 | 2016-07-27 | 南京汽车集团有限公司 | Vehicle-mounted emission testing device and testing method thereof |
| CN110209990A (en) * | 2019-05-22 | 2019-09-06 | 中山大学 | A kind of single vehicle trajectory of discharge calculation method based on vehicle identification detection data |
| KR102238001B1 (en) * | 2020-10-30 | 2021-04-08 | 재단법인 서울특별시 서울기술연구원 | Device and method for calculating air pollutant emissions from road mobile source |
| CN116070839A (en) * | 2022-12-23 | 2023-05-05 | 成都市机动车排气污染防治技术保障中心 | Multi-source data road micro-environment air quality contribution tracing method |
| CN116307837A (en) * | 2023-02-10 | 2023-06-23 | 中汽研汽车检验中心(昆明)有限公司 | Regional motor vehicle emission evaluation method and system based on multiple monitoring technologies |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117112969A (en) | 2023-11-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Delgado et al. | Fuel efficiency technology in european heavy-duty vehicles: Baseline and potential for the 2020–2030 time frame | |
| Ang-Olson et al. | Energy efficiency strategies for freight trucking: potential impact on fuel use and greenhouse gas emissions | |
| Sun et al. | Energy consumption simulation and economic benefit analysis for urban electric commercial-vehicles | |
| Tu et al. | Real-world emissions and fuel consumption of gasoline and hybrid light duty vehicles under local and regulatory drive cycles | |
| Weller et al. | Real world fuel consumption and emissions from LDVs and HDVs | |
| Guillou et al. | Fuel consumption testing to verify the effect of tire rolling resistance on fuel economy | |
| CN112486962A (en) | Extraction and combination short segment calculation heavy-duty diesel vehicle NOxMethod of discharging | |
| Delucchi et al. | Electric and gasoline vehicle lifecycle cost and energy-use model | |
| CN108052791A (en) | The evaluation method that following vehicle structure adjustment influences traffic emission source inventory | |
| CN114088878A (en) | Vehicle carbon emission monitoring system based on 5G technology and monitoring method thereof | |
| CN117112969B (en) | Vehicle pollutant emission list accounting method | |
| CN115876484A (en) | System and method for actual road condition test and working condition simulation test of heavy whole vehicle | |
| WO2018214788A1 (en) | Method and system for maintaining qualified tail gas emission of automobile by using big data | |
| WO2014062127A1 (en) | Determination of consumption of energy for a vehicle | |
| CN117114654A (en) | Configured vehicle maintenance implementation method, storage medium and device | |
| CN112945567A (en) | Method and system for predicting vehicle-mounted emission of low-temperature diesel engine | |
| CN114038085B (en) | Vehicle maintenance management system and method | |
| CN117093807B (en) | Vehicle exhaust emission list calculation method based on multi-dimension Bin interval | |
| CN113283506A (en) | Vehicle-mounted TBOX load state identification method based on road working conditions | |
| CN105134397B (en) | Diagnostic method, system and the vehicle with the system of lambda sensor | |
| Chandra et al. | A Study on Vehicle Tire Inflation and Fuel Consumption | |
| Li et al. | Analysis of Nitrogen Oxide Emission Characteristics of Heavy-Duty Refrigerated Trucks Based on Remote Monitoring Data | |
| Farkas | Fuel efficiency in truck industry | |
| Ragon et al. | emissions from trucks in the EU: An analysis of the heavy-duty CO | |
| CN117072295A (en) | Early warning method and device for frequent regeneration of DPF |
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 |