[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN104815545A - Real-time desulfurizing agent and by product flow calculation method for wet desulphurization system - Google Patents

Real-time desulfurizing agent and by product flow calculation method for wet desulphurization system Download PDF

Info

Publication number
CN104815545A
CN104815545A CN201510198049.7A CN201510198049A CN104815545A CN 104815545 A CN104815545 A CN 104815545A CN 201510198049 A CN201510198049 A CN 201510198049A CN 104815545 A CN104815545 A CN 104815545A
Authority
CN
China
Prior art keywords
real
fgdout
fgdin
time
yan
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.)
Granted
Application number
CN201510198049.7A
Other languages
Chinese (zh)
Other versions
CN104815545B (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.)
Xian TPRI Power Station Information Technology Co Ltd
Original Assignee
Xian TPRI Power Station Information 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 Xian TPRI Power Station Information Technology Co Ltd filed Critical Xian TPRI Power Station Information Technology Co Ltd
Priority to CN201510198049.7A priority Critical patent/CN104815545B/en
Publication of CN104815545A publication Critical patent/CN104815545A/en
Application granted granted Critical
Publication of CN104815545B publication Critical patent/CN104815545B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Treating Waste Gases (AREA)

Abstract

The invention discloses a real-time desulfurizing agent and by product flow calculation method for a wet desulphurization system. According to the method, flue gas desulphurization (FGD) tower inlet work condition SO2 concentration SO [2] [fgdIn], FGD tower inlet glue gas V [yan, fgdIn] produced in each kilogram coal work condition, FGD tower outlet work condition SO2 concentration SO [2] [fgdOut], FGD tower outlet glue gas V [yan, fgdOut] produced in each kilogram coal work condition, fired coal real-time caloric value Qar and generator set active power Ne and power generation coal consumption Bf are subjected to real-time calculation to obtain the desulfurizing agent and by product flow of the wet desulphurization system. The basis for desulfurizing agent and by product calculation is provided by obtaining the desulfurizing agent and by product flow through calculation.

Description

