CN101608996A - The assay method of inland non-industrial atmospheric corrosiveness - Google Patents

Abstract

The present invention relates to a kind of assay method of inland non-industrial atmospheric corrosiveness, it is characterized in that following steps are arranged: 1) determine the monitoring point, extract the nearest 10 years diary temperature of this area's meteorological station and the time Value Data of relative humidity; 2) temperature, humidity data arrangement; 3) calculate the wetting time value of each year according to the wetting time formula; 4), calculate the average and the standard deviation of wetting time according to the above-mentioned gained wetting time value of each year; 5) determine wetting time grade according to average, determine the variation tendency of wetting time according to standard deviation; 6) determine the atmospheric corrosiveness grade of monitoring point according to wetting time grade.The Value Data when temperature that adopts the inventive method can make full use of observatory to record, humidity, the non-industrial atmospheric result who is surveyed is accurate, and method is simple.

Description

The assay method of inland non-industrial atmospheric corrosiveness

Technical field

The present invention be more particularly directed to a kind of method of measuring atmospheric environment to effect of metal corrosion, particularly a kind of assay method of inland non-industrial atmospheric corrosiveness.

Background technology

When relative humidity reached a certain value in the air, the metal surface can form certain thickness liquid film, and liquid film has constituted the necessary electrolyte of metal erosion, and metal can corrode under this condition.Because corrosive medium dissolves in the airborne pollutant that can quicken metal erosion in the liquid film of metal surface, as airborne chloride and sulfate, therefore Important Project need be used a large amount of metals and coat of metal member usually, and different atmospheric environments has a significant impact the corrosion of metal and coat of metal member, therefore, measuring local atmospheric corrosiveness is one of capsule information of Important Project design.At present, atmospheric corrosiveness assessment both at home and abroad mainly contains two kinds of methods: standard metal method and environmental factor method.Wherein, the standard metal method is to be standard sample with aluminium, zinc, copper simple metal and Q235 steel, carries out on the spot 1 year by a definite date or the outdoor exposure of longer time, determines local atmospheric corrosiveness according to the annual corrosion rate of metal; The environmental factor method is according to SO in the atmosphere ₂Determine local atmospheric corrosiveness with year subsidence rate, year wetting time TOW (time of wetness) of Cl.Year wetting time in the environmental factor method be meant temperature in 1 year greater than 0 ℃, relative humidity greater than T.T. of 80%, wetting time will manually be read according to humiture diary curve usually.Existing method has following shortcoming:

(1) length consuming time needs the time more than 1 year

The standard metal method requires standard metal is carried out the exposure in 1 year by a definite date, and the zero-time of exposure required in corrosivity the strongest season; The environmental factor method requires to use the humiture registering instrument local humiture diary curve that draws uninterruptedly in 1 year.Therefore, rely on prior art assessment atmospheric corrosion, 1 year often not enough.

(2) need veteran professional on duty

Whole processes such as the removing of the collection of the preparation of the metal sample of standard metal method, the design and fabrication of test fixture, raw data, the daily management between exposure period, test back sample corrosion product, the collection of test back related data, experimental data processing all need the professional to be responsible for, and especially the removing of corrosion product requires operating personnel to possess rich experience.The environmental factor method need be settled humiture registering instrument, recording chart, needs draw temperature, moisture curve every day, needs every day and changes the humiture recording chart, and need proofread and correct and reading of data every day.This process also needs the special messenger and is responsible for.

(3) the conclusion confidence level that only obtains according to the detection data in 1 year is lower

The atmospheric corrosiveness of environment is in ceaselessly fluctuation, and the detection data in 1 year can not be assessed the variation tendency of atmospheric corrosiveness.And the data that prior art was only gathered 1 year are usually assessed, so the assessment result confidence level that prior art obtains is lower.

(4) can not obtain to support from meteorological station

There are nearly 700 meteorological stations in the whole nation, and all there is weather station each counties and cities, but they can not provide wetting time, and the humiture diary curve of every day can not be provided, so prior art can not make full use of the resource that the weather station obtains.

Summary of the invention