A kind of method of real-time calculating wet desulfurization system desulfurizing agent and byproduct flow
Technical field:
The present invention relates to the fields such as heat energy, chemistry, information technology, be specifically related to a kind of method of real-time calculating wet desulfurization system desulfurizing agent and byproduct flow.
Background technology:
At present, thermal power plant's lime stone (lime)-gypsum wet desulphurization system sorbent consumption amount adopts the method for belt conveyer scale or coal powder silo with level meter, and its accuracy and real-time are difficult to ensure, and gypsum limits due to production technology, measures difficulty in real time larger.
Summary of the invention:
The object of the present invention is to provide a kind of method of real-time calculating wet desulfurization system desulfurizing agent and byproduct flow, utilize FGD absorption tower entrance condition condition SO 2concentration SO2 fgdIn, the flue gas volume V that produces of FGD absorption tower entrance every kilogram of coal working condition yan, fgdIn, FGD absorption tower outlet condition condition SO 2concentration SO2 fgdOut, FGD absorption tower exports the flue gas volume V that every kilogram of coal working condition produces yan, fgdOut, the real-time calorific value Q of as-fired coal ar, generating set active power N eand unit generation coal consumption B f, calculate sorbent consumption amount and byproduct generation in real time, for sorbent consumption amount and byproduct generation statistics provide foundation.
For achieving the above object, the present invention adopts following technical scheme to be achieved:
A method for real-time calculating wet desulfurization system desulfurizing agent and byproduct flow, gathers FGD absorption tower entrance condition condition SO 2concentration SO 2fgdIn, the flue gas volume V that produces of FGD absorption tower entrance every kilogram of coal working condition yan, fgdIn, FGD absorption tower outlet condition condition SO 2concentration SO 2fgdOut, FGD absorption tower exports the flue gas volume V that every kilogram of coal working condition produces yan, fgdOut, the real-time calorific value Q of as-fired coal ar, generating set active power N eand unit generation coal consumption B f, utilizing formula (1) to obtain take lime stone as the real-time traffic of desulfurizing agent:
CaCO 3 ( m ) = 45.794 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) × Ca / S η CaCO 3 - - - ( 1 )
In formula: N e---generating set active power, MW;
B f---the real-time gross coal consumption rate of unit, g/kWh;
V yan, fgdIn---the flue gas volume that entrance every kilogram of coal working condition in FGD absorption tower produces, m 3/ kg;
V yan, fgdOut---FGD absorption tower exports the flue gas volume that every kilogram of coal working condition produces, m 3/ kg;
SO 2fgdIn---FGD absorption tower entrance condition condition SO 2concentration, mg/m 3;
SO 2fgdOut---FGD absorption tower outlet condition condition SO 2concentration, mg/m 3;
Q ar---the real-time calorific value of as-fired coal, kJ/kg;
Ca/S---calcium to sulphur mole ratio ,/;
---lime stone purity ,/;
CaCO 3(m)---the real-time calculated mass flow of lime stone, t/h;
Utilizing formula (2) to obtain take lime as the real-time traffic of desulfurizing agent:
CaO ( m ) = 25.665 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) × Ca / S η CaO - - - ( 2 )
In formula: η caO---lime purity ,/;
CaO (m)---the real-time calculated mass flow of lime, t/h;
Formula (3) is utilized to obtain desulfuration byproduct real-time traffic:
CaSO 4 · 2 H 2 O ( m ) = 78.765 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) / η sg - - - ( 3 )
In formula: η sg---gypsum purity ,/;
CaSO 42H 2o (m)---the plaster of paris calculates generation, t/h.
The present invention further improves and is,
When gateway, FGD absorption tower have the real-time measuring point of flue gas total flow or by the total flue gas flow calculated can meet precision be greater than 95% time, then:
Utilizing formula (4) to obtain take lime stone as the real-time traffic of desulfurizing agent:
CaCO 3 ( m ) = 1.562 × 10 - 3 × ( V Y , fgdIn × SO 2 fgdIn - V Y , fgdOut × SO 2 fgdOut ) η NH 3 × × Ca / S η CaCO 3 - - - ( 4 )
Utilizing formula (5) to obtain take lime as the real-time traffic of desulfurizing agent:
CaO ( m ) = 8.75 × 10 - 4 × ( V Y , fgdIn × SO 2 fgdIn - V Y , fgdOut × SO 2 fgdOut ) η NH 3 × × Ca / S η CaO - - - ( 5 )
Formula (6) is utilized to obtain desulfuration byproduct real-time traffic:
CaSO 4·2H 2O(m)=2.688×10 -3×(V Y,fgdIn×SO 2fgdIn-V Y,fgdOut×SO 2fgdOut)/η sg(6)
In formula: V y, fgdIn---porch, FGD absorption tower flue gas volume flow, km 3/ h;
V y, fgdOut---exit, FGD absorption tower flue gas volume flow, km 3/ h.
Relative to prior art, the present invention has the following advantages: a kind of method of real-time calculating wet desulfurization system desulfurizing agent and byproduct flow, utilizes FGD absorption tower entrance condition condition SO 2concentration SO 2fgdIn, the flue gas volume V that produces of FGD absorption tower entrance every kilogram of coal working condition yan, fgdIn, FGD absorption tower outlet condition condition SO 2concentration SO 2fgdOut, FGD absorption tower exports the flue gas volume V that every kilogram of coal working condition produces yan, fgdOut, the real-time calorific value Q of as-fired coal ar, generating set active power N eand unit generation coal consumption B f, calculate the real-time consumption of desulfurizing agent and byproduct generation in real time; Result of calculation of the present invention may be used for desulfurizing agent, byproduct statistics provides a kind of method; Realization of the present invention obtains the real-time generation of desulfuration byproduct first.
Detailed description of the invention:
The method of a kind of real-time calculating wet desulfurization system desulfurizing agent of the present invention and byproduct flow, gathers FGD absorption tower entrance condition condition SO 2concentration SO 2fgdIn, the flue gas volume V that produces of FGD absorption tower entrance every kilogram of coal working condition yan, fgdIn, FGD absorption tower outlet condition condition SO 2concentration SO 2fgdOut, FGD absorption tower exports the flue gas volume V that every kilogram of coal working condition produces yan, fgdOut, the real-time calorific value Q of as-fired coal ar, generating set active power N eand unit generation coal consumption B f, calculate desulfurizing agent flow in real time.
1) lime stone flow rate calculation
Utilizing formula (1) to obtain take lime stone as the real-time traffic of desulfurizing agent:
CaCO 3 ( m ) = 45.794 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) × Ca / S η CaCO 3 - - - ( 1 )
In formula: N e---generating set active power, MW;
B f---the real-time gross coal consumption rate of unit, g/kWh;
V yan, fgdIn---the flue gas volume that entrance every kilogram of coal working condition in FGD absorption tower produces, m 3/ kg;
V yan, fgdOut---FGD absorption tower exports the flue gas volume that every kilogram of coal working condition produces, m 3/ kg;
SO 2fgdIn---FGD absorption tower entrance condition condition SO 2concentration, mg/m 3;
SO 2fgdOut---FGD absorption tower outlet condition condition SO 2concentration, mg/m 3;
Q ar---the real-time calorific value of as-fired coal, kJ/kg;
Ca/S---calcium to sulphur mole ratio ,/;
η caCO3---lime stone purity ,/;
CaCo 3(m)---the real-time calculated mass flow of lime stone, t/h;
2) lime flow rate calculation
Utilizing formula (2) to obtain take lime as the real-time traffic of desulfurizing agent:
CaO ( m ) = 25.665 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) × Ca / S η CaO - - - ( 2 )
In formula: η caO---lime purity ,/;
CaO (m)---the real-time calculated mass flow of lime, t/h;
3) gypsum generation calculates
Formula (3) is utilized to obtain desulfuration byproduct gypsum real-time traffic:
CaSO 4 · 2 H 2 O ( m ) = 78.765 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) / η sg - - - ( 3 )
In formula: η sg---gypsum purity ,/;
CaSO 42H 2o (m)---the plaster of paris calculates generation, t/h;
4) related description
In view of SO 2real-time concentration is measured under being working condition, therefore with V in above formula (1), formula (2), formula (3) yan, fgdIn, V yan, fgdOutthe dry flue gas volume flow produced by every kilogram of coal mark condition condition, inleakage and water vapour content three sum, then be converted to the flue gas volume under working condition, its specific algorithm can with reference to fired coal combustion equation and exhaust gas volumn calculating etc.
When gateway, FGD absorption tower have the real-time measuring point of flue gas total flow or by the total flue gas flow calculated can meet precision be greater than 95% time, then:
Lime stone mass flow calculation can adopt:
CaCO 3 ( m ) = 1.562 × 10 - 3 × ( V Y , fgdIn × SO 2 fgdIn - V Y , fgdOut × SO 2 fgdOut ) η NH 3 × × Ca / S η CaCO 3 - - - ( 4 )
Lime quality flow rate calculation can adopt:
CaO ( m ) = 8.75 × 10 - 4 × ( V Y , fgdIn × SO 2 fgdIn - V Y , fgdOut × SO 2 fgdOut ) η NH 3 × × Ca / S η CaO - - - ( 5 )
Plaster of paris mass flow calculation can adopt:
CaSO 4·2H 2O(m)=2.688×10 -3×(V Y,fgdIn×SO 2fgdIn-V Y,fgdOut×SO 2fgdOut)/η sg(6)
In formula: V y, fgdIn---porch, FGD absorption tower flue gas volume flow, km 3/ h;
V y, fgdOut---exit, FGD absorption tower flue gas volume flow, km 3/ h.
5) stability calculated and accuracy
In order to ensure the accurately fixed of result of calculation and stability, FGD absorption tower entrance condition condition SO 2concentration, FGD absorption tower outlet condition condition SO 2the real time datas such as the real-time gross coal consumption rate of concentration, unit, the real-time calorific value of as-fired coal adopt the mean value in certain hour to calculate, and Real-time Collection needs verification, then substitutes into formulae discovery.Real-time result of calculation may be used for being analyzed with measured data and providing a kind of method for grasping FGD system byproduct real-time traffic.