The assay method that the purpose of this invention is to provide a kind of inland non-industrial atmospheric corrosiveness.The Value Data when temperature that adopts the inventive method can make full use of observatory to record, humidity, the non-industrial atmospheric result who is surveyed is accurate, and method is simple.

Technical scheme of the present invention is:

1) determines the monitoring point, extract the nearest 10 years diary temperature of this area's meteorological station and the time Value Data of relative humidity;

2) humiture duration data preparation

Judge whether the temperature, the humidity data that are extracted have duplicate record, invalid record, disappearance record, and put temperature, the humidity data that is extracted in order, guarantee that temperature is consistent with the date of relative humidity data, calculates wetting time again;

3) the wetting time computing formula is as follows:

$\mathrm{TOW}=\underset{i=1}{\overset{365}{\mathrm{\Σ}}}\underset{j}{\mathrm{\Σ}}\mathrm{TOWd}(i,j),(j=\mathrm{2,8,14,20})...\left(1\right)$

In the formula

TOW is a year wetting time (h)

TOWd (i is 1 annual control i days j wetting times (h) constantly j), and its value is 6 or 0, when temperature＞0 ℃ and humidity＞80%, TOWd (i, j)=6, otherwise TOWd (i, j)=0, formula (1) is arranged in the Computing program;

4), calculate the average and the standard deviation of wetting time according to the above-mentioned gained wetting time value of each year;

5) determine wetting time grade according to average, determine the variation tendency of wetting time according to standard deviation;

6) determine the atmospheric corrosiveness grade of monitoring point according to wetting time grade.

Above-mentioned steps 2)～6) all in computing machine, handle.

The temperature of the monitoring over the years of the meteorological station above-mentioned steps 1) and relative humidity are the temperature and the relative humidity of every days 2 point, 8 points, 14 points, 20 records.

Described disappearance record is meant should have and do not have humiture data on certain period date.

The method of data preparation step 2) is: be criterion with the date, deletion duplicate record and invalid record, disappearance record is repaired, repaired principle and be: lack the time Value Data of 1,2 time point in one day, remedy according to the time Value Data of adjacent time point on the same day; The time Value Data that lacks within a couple of days is remedied according to the time Value Data of residing season, front and back a couple of days; The time Value Data of the more records of disappearance then directly calculates wetting time and will not remedy, and multiply by time coefficient τ, the actual fate in described time coefficient τ=365/ when calculating wetting time.

Atmospheric corrosiveness described in the step 4) is determined at the rate of corrosion of corresponding wetting time scope according to standard metal steel, Al, Zn, Cu.

Because every day repeatedly monitored local temperature and humidity by each meteorological station, accumulated the data of many decades.At present, meteorological station spreads all over the country, and the required distance that atmospheric corrosiveness detects only at a distance of tens of kilometers, can be satisfied fully in each monitoring point, and therefore, the present invention utilizes the data of meteorological station accumulation to record wetting time can solve the time-consuming difficult problem of taking a lot of work of prior art.

Suitable environment of the present invention: the outer and distance off sea of industrial pollution source 20Km is greater than the environment of 2Km.

SO in atmosphere ₂And Cl ^-Year subsidence rate respectively at 35mg/ (m ².d) and 60mg/ (m ².d) in the time of in, atmospheric corrosiveness is only relevant with wetting time.For SO ₂, have only usually industrial atmosphere year subsidence rate greater than 35mg/ (m ²And the shadow of industrial atmosphere is the zones of 20Km around the pollution source to scope .d); For Cl ^-, China's terrestrial environment has only the zone of limit 2Km off sea just may be higher than 60mg/ (m ².d).Therefore, the present invention is applicable to China overwhelming majority area.And, also can be and detect SO ₂And Cl ^-The atmospheric corrosiveness in seriously polluted area provides the wetting time data.

The applicant adopts the inventive method to experimentize on China's southeastern coast, southwest, northeast, northwest and other places, its result compares with current methods, less than 5%, Di Qu error is 5%～10% northeastward in the error of China's southeastern coast, southwest, and the Northwest is 10%～15%.

The explanation of nouns that the present invention relates to

Atmospheric corrosion

In atmospheric environment, metal, alloy, the coat of metal can be subjected to corrosion, the type of corrosion and degree depend on the electrolytical character that the wetted with product surface forms, it is generally acknowledged that the factor that influences electrolyte property has the kind of corrosive medium in the air and concentration and the corrosive medium time at surperficial continuous action.

The atmospheric corrosion grade

Characterize the severity of atmospheric corrosiveness, can divide, also can divide according to the rate of corrosion of standard metal outdoor exposure according to the size of wetting time and corrosive medium.

Standard metal

Be meant that chemical constitution meets the carbon steel of ISO9223 standard-required, fine aluminium, pure zinc, 4 kinds of metals of fine copper.

The ISO9223 standard

The international standard of assessment atmospheric corrosion.The corrosive medium that this standard is considered is SO ₂And Cl ^-, be wetting time TOW action time.

Sulphuric dioxide (SO ₂)

One of main corrosive medium of metal atmospheric corrosion, it has reflected the content of sulfate in the air, often represents with subsidence rate, unit is mg/ (100cm ².d).Discharged flue gas when being mainly derived from coal burning.

Chlorion (Cl ^-)

One of main corrosive medium of metal atmospheric corrosion, it has reflected muriatic content in the air, often represents with subsidence rate, unit is mg/ (100cm ².d).Cl in the marine atmosphere ^-Be mainly derived from seawater, Cl in the salt marsh district atmosphere ^-Be mainly derived from soil.

Critical temperature

Temperature is low more, and the metal atmospheric corrosion is slow more, when temperature is lower than a certain value, can think that metal atmospheric corrosion can not take place, and this temperature is critical temperature.The ISO9223 standard thinks that critical temperature is 0 ℃.

The necessary condition of critical relative moisture---metal atmospheric corrosion is that the metal surface can form certain thickness electrolyte liquid film, and when relative humidity during greater than a certain value, the electrolyte liquid film could form, and this relative humidity is critical relative moisture.The ISO9223 standard thinks that critical relative moisture is 80%.

Wetting time---temperature is higher than the T.T. that critical temperature, relative humidity are higher than critical relative moisture in the whole year, is wetting time.

Humiture diary curve---adopt the continuous temperature of drawing of humiture self-recording device and the curve of relative humidity real-time change.

Metal erosion rate---the physical quantity of characterizing metal extent of corrosion often represents that with the weightlessness that metal internal corrosion in a year causes unit is mg/ (m ².a) or mm/a.

Description of drawings

Fig. 1 is the inventive method process flow diagram;

Fig. 2 is wetting time computing method synoptic diagram;

Fig. 3 is type area's year wetting time.

Embodiment

The selection of the assay method of inland non-industrial atmospheric corrosiveness.

The critical temperature of wetting time and critical relative moisture are respectively 0 ℃ and 80%, overwhelming majority meteorological station is monitored temperature, the humidity of 4 times every day, time point is 2 points, 8 points, 14 points, 20 points, the data of adjacent 2 time points simultaneously greater than, simultaneously less than, 1 calculating greater than 1 wetting time during less than critical value:

Assumed temperature has designed 2 kinds computing method at 2 greater than 0 ℃, adjacent relative humidity as shown in Figure 2 in research process.Method one: the relative humidity that consideration A, B are ordered and the size of critical value difference, calculate wetting time in proportion, i.e. C ₁C ₄Method two: do not consider the relative humidity that A, B are ordered and the difference of critical value, to calculate wetting time, i.e. C 1/2nd of interval time ₁C ₃Among Fig. 2, the positive law of wetting time is C ₁C ₂, the error of method one is C ₂C ₄, its size was both relevant with positive law, and also the relative humidity of ordering with A, B is relevant with the difference between critical value, and the error of method two is C ₂C ₃, its size is only relevant with positive law.Since all multifactor influences, C ₂, C ₃, C ₄Relative position be not changeless, promptly the error of 2 kinds of methods has size mutually.Studies show that these 2 kinds of new method results and positive law do not have significant difference, easy in view of the 2nd kind of method calculating, the present invention adopts the 2nd kind of method to measure.