Claims (2)

1. calculate a method for wet desulfurization system desulfurizing agent and byproduct flow in real time, it is characterized in that, gather FGD absorption tower entrance condition condition SO 2concentration SO 2fgdIn, the flue gas volume V that produces of FGD absorption tower entrance every kilogram of coal working condition yan, fgdIn, FGD absorption tower outlet condition condition SO 2concentration SO 2fgdOut, FGD absorption tower exports the flue gas volume V that every kilogram of coal working condition produces yan, fgdOut, the real-time calorific value Q of as-fired coal ar, generating set active power N eand unit generation coal consumption B f, utilizing formula (1) to obtain take lime stone as the real-time traffic of desulfurizing agent:
CaCO 3 ( m ) = 45.797 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) × Ca / S η CaC O 3 - - - ( 1 )
In formula: N e---generating set active power, MW;
B f---the real-time gross coal consumption rate of unit, g/kWh;
V yan, fgdIn---the flue gas volume that entrance every kilogram of coal working condition in FGD absorption tower produces, m 3/ kg;
V yan, fgdOut---FGD absorption tower exports the flue gas volume that every kilogram of coal working condition produces, m 3/ kg;
SO 2fgdIn---FGD absorption tower entrance condition condition SO 2concentration, mg/m 3;
SO 2fgdOut---FGD absorption tower outlet condition condition SO 2concentration, mg/m 3;
Q ar---the real-time calorific value of as-fired coal, kJ/kg;
Ca/S---calcium to sulphur mole ratio ,/;
---lime stone purity ,/;
CaCO 3(m)---the real-time calculated mass flow of lime stone, t/h;
Utilizing formula (2) to obtain take lime as the real-time traffic of desulfurizing agent:
CaO ( m ) = 25.665 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) × Ca / S η CaO - - - ( 2 )
In formula: η caO---lime purity ,/;
CaO (m)---the real-time calculated mass flow of lime, t/h;
Formula (3) is utilized to obtain desulfuration byproduct real-time traffic:
CaSO 4 · 2 H 2 O ( m ) = 78.765 × 10 - 6 × N e × B f Q ar × ( V yan , fgdIn × SO 2 fgdIn - V yan , fgdOut × SO 2 fgdOut ) / η sg - - - ( 3 )
In formula: η sg---gypsum purity ,/;
CaSO 42H 2o (m)---the plaster of paris calculates generation, t/h.
2. the method for a kind of real-time calculating wet desulfurization system desulfurizing agent according to claim 1 and byproduct flow, is characterized in that,
When gateway, FGD absorption tower have the real-time measuring point of flue gas total flow or by the total flue gas flow calculated can meet precision be greater than 95% time, then:
Utilizing formula (4) to obtain take lime stone as the real-time traffic of desulfurizing agent:
CaCO 3 ( m ) = 1.562 × 10 - 3 × ( V Y , fgdIn × SO 2 fgdIn - V Y , fgdOut × SO 2 fgdOut ) η NH 3 × × Ca / s η CaCO 3 - - - ( 4 )
Utilizing formula (5) to obtain take lime as the real-time traffic of desulfurizing agent:
CaO ( m ) = 8.75 × 10 - 4 × ( V Y , fgdIn × SO 2 fgdIn - V Y , fgdOut × SO 2 fgdOut ) η NH 3 × × Ca / s η CaO - - - ( 5 )
Formula (6) is utilized to obtain desulfuration byproduct real-time traffic:
CaSO 4·2H 2O(m)=2.688×10 -3×(V Y,fgdIn×SO 2fgdIn-V Y,fgdOut×SO 2fgdOut)/η sg(6)
In formula: V y, fgdIn---porch, FGD absorption tower flue gas volume flow, km 3/ h;
V y, fgdOut---exit, FGD absorption tower flue gas volume flow, km 3/ h.
CN201510198049.7A 2015-04-23 2015-04-23 A kind of real-time method for calculating wet desulfurization system desulfurizing agent and byproduct flow Expired - Fee Related CN104815545B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510198049.7A CN104815545B (en) 2015-04-23 2015-04-23 A kind of real-time method for calculating wet desulfurization system desulfurizing agent and byproduct flow

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510198049.7A CN104815545B (en) 2015-04-23 2015-04-23 A kind of real-time method for calculating wet desulfurization system desulfurizing agent and byproduct flow