Assay method of the present invention is as follows: (referring to Fig. 1)

1) Value Data when local meteorological station is collected the humiture of every day over the years

The various places meteorological station all will be monitored local temperature and humidity every day, and writes down 2 points, 8 points, 14 points, 20 data.At least collect nearest 10 years data.

2) inavid date record deletion

May there be the data recording of inavid date in collected data, as February 30, April 31 etc.; Before calculating wetting time, it need be rejected.

3) repairing of missing data

Do not require when filling up very accurate, only need to judge temperature whether greater than 0 ℃ (critical temperature), relative temperature whether greater than 80% (critical relative moisture).In view of the above, a lot of missing datas can be remedied.

The data that lack 1,2 point in one day can be remedied according to other data on the same day; Data within a couple of days disappearance can be remedied according to the data of residing season, front and back a couple of days.When remedying when having no basis, just abandon remedying, at this moment when calculating wetting time, multiply by 1 time coefficient τ, the actual fate in τ=365/.

Data computation wetting time after employing is remedied should be indicated the fate of remedying.

4) rejecting of duplicate record

May have the record on repetition date from the data of meteorological station copy, unnecessary record should be rejected.

5) humiture data recording consistency check

Because missing data is difficult to remedy sometimes, thereby what cause temperature data record and relative humidity data recording is not corresponding one by one in time, and corresponding record is rejected or remedied to not corresponding record from data.

The wetting time computing formula is as follows:

$\mathrm{TOW}=\underset{i=1}{\overset{365}{\mathrm{\Σ}}}\underset{j}{\mathrm{\Σ}}\mathrm{TOWd}(i,j),(j=\mathrm{2,8,14,20})...\left(1\right)$

In the formula

TOW is a year wetting time (h)

TOWd (i is 1 annual control i days j wetting times (h) constantly j), and its value is 6 or 0, when temperature＞0 ℃ and humidity＞80%, TOWd (i, j)=6, otherwise TOWd (i, j)=0, formula (1) is arranged in the Computing program;

Calculate nearly 10 years annual wetting times, calculate its average and standard deviation again.

Above-mentioned steps all adopts computing machine to handle.

6) the atmospheric corrosion grade determines

After obtaining wetting time, can press the scope that the data of tabulation in the lattice are given the atmospheric corrosion rate of tapping, aluminium, zinc, copper.

Table 1 typical metal is at the rate of corrosion (μ m/a) of different wetting time

??TOW	??≤250h	??250～2500h	??2500～5500h	??＞5500h
					The TOW grade	??τ1、τ2	??τ3	??τ4	??τ5
The atmospheric corrosion grade	??C1	??C2	??C3	??C4
					??Q235	??≤1.3	??1.3～25	??25～50	??50～80
??Zn	??≤0.7	??0.7～2.1	??2.1～4.2	??4.2～8.4
					??Cu	??≤0.1	??0.1～0.6	??0.6～1.3	??1.3～2.8
??Al	??≤0.2	??0.2～0.7	??0.7～2	??2～3.7

For example: Bangbu calendar year 2001 TOW is 3324h, wetting time grade be 4 grades of τ on the upper side, in view of the above, the atmospheric corrosion grade in decidable Bangbu is the C3 level, and the rate of corrosion of typical metal is: Q235 25～50 μ m/a, Zn 2.1～4.2 μ m/a, Cu 0.6～1.3 μ m/a, Al 0.7～2 μ m/a.

With 5 places such as Chongqing, Haikous is example, and the wetting time in 1991 to 2000 in various places is seen Fig. 3, and the average of wetting time and standard deviation and corrosion class see Table 2.As seen the atmospheric corrosion grade in Haikou, Chongqing, Mo River, Lhasa, Dunhuang is respectively C4 level, C3 level, C2 level, C2 level, C1 level.

Table 2 type area atmospheric corrosion assessment result

The place	Average	Standard deviation	Wetting time grade	The atmospheric corrosion grade
					Mo River	??1796	??175	??τ3	??C2
Dunhuang	??181	??63	??τ2	??C1
					Lhasa	??809	??117	??τ3	??C2
Chongqing	??5096	??188	??τ4	??C3
					The Haikou	??5771	??200	??τ5	??C4

Embodiment

The example that is determined as with certain airport construction atmospheric corrosiveness describes the mensuration process that adopts the inventive method in detail:

1) meteorological station is selected: certain location, airport observatory, station index 57516.

2) nearest 10 annual datas

From observatory on Dec 31,1 day to 2005 January in 1996, Value Data when Value Data and relative humidity respectively has 3653 records during temperature.Only with in table 4, the table 5 surplus in the of 20 during humiture Value Data illustrate.

Table 4 57516 station temperature data examples

Sequence number	Station index	Date	2 points	8 points	14 points	20 points
							??1	??57516	??20050221	??59	??-52	??104	??-79
??2	??57516	??20050222	??68	??-41	??120	??97
							??3	??57516	??20050223	??64	??58	??100	??95
??4	??57516	??20050224	??80	??-50	??157	??140
							??5	??57516	??20050225	??110	??116	??114	??104
??6	??57516	??20050226	??98	??89	??100	??98
							??7	??57516	??20050227	??79	??73	??93	??90
??8	??57516	??20050228	??83	??79	??130	??122
							??9	??57516	??20050229	??9999	??9999	??9999	??9999
??10	??57516	??20050230	??9999	??9999	??9999	??9999
							??11	??57516	??20050231	??9999	??9999	??9999	??9999
??12	??57516	??2005031	??104	??95	??123	??112
							??13	??57516	??2005032	??100	??92	??140	??127
??14	??57516	??2005033	??113	??97	??8888	??101
							??15	??57516	??2005034	??97	??90	??8888	??128
??16	??57516	??2005036	??98	??80	??196	??170
							??17	??57516	??2005037	??125	??95	??200	??180
??18	??57516	??2005038	??146	??132	??222	??195
							??19	??57516	??2005039	??159	??138	??168	??168
??20	??57516	??20050310	??154	??157	??216	??200
							??21	??57516	??20050311	??167	??157	??172	??154

Table 5 57516 station relative humidity data instance

Sequence number	Station index	Date		2 points	8 points	14 points	20 points
								??1	??57516	??20050221	??81	??89	??53	??69
??2	??57516	??20050222	??79	??94	??54	??65
							??3	??57516	??20050223	??85	??89	??65	??68

??4	??57516	??2005024	??80	??93	??48	??59
							??5	??57516	??20050225	??73	??72	??79	??86
??6	??57516	??20050226	??92	??85	??73	??73
							??7	??57516	??20050227	??91	??91	??83	??86
??8	??57516	??20050228	??91	??93	??65	??70
							??9	??57516	??20050229	??999	??999	??999	??999
??10	??57516	??20050230	??999	??999	??999	??999
							??11	??57516	??20050231	??999	??999	??999	??999
??12	??57516	??2005031	??84	??90	??72	??75
							??13	??57516	??2005034	??90	??86	??25	??35
??14	??57516	??2005035	??76	??88	??29	??35
							??15	??57516	??20050225	??73	??72	??79	??86
??16	??57516	??20050226	??92	??85	??73	??73
							??17	??57516	??20050227	??91	??91	??83	??86
??18	??57516	??20050228	??91	??93	??65	??70
							??19	??57516	??2005036	??66	??78	??28	??40
??20	??57516	??2005037	??65	??85	??45	??47
							??21	??57516	??2005038	??67	??888	??888	??54
??22	??57516	??2005039	??71	??92	??888	??80
							??23	??57516	??20050310	??88	??85	??66	??71

3) data declaration

Table 4, table 5 are through insertion, deletion, replication processes, the data on this section date have just had duplicate record (as data as described in the sequence number in the table 5 15,16,17,18), inavid date record (as data as described in the sequence number in the table 49,10,11, sequence number 8,9,10 described data in the table 5), disappearance record (referring to lack all data recording all day), missing data defectives such as (Value Datas when referring to lack certain 1,2 in a day), above-mentioned defective is some intrinsic forms of data that meteorological station provides.