Publications (2)

Publication Number Publication Date
CN104815545A true CN104815545A (en) 2015-08-05
CN104815545B CN104815545B (en) 2017-03-15

Family

ID=53726168

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510198049.7A Expired - Fee Related CN104815545B (en) 2015-04-23 2015-04-23 A kind of real-time method for calculating wet desulfurization system desulfurizing agent and byproduct flow

Country Status (1)

Country Link
CN (1) CN104815545B (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105536490A (en) * 2015-12-23 2016-05-04 北京首钢自动化信息技术有限公司 Automatic desulfurization system and method for sintering flue gas
CN105536488A (en) * 2015-12-11 2016-05-04 新疆生产建设兵团农八师天山铝业有限公司 Control method for SO2 during semi-dry desulfurization
CN107308803A (en) * 2017-08-16 2017-11-03 华能白山煤矸石发电有限公司 A kind of circulating fluidized bed desulfurization tower and control method
CN108305105A (en) * 2018-03-13 2018-07-20 广州华润热电有限公司 A kind of consumptive material predictor method and processing equipment

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005160171A (en) * 2003-11-25 2005-06-16 Mitsubishi Electric Corp Industrial energy management system
US20050186605A1 (en) * 1996-06-19 2005-08-25 Short Jay M. Phosphatase -encoding nucleic acids and methods of making and using them
CN101082415A (en) * 2006-11-30 2007-12-05 长沙市朝槿科技有限公司 Method for real time centralized controlling electric generating set based on SO2 discharging
CN101122590A (en) * 2007-08-17 2008-02-13 东南大学 Limestone activity test method and test device for wet process of FGD
CN103390098A (en) * 2013-06-28 2013-11-13 武汉钢铁(集团)公司 Desulfurization computing method and system applying same

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050186605A1 (en) * 1996-06-19 2005-08-25 Short Jay M. Phosphatase -encoding nucleic acids and methods of making and using them
JP2005160171A (en) * 2003-11-25 2005-06-16 Mitsubishi Electric Corp Industrial energy management system
CN101082415A (en) * 2006-11-30 2007-12-05 长沙市朝槿科技有限公司 Method for real time centralized controlling electric generating set based on SO2 discharging
CN101122590A (en) * 2007-08-17 2008-02-13 东南大学 Limestone activity test method and test device for wet process of FGD
CN103390098A (en) * 2013-06-28 2013-11-13 武汉钢铁(集团)公司 Desulfurization computing method and system applying same

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105536488A (en) * 2015-12-11 2016-05-04 新疆生产建设兵团农八师天山铝业有限公司 Control method for SO2 during semi-dry desulfurization
CN105536488B (en) * 2015-12-11 2018-04-27 新疆生产建设兵团农八师天山铝业有限公司 A kind of control method using semi-dry desulphurization SO2
CN105536490A (en) * 2015-12-23 2016-05-04 北京首钢自动化信息技术有限公司 Automatic desulfurization system and method for sintering flue gas
CN105536490B (en) * 2015-12-23 2018-02-27 北京首钢自动化信息技术有限公司 The automatic desulphurization system of one kind sintering flue gas and its method
CN107308803A (en) * 2017-08-16 2017-11-03 华能白山煤矸石发电有限公司 A kind of circulating fluidized bed desulfurization tower and control method
CN108305105A (en) * 2018-03-13 2018-07-20 广州华润热电有限公司 A kind of consumptive material predictor method and processing equipment

Also Published As

Publication number Publication date
CN104815545B (en) 2017-03-15

Similar Documents

Publication Publication Date Title
CN103697946B (en) A kind of computing method of coal fired boiler of power plant flue gas flow and the control method of pollutant discharge amount
CN104815545A (en) Real-time desulfurizing agent and by product flow calculation method for wet desulphurization system
CN104215546B (en) A kind of power station boiler air pre-heater stifled ash monitoring system and method for work thereof
CN101865867B (en) Method for calculating coal elements and industrial components in real time
CN103674103A (en) Performance testing method for SCR (Selective Catalytic Reduction) flue gas denitration system
CN105920997B (en) A kind of coal-burning boiler denitrating system that burnout degree is coupled with SNCR and method
CN104699937A (en) Boiler efficiency self-correction computing method based on flue gas testing
CN103699780B (en) Ature of coal parameter is in the chaos optimization method of line computation
CN103955751A (en) Method for predicting multi-element nonlinear limestone-gypsum wet desulphurization efficiency
CN105242000A (en) Method for accurate measurement of carbon emission from coal-fired power plant
CN105457465A (en) Design method of SNCR-SCR combined flue gas denitrification system process
CN105536490B (en) The automatic desulphurization system of one kind sintering flue gas and its method
CN105138718B (en) A kind of projectional technique and its subsidiary engine load adjustment method of combination UDF Novel fires desulfurizing tower of electric power plant desulfuration efficiency
CN105160479A (en) Method for measuring boiler dry flue gas volume based on coal caloric values and exhaust fume oxygen amount
CN101229483B (en) Method of preventing serosity toxicosis through determining and absorbing free F- or Al3+ content in serosity
CN103674600B (en) A kind of monitoring method of coal-burning boiler operational efficiency
CN104809341B (en) A kind of real-time method for calculating SCR denitration system reduction agent flux
CN102721798B (en) Method for calculating coal quality sulfur component in real time
CN105964123B (en) A kind of ammonia process of desulfurization spray Deethanizer design and the computation model and optimization regulating method of operation
Huang et al. Aerosol formation by heterogeneous reactions in ammonia-based WFGD systems
CN105092827B (en) Method for evaluating grade of flue gas desulfurization gypsum
CN114626734B (en) Carbon emission evaluation method for garbage incineration power plant
CN101780364A (en) On-line calculation method for desulfurization efficiency of circulating fluidized bed boiler
CN105043799A (en) Emission performance test method of coal-fired boiler unit
Yoon et al. Characteristics of CO2 capture system using KIERSOL in the LNG flue gas

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
EXSB Decision made by sipo to initiate substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170315

Termination date: 20180423