4) adopt the database language programmed process

Rough handling can be deleted duplicate record, inavid date record, and the character type date is converted to the type date on date, and all temperature datas are divided by 10.Temperature, the result behind the humidity data edit see Table 6 respectively, table 7.It should be noted that in deletion during repeating data, should confirm whether the legacy data record exists and effectively, delete in order to avoid miss.

Table 6 temperature data edit result

Table 7 humidity data edit result

Sequence number	Station index	Date	2 points	8 points	14 points	20 points
							??1	??57516	??2005-2-21	??81	??89	??53	??69
??2	??57516	??2005-2-22	??79	??94	??54	??65
							??3	??57516	??2005-2-23	??85	??89	??65	??68
??4	??57516	??2005-2-24	??80	??93	??48	??59
							??5	??57516	??2005-2-25	??73	??72	??79	??86
??6	??57516	??2005-2-26	??92	??85	??73	??73
							??7	??57516	??2005-2-27	??91	??91	??83	??86
??8	??57516	??2005-2-28	??91	??93	??65	??70
							??9	??57516	??2005-3-1	??84	??90	??72	??75
??10	??57516	??2005-3-4	??90	??86	??25	??35
							??11	??57516	??2005-3-5	??76	??88	??29	??35
??12	??57516	??2005-3-6	??66	??78	??28	??40
							??13	??57516	??2005-3-7	??65	??85	??45	??47
??14	??57516	??2005-3-8	??67	??888	??39	??54
							??15	??57516	??2005-3-9	??71	??92	??888	??80
??16	??57516	??2005-3-10	??88	??85	??66	??71

5) inspection of disappearance record

Adopt the database language programming to check.Implementation step is: the criterion date is set, opens and select tables of data, on the date of the 1st record of reading of data table, relatively whether record date equated with the criterion date.If equate, illustrate that data recording is normal, the criterion date adds 1 day then, simultaneously the date of the 2nd record of reading of data table, compares judgement again; If the criterion date, has illustrated the record disappearance less than record date, the criterion date is exactly the record of disappearance to the record between the record date, disappearance information is noted, compose after record date added 1 day and give the criterion date, next bar record of reading of data table compares judgement more simultaneously; If the criterion date greater than record date, has illustrated duplicate record, this situation is the 4th) handle in the step.The rest may be inferred, can find out all the disappearance records in the selected data table.The temperature data and the humidity data that extract from observatory all will carry out this inspection.

Check to be example explanation checking process with table 6 temperature data record disappearance.If the criterion date is on February 21st, 2005, the date of reading in the 1st record in the table 6 is on February 21st, 2005, equate with the criterion date, illustrate that data recording is normal, the criterion date is on February 22nd, 2005 after adding 1 day, the 2nd record date of reading of data table is on February 22nd, 2005, equates with the criterion date.So repeatedly, the criterion date is on March 5th, 2005 during to the 13rd record, record date is on March 6th, 2005, and record date is bigger by 1 than the criterion date, and 1 record of disappearance is described, disappearance is recorded as on March 5th, 2005, behind disappearance information record, record date is added 1 day promptly composed on March 7th, 2005 that the date of reading the 14th record again is on March 7th, 2005 to the criterion date, the two equates, illustrates that data recording is normal.The rest may be inferred, can find out all the disappearance records in the table 6.The humidity data of table 7 all will carry out this inspection.

6) repairing of disappearance record

If the front and back record that shows the disappearance record the record data of disappearance a couple of days, can be remedied according to the data of residing season, front and back a couple of days.As lacking the thermograph on March 5th, 2005 in the table 6, and the temperature before and after March 5 all is higher than 0 ℃, therefore, and can be with the temperature on March 4 temperature as March 5.The humidity data that and for example lacks March 2,3 in 2005 in the table 7, and point, 8 humidity all were higher than 80%, 14 point, at 20 and all were lower than 80% on March 1,42, therefore, can be with the humidity on March 1 humidity as March 2,3.When being difficult to determine, that is not just repaired, multiply by time coefficient at last and get final product, this situation with year average repair.

7) inspection of missing data in a day

Temperature in one day, humidity raw data can be designated as 8888 and 888 respectively usually if disappearance is arranged, and are respectively 888.8 and 888 behind the edit, but blank, nonnumeric situation also can occur sometimes.For this reason, the criterion of setting when checking missing data can be: 1. numerical value shape, non vacuum data; 2. temperature range can be set to-60 ℃～60 ℃, and humidity range is 0～100%, according to above-mentioned 1., 2. whether described criterion inspection has the data of disappearance in one day.

8) repairing of disappearance intermediate data

Because when repairing data, whether temperature data only judge more than or equal to 0 ℃, and therefore whether humidity data is only judged, often can repair scarce data according to the front and back data more than or equal to 80%.As in the table 6 March 3,4 in 2005 14 temperature missing data decidable greater than 0 ℃, and should be about 14 ℃; The humidity missing data in March 8,9 in 2005 in the table 7 and for example, according to the existing humidity data decidable in March 1 to March 10: March 88 relative humidity to be higher than 80% probability very big, March 9,14 relative humidity was lower than 80% certainly.Judge in view of the above and can repair scarce data.

9) when carrying out wetting time calculating, must guarantee that the date in table 6, the table 7 is corresponding one by one, as: the date of the 5th record is on February 25th, 2005 in the table 6, and the date of the 5th record also is necessary on February 25th, 2005 in the table 7, and the rest may be inferred.And writing down when repairing, may there be the situation that can't repair, thereby causes the date of each record in two tables correspondence that differs, therefore, must carry out the date and check one to one, another table has if certain table lacks a record, and the respective record of then deleting another table gets final product.After all record checks finished, rechecking once guaranteed that respective record is corresponding one by one.

10) carry out duplicate record deletion, inavid date record deletion, disappearance record and data modification, one by one correspondence proving with repair after, can carry out the calculating of wetting time.Temperature, the final reduced data of humidity see Table 8 respectively, table 9, carry out the calculating of wetting time according to the data of formula of the present invention (1) his-and-hers watches 8, table 9.

Table 8 temperature data is finally put the result in order

Table 9 humidity data is finally put the result in order

Sequence number	Station index	Date	2 points	8 points	14 points	20 points
							??1	??57516	??2005-2-21	??81	??89	??53	??69
??2	??57516	??2005-2-22	??79	??94	??54	??65
							??3	??57516	??2005-2-23	??85	??89	??65	??68
??4	??57516	??2005-2-24	??80	??93	??48	??59
							??5	??57516	??2005-2-25	??73	??72	??79	??86
??6	??57516	??2005-2-26	??92	??85	??73	??73
							??7	??57516	??2005-2-27	??91	??91	??83	??86
??8	??57516	??2005-2-28	??91	??93	??65	??70
							??9	??57516	??2005-3-1	??84	??90	??72	??75
??10	??57516	??2005-3-1	??84	??90	??72	??75
							??11	??57516	??2005-3-1	??84	??90	??72	??75
??12	??57516	??2005-3-4	??90	??86	??25	??35
							??13	??57516	??2005-3-5	??76	??88	??29	??35
??14	??57516	??2005-3-6	??66	??78	??28	??40
							??15	??57516	??2005-3-7	??65	??85	??45	??47
??16	??57516	??2005-3-8	??67	??85	??39	??54
							??17	??57516	??2005-3-9	??71	??92	??39	??80
??18	??57516	??2005-3-10	??88	??85	??66	??71
							??19	??57516	??2005-3-11	??88	??85	??66	??71

11) calculate wetting time

As: the data on February 21st, 2005 (sequence number 1) in table 8, the table 9, greater than 80%, wetting time then adds 6h greater than 0 ℃ of relative humidity for 2 temperature, and 8 temperature is less than 0 ℃, disregard wetting time, 14 temperature less than 80%, are disregarded wetting time greater than 0 ℃, relative humidity, and 20 temperature is less than 0 ℃, also disregard wetting time, therefore, on February 21st, 2005, the wetting time on certain airport is 6h.The rest may be inferred, calculates annual wetting time.

12) can calculate according to the method described above, the wetting time of each year

The year wetting time (h) in 1996 to 2005 in certain airport is respectively 5298,5010,4830,5274,5214,4764,5184,5118,4638,4698.

13) calculate the average and the standard deviation of year wetting time according to the wetting time of each year, standard deviation calculation formulas is:

$\sqrt{\frac{\mathrm{n\Σ}{x}^2-{\left(\mathrm{\Σx}\right)}^2}{n(n-1)}}$

N is the quantity in the year of mensuration wetting time in the formula, and x is the mean value of wetting time.

Above-mentioned steps all adopts computing machine to handle.

The average and the standard deviation of certain airport wetting time are respectively: 5002h, 250h.

14) wetting time and atmospheric corrosiveness grade are judged

The probable range of knowing certain airport wetting time 99.7% by the average and the standard deviation of wetting time is 4252h～5752h, knows according to table 1, table 2, and certain airport wetting time grade is 4 grades of upper limits of τ and τ 5 lower limits, and the atmospheric corrosiveness grade is the C3 level.The rate of corrosion that further can determine each standard metal should be near the higher limit of C3 level.

Claims (6)

1. the assay method of an inland non-industrial atmospheric corrosiveness is characterized in that following steps are arranged:

1) determines the monitoring point, extract the nearest 10 years diary temperature of this area's meteorological station and the time Value Data of relative humidity;

2) temperature, humidity data arrangement

Judge whether the temperature, the humidity data that are extracted have duplicate record, invalid record, disappearance record, and put temperature, the humidity data that is extracted in order, guarantee that temperature is consistent with the date of relative humidity data, calculates wetting time again;

3) the wetting time computing formula is as follows:

$\mathrm{TOW}=\underset{i=1}{\overset{365}{\mathrm{\Σ}}}\underset{j}{\mathrm{\Σ}}\mathrm{TOWd}(i,j)$ (j＝2、8、14、20)................(1)

In the formula

TOW is a year wetting time (h)

TOWd (i is 1 annual control i days j wetting times (h) constantly j), and its value is 6 or 0, when temperature＞0 ℃ and humidity＞80%, TOWd (i, j)=6, otherwise TOWd (i, j)=0;

4), calculate the average and the standard deviation of wetting time according to the above-mentioned gained wetting time value of each year;

5) determine wetting time grade according to average, determine the variation tendency of wetting time according to standard deviation;

6) determine the atmospheric corrosiveness grade of monitoring point according to wetting time grade.

2. the assay method of inland non-industrial atmospheric corrosiveness according to claim 1 is characterized in that the time Value Data described in the step 1) is every days 2 point, 8 points, 14 points, 20 temperature and relative humidity data.

3. the assay method of inland non-industrial atmospheric corrosiveness according to claim 1 is characterized in that: step 2) described in disappearance record be meant have and do not have humiture data on certain or certain period date.

4. the assay method of inland non-industrial atmospheric corrosiveness according to claim 1, it is characterized in that step 2) described in the method for data preparation be: be criterion with the date, deletion duplicate record and invalid record, the disappearance record is repaired, the repairing principle is: for the temperature data of disappearance need judge its whether more than or equal to 0 ℃, judge then that for the humidity data of disappearance whether it is greater than 80%, if lack the time Value Data of 1,2 time point in one day, remedy according to the time Value Data of adjacent time point on the same day; The time Value Data that lacks within a couple of days is remedied according to the time Value Data of residing season, front and back a couple of days; The time Value Data of the more records of disappearance then directly calculates wetting time and will not remedy, and multiply by time coefficient τ when calculating wetting time.

5. the assay method of inland non-industrial atmospheric corrosiveness according to claim 4 is characterized in that: the actual fate in described time coefficient τ=365/.

6. the assay method of inland non-industrial atmospheric corrosiveness according to claim 1 is characterized in that the atmospheric corrosiveness described in the step 4) determines at the rate of corrosion of corresponding wetting time scope according to standard metal steel, Al, Zn, Cu.

Images