CN106777830A - A kind of fast appraisement method of pipeline hydraulic frictional resistance - Google Patents

A kind of fast appraisement method of pipeline hydraulic frictional resistance Download PDF

Info

Publication number
CN106777830A
CN106777830A CN201710072656.8A CN201710072656A CN106777830A CN 106777830 A CN106777830 A CN 106777830A CN 201710072656 A CN201710072656 A CN 201710072656A CN 106777830 A CN106777830 A CN 106777830A
Authority
CN
China
Prior art keywords
value
test
appraisal
sample length
length
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
CN201710072656.8A
Other languages
Chinese (zh)
Other versions
CN106777830B (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN201710072656.8A priority Critical patent/CN106777830B/en
Publication of CN106777830A publication Critical patent/CN106777830A/en
Application granted granted Critical
Publication of CN106777830B publication Critical patent/CN106777830B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Sampling And Sample Adjustment (AREA)

Abstract

The present invention relates to a kind of fast appraisement method of pipeline hydraulic frictional resistance, including:Preresearch estimates;Choose test and appraisal sample length and test and appraisal evaluation length;Test and appraisal detection;Whether assessment sample length is in proper range;Select shorter sample length;Determine sample length;General measure;Calculate average value;Calculate hydraulic simulation experiment.Surface roughness measuring instrument during method of the present invention is processed using existing machinery measures length by adjustment repeatedly, and determines the maximum height of the profile of surface roughnessR z Can be as equivalent roughnessk, bring Ke Erboluoke-White's formula into and calculate hydraulic simulation experimentλ.Make the measurement of hydraulic simulation experiment and calculate more convenient, fast, simplified, improve pipeline construction quality, reduce the operating cost of pipeline engineering.

Description

A kind of fast appraisement method of pipeline hydraulic frictional resistance
Technical field
It is a kind of assessment method of water conservancy project parameter the present invention relates to a kind of fast appraisement method of pipeline hydraulic frictional resistance, is A kind of method to the test and appraisal of hydraulic pipe inwall frictional resistance.
Background technology
The pipeline hydraulic coefficient of friction resistance is one of important parameter of fluid conveying engineering design, and the result of its value is direct Influence the precision of water force achievement.The accurate measurement and reasonable value of hydraulic simulation experiment are for long distance water transfer project, city City's water supply and sewage work, agricultural irrigation engineering, and the total arrangement of the pipeline engineering such as oil transportation, design scale, lectotype selection, Operational management and energy-saving and emission-reduction have great importance.
The conventional hydraulic calculation formula of engineering is darcy-Wei Si Baha(Darcy-Weisbach)Formula,
In formula:h f - And Research For The Loss Head Along Pipeline(m);λ- hydraulic simulation experiment, or friction factor of head loss;DIn-pipeline Footpath (m);L- length of pipe section(m);V- pipeline mean flow rate(m/s);g- acceleration of gravity(m/s2).
Hydraulic simulation experiment λ is calculated and is generally used Ke Erboluoke-White(Colebrook-White)Formula
In formula:k- equivalent roughness (m);Re- Reynolds number, computing formula isνThe kinematic coefficient of viscosity of-water (m2/s)。
The applicable reynolds number range of the formula is wide, and it is good to be coincide with the resistance test result of actual commercialization pipeline, by engineering Master's master is considered to calculate the normalized form of hydraulic simulation experiment.
Equivalent roughness in Ke Shi formulakThe shadow of concentrated expression inner-walls of duct rough grain and various factors to linear loss Ring, all the time the main linear loss that pipeline under high reynolds number is determined by fluid test, and by hydraulic simulation experimentλConversion Be the rough grain of Ni Gulazi rough pipes of same diameter highly, it is believed that rough grain height is flat with the irregular coarse infinitesimal of inner-walls of duct It is closely sized to, this is equivalent roughnessk, but fluid test needs fixed experimental site and detection platform, and expend more Time and expense, be not suitable for the detection of large diameter pipeline.
The content of the invention
In order to overcome problem of the prior art, the present invention to propose a kind of fast appraisement method of pipeline hydraulic frictional resistance.Institute The various contact pin type surface roughness instruments devices that the method stated is used in being processed using existing machinery, by the measurement transformed Journey is measured to the roughness of inner-walls of duct, the quick hydraulic simulation experiment for obtaining pipeline.
The object of the present invention is achieved like this:A kind of fast appraisement method of pipeline hydraulic frictional resistance, the step of methods described It is rapid as follows:
The step of preresearch estimates:Estimated value for estimating the arithmetic average deviation of the roughness profile of tested pipelineR aEitimate 、 The estimated value of the maximum height of roughness profileR zEitimate
The step of choosing test and appraisal sample length and test and appraisal evaluation length:For basisR aEitimate R zEitimate Or roughness profile Arithmetic average deviation test and appraisal valueR ae With the test and appraisal value of the maximum height of roughness profileR ze Choose test and appraisal sample lengthlr e Or Test and appraisal measurement lengthln e
The step of test and appraisal detection:For according to test and appraisal sample lengthlr e Tested pipeline inwall is surveyed along tested pipeline axis Detection is commented, one group is obtainedR ae R ze Value;
Assessment sample length whether the step in proper range:For to the groupR ae WithR ze Numerical value be estimated, its numerical value Whether in test and appraisal sample lengthlr e Or test and appraisal measurement lengthln e The maximum height value of corresponding roughness profileR a And roughness The arithmetic average deviation value of profileR z In the range of value;The group is determined if "Yes"R ae WithR ze Numerical value be pre-selectionR a WithR z Value, and enter next step, returned to if "No" " the step of choosing test and appraisal sample length and test and appraisal evaluation length ";
The step of selecting shorter sample length:For choosing shorter test and appraisal sample lengthlr s , and one group of roughness profile of acquisition The shorter test and appraisal value of arithmetic average deviationR as With the shorter test and appraisal value of the maximum height of roughness profileR zs , described shorter survey Comment sample lengthlr s Refer to:Than in the most short survey selected by " the step of choosing test and appraisal sample length and test and appraisal evaluation length " Sample length is commented also to want short test and appraisal length;
The step of determining sample length:IfR as R zs Numerical value meet corresponding sampling scope, it is determined thatR as R zs It is corresponding Shorter test and appraisal sample lengthlr s Sample length used in general measurelr, otherwise select " assessment sample length be The no step in proper range " pre-selectionR a R z The corresponding test and appraisal sample length of valuelr e Taking used in general measure Sample lengthlr
The step of general measure:For using sample lengthlrCollection tested pipeline inwall multiple positionR ai WithR zi Value, its In:i=1,2 ... ...,n,nIt is positive integer;
The step of calculating average value:For calculatingR ai WithR zi Value average value:
,
The step of calculating hydraulic simulation experiment:For by measureAskValue is substituted into Ke Shi formula:
,
The hydraulic simulation experiment λ of tested pipeline is calculated, and draws corresponding hydraulic calculation chart;
In formula:λ- hydraulic simulation experiment;D- internal diameter of the pipeline;k- equivalent roughness, usesInstead of;Re- Reynolds number, meter Calculating formula isνThe kinematic coefficient of viscosity of-water;V- pipeline mean flow rate;
DescribedR aEitimate R zEitimate R an R zn R ae R ze R as R zs R ai R zi It is the different phase in detection process 'sR a WithR z , changed to show difference with footnote;Describedlr e ln e lr s It is the different phase in detection processlrln, use Footnote changes to show difference.
Further, described sample lengthlrOr measurement lengthlnCorrespondingR a WithR z Scope it is as follows:
R a ≤ 10 microns,lrIt is 2.5 millimeters,lnIt is 12.5 millimeters;10 microns of <R a ≤ 80 microns,lrIt is 8 millimeters,lnIt is 40 millis Rice;
R z ≤ 50 microns,lrIt is 2.5 millimeters,lnIt is 12.5 millimeters;50 microns of <R z ≤ 200 microns,lrIt is 8 millimeters,lnIt is 40 Millimeter.
Further, the step of described calculating average value after set:The step of calculating standard variance and judgement are managed Road quality step:
The step of calculating standard variance:
,
The step of judging pipeline quality:For being calculated by formulaR a WithR z Coefficient of dispersionCV Ra CV Rz
,
And judgeCV Ra CV Rz Whether coefficient of dispersion is less thanCVThreshold value, illustrated if "Yes" measured pipeline be it is uniform, Otherwise illustrate the non-material liner of the same race of measured pipe fitting, or the partial tube liner quality of production is not up to standard.
Further, described coefficient of dispersionCVThreshold value be:0.3.
The beneficial effect comprise that:Method of the present invention uses the surface roughness in existing machinery processing Measuring instrument determines the maximum height of the profile of surface roughness by adjustment measurement length repeatedlyR z Can be thick as equivalent Rugosityk, bring Ke Erboluoke-White into(Colebrook-White)Formula calculates hydraulic simulation experimentλ.Make flow resistance system Several measurements and calculating more facilitate, simplify, and improve construction quality, reduce the operating cost of hydraulic pipe.
Brief description of the drawings
The invention will be further described with reference to the accompanying drawings and examples.
Fig. 1 is the flow chart of the methods described of embodiments of the invention one;
Fig. 2 is the maximum height value of roughness profileR z Schematic diagram;
Fig. 3 is the flow chart of the methods described of embodiments of the invention three.
Specific embodiment
Embodiment one:
The present embodiment is a kind of fast appraisement method of pipeline hydraulic frictional resistance.The basic ideas of methods described are by machining 's(Machining)Surface roughness, is incorporated into the inner wall surface detection of water pipeline, as the pipeline hydraulic coefficient of friction resistancekMake With.
Existing method be withR q Value represents equivalent roughness value, and givesk≈5R q k≈3R q k≈1.6R q Deng not Same calculation.
The present embodiment has then used the maximum height value of roughness profileR z , makeR z k.Table 1 is givenR a R z R q Experiment Result is compareed.From 1-5 groups data in table,R z k, with 5R q Compare,R z Value is closerkValue, andR z Definition and equivalent RoughnesskPhysical significance it is closer, can more represent influence of the inner-walls of duct degree of roughness to current friction loss.
The result of the test of the surface roughness parameter of table 1 and equivalent roughness k relations
The equipment that the detection of surface roughness is used, it is possible to use contact pin type surface roughometer, the instrument is generally used for machine Tool processes the surface roughness of vacuum metrics design of part(Or surface microscopic unevenness).Contact pin type surface roughometer is using golden steel Stone contact pilotage carrys out the roughness on tracking measurement surface, rides on workpiece contact pilotage during measurement, and measured surface perpendicular contact, using drive With certain speed tow sensor, because measured surface profile peak valley rises and falls, contact pilotage will be produced dynamic device when measured surface is slided Life is moved up and down, and the length travel of contact pilotage is converted into electric signal, and this movable signal is amplified by electronic installation, To be felt the roughness found so as to measure human eye, this kind of equipment has good lateral resolution, and accuracy of detection can reach To secondary nanophase.
The major parameter of its measurement includes:
Sample length(Sampling length)lr:For differentiating the length being assessed on the X axis of the atypical characteristic of profile Degree.Regulation and selection sample length are the influences in order to limit and weakening surface waviness to the measurement result of surface roughness, Surface is more coarse, and sample length is bigger.
Evaluation length(evaluation length)ln:For differentiating the length being assessed in the X-direction of profile, can Comprising one or several sample lengths.It is past on a sample length because the roughness of surface each several part is not necessarily very uniform It is past reasonably to reflect a certain surface roughness characteristics, therefore need to several sample lengths be taken on the surface to evaluate rough surface Degree.
The maximum height of roughness profile(maximum height of roughness profile)R z :Taken at one In sample length, largest contours peak heightZ p With largest contours paddy depthZ v The height of sum, as shown in Figure 2.
In Fig. 2mIt is mean profile line, with mean profile line as line of demarcation, measures crest and trough.In Fig. 2Z P1~Z P6 It is six crests, wherein maximum crest isZ P6, as roughness concentration lengthlrInteriorZ Pmax ,Z V1~Z V6It is six troughs, its Middle maximum trough isZ V6, as roughness concentration lengthlrInteriorZ Vmax
The arithmetic average deviation of roughness profile(arithmetical meandeviation of roughness profile)R a :The arithmetic mean of instantaneous value of ordinate value Z (x) absolute value in a sample length
The present embodiment is by roughness profile maximum heightR z ≈ equivalent roughnessk, can be by practiceR z Value is used as equivalent roughnessk Value substitutes into Ke Shi formula:
In,
To the hydraulic simulation experiment of pipelineλCalculated and evaluated.
The present embodiment methods described is comprised the following steps that:
(One)The step of preresearch estimates:Estimated value for estimating the arithmetic average deviation of the roughness profile of tested pipelineR aEitimate , roughness profile maximum height estimated valueR zEitimate .The essential step of inner-walls of duct roughness measurement is true Determine sample length, or measurement length.The detection of pipe roughness and the measurement of machining surface roughness, its is important Difference is, the formation of machining surface be it is relevant with process, for example:Vehicle Processing, planing operation etc. are rough machined Cutter trade is fairly obvious, and the cutter trade of the retrofit such as grinding is finer and closely woven, and both very easily distinguish, because This, how as a rule need not be excessive take notice of determines sample length and measurement length, directly with the naked eye visually can be really Determine sample length or measurement length.Test for inner-walls of duct roughness is then entirely different.Mainly due to the material of pipeline The differences such as matter, processing mode, and form the diversity of inner-walls of duct roughness.Accordingly, it would be desirable to measurement length or sampling length Degree carries out assessment at least one times, to find appropriate measurement length or sample length.
This step is the initial of measurement roughness, first according to the surface roughness for visually observing tested pipeline, according to sight The result examined, estimates the roughness profile parameter of test specimenR a R z Numerical value, for the ease of difference, this step is thick by thisR a R z Numerical value be defined as roughness profile arithmetic average deviation estimated valueR aEitimate , roughness profile maximum height Estimated valueR zEitimate
(Two)The step of choosing test and appraisal sample length and test and appraisal evaluation length:For basisR aEitimate R zEitimate Or it is thick The test and appraisal value of the arithmetic average deviation of rugosity profileR ae With the test and appraisal value of the maximum height of roughness profileR ze Choose test and appraisal sampling Lengthlr e Or test and appraisal measurement lengthln e
This step is to choose test and appraisal sample length in both caseslr e Or test and appraisal measurement lengthln e .A kind of situation is just Begin measurement, in the case of initial measurement, using according to a preliminary estimateR a R z Value, carrys out lower selection test and appraisal sample lengthlr e Or test and appraisal Measurement lengthln e .ShouldR a R z It is defined as in previous step:R aEitimate R zEitimate .Also a kind of situation is:Have been carried out Preliminary surveying is measured several times, but without the requirement for reaching determination sample length or measurement length, also further to be adjusted It is whole.In this case, the simply experiment test value for being used, without what is be to determineR a R z Value, for ease of difference, the present embodiment Define the test and appraisal value of the arithmetic average deviation that the value is roughness profileR ae With the test and appraisal value of the maximum height of roughness profileR ze , to show difference.Use simultaneouslyR aEitimate R zEitimate Value orR ae R ze Value is identified to measure length and sample length also not The measurement length that must be used in the formal measurement that must be futurelrAnd sample lengthln, taken using test and appraisal for the sake of difference Sample lengthlr e Or test and appraisal measurement lengthln e , to show difference.
To make measurement reach certain precision, it is necessary to select appropriate measurement lengthlnOr sample lengthlr.AndlnlrReally Surely it is basisR a R z The size of value, andR a R z The size of value exactly needs the numerical value of measurement, therefore, accurately to be measuredR a R z Numerical value, first has to experimentally be measured, and finds optimal sample lengthlnOr measurement lengthlr
(Three)The step of test and appraisal detection:For according to test and appraisal sample lengthlr e Along tested pipeline axis in tested pipeline Wall carries out test and appraisal detection, obtains one groupR ae R ze Value.
The measurement of this step can be that roughmeter enters detection using probe.By roughmeter probe along conduit axis level Place, make measurement direction parallel to tested pipeline inwall axial direction, press(Two)The test and appraisal sample length of middle pre-selectionlr e Complete once pre- Measurement.
(Four)Assessment sample length whether the step in proper range:For to the groupR ae WithR ze Numerical value commented Estimate, whether its numerical value is in test and appraisal sample lengthlr e Or test and appraisal measurement lengthln e The maximum height value of corresponding roughness profileR a With the arithmetic average deviation value of roughness profileR z In the range of value;The group is determined if "Yes"R ae WithR ze Numerical value for pre-selection 'sR a WithR z Value, and enter next step, returned to if "No" and " choose the step of test and appraisal sample length and test and appraisal evaluation length Suddenly ".
By what is measuredR ae WithR ze Numerical value, withR a WithR z Measurement length in the sampling scope and the measurement range table of comparisons of value With corresponding to sample lengthR a WithR z Number range compare.If measured value is beyond the corresponding number of pre-selection sample length Value scope, then should be set by the corresponding sample length of measured value, i.e., instrument is adjusted to corresponding higher or lower sampling Length.Then measure one group of parameter value using the sample length of this adjustment, and again withR a WithR z The sampling scope of value and measurement In the scope table of comparisons numeric ratio compared with.This process can be with repeated multiple times, to reachR a WithR z Value sampling scope and measurement range The combination of the measured value and sample length of control tabular value suggestion.
(Five)The step of selecting shorter sample length:For choosing shorter test and appraisal sample lengthlr s , and one group of acquisition is coarse Spend the shorter test and appraisal value of the arithmetic average deviation of profileR as With the shorter test and appraisal value of the maximum height of roughness profileR zs , it is described Shorter test and appraisal sample lengthlr s Refer to:Than selected by " the step of choosing test and appraisal sample length and test and appraisal evaluation length " Most short test and appraisal sample length also wants short test and appraisal length;
By after a test and appraisal detection twice, having had certain ranging concept to measurement length or sample length, but to make survey Amount is more accurate, can also be to sample lengthlrOr measurement lengthlnSelect again once.This step was selected in former survey Commenting in step does not have to use shorter sample lengthlr s , then sample length is adjusted to one group of more shorter acquisitionR a R z Numerical value, This step is defined as:R ae R ze , to show difference.Measured by checkingR ae R ze Numerical value and sample length combination it is whether also full FootR a WithR z The sampling scope of value and the regulation of the measurement range table of comparisons.
(Six)The step of determining sample length:IfR as R zs Numerical value meet corresponding sampling scope, it is determined thatR as R zs Corresponding shorter test and appraisal sample lengthlr s Sample length used in general measurelr, otherwise select in " assessment sampling Length whether the step in proper range " pre-selectionR a R z The corresponding test and appraisal sample length of valuelr e By making in general measure Sample lengthlr
Up to the present two groups have been obtainedR a R z Value, one group is:R as R zs Numerical value, another group isR ae R ze Number Value.As long asR as R zs Numerical value meet sampling scope and the measurement range table of comparisons, then select:R as R zs Correspondinglr s ln s For in general measurelrOrln.Otherwise just selectR ae R ze Corresponding to valuelr e ln e For in general measurelrOrln
(Seven)The step of general measure:For using sample lengthlrCollection tested pipeline inwall multiple positionR ai WithR zi Value, wherein:i=1,2 ... ...,n,nIt is positive integer.
With the sample length selected in above-mentioned stepslrThe multiple measurement of different parts is carried out to being detected inner-walls of duct, it is complete Into the detection of required parameter, obtain multigroupR ai WithR zi Value.
(Eight)The step of calculating average value:For calculatingR ai WithR zi Value average value:
,
(Nine)The step of calculating hydraulic simulation experiment:For by measureAskValue is substituted into Ke Shi formula:
,
The hydraulic simulation experiment λ of tested pipeline is calculated, and draws corresponding hydraulic calculation chart;
In formula:λ- hydraulic simulation experiment;D- internal diameter of the pipeline;k- equivalent roughness, usesInstead of;Re- Reynolds number, meter Calculating formula isνThe kinematic coefficient of viscosity of-water;V- pipeline mean flow rate.
DescribedR aEitimate R zEitimate R an R zn R ae R ze R as R zs R ai R zi It is in detection process different StageR a WithR z , changed to show difference with footnote;Describedlr e ln e lr s ln s It is the different phase in detection processlrln, changed to show difference with footnote.
Embodiment two:
The present embodiment is the improvement of embodiment one, is embodiment one on sample lengthlrOr measurement lengthlnCorrespondingR a WithR z Scope refinement.Sample length described in the present embodimentlrOr measurement lengthlnCorrespondingR a WithR z Scope it is as follows:
R a ≤ 10 microns,lrIt is 2.5 millimeters,lnIt is 12.5 millimeters;10 microns of <R a ≤ 80 microns,lrIt is 8 millimeters,lnIt is 40 millis Rice;
R z ≤ 50 microns,lrIt is 2.5 millimeters,lnIt is 12.5 millimeters;50 microns of <R z ≤ 200 microns,lrIt is 8 millimeters,lnIt is 40 Millimeter.
It is expressed as with form:
Table 2R a The roughness sample length of value
Table 3R z The roughness sample length of value
Table 4 and table 5 are the specification of machining surface roughness, CNS(GB/T 10610 2009)Middle recommendation Roughness sample length.The selection of sample length has large effect to roughness parameter measurement result.6-9 groups data in table 1R a > 10umOrR z > 50um, now sample length should belr=8.0mm, and its sample length is stilllr=2.5mm, sample lengthlrIt is less than normal to cause pipelineR z Value is less thankThe reason for value.Equally, sample length is too small, such aslr≤ 0.8mm, measuresR z Value is inclined It is small, can not well represent influence of the inner wall roughness to flow turbulence.
During table 4 is machinedR a The roughness sample length of value
During table 5 is machinedR z The roughness sample length of value
Embodiment three:
The present embodiment is the improvement of above-described embodiment, is above-described embodiment on using calculatingR a R z Value is evaluated in pipeline The method of wall quality, flow is as shown in Figure 3:
Set after the step of described calculating average value:Calculate standard variance the step of and judge pipeline quality step:
The step of calculating standard variance:
,
The step of judging pipeline quality:For being calculated by formulaR a WithR z Coefficient of dispersionCV Ra CV Rz
,
And judgeCV Ra CV Rz Whether coefficient of dispersion is less thanCVThreshold value, illustrated if "Yes" measured pipeline be it is uniform, Otherwise illustrate the non-material liner of the same race of measured pipe fitting, or the partial tube liner quality of production is not up to standard.
The present embodiment is by comparing the coefficient of dispersion of measurement result(The ratio of standard deviation and its average), make Coefficient of dispersion CV < 0.3, so that it is guaranteed that the uniformity of measured pipeline, is otherwise considered as in the non-material of the same race of measured pipe fitting Lining, or the partial tube liner quality of production is not up to standard.The material and quality of hydraulic pipe are may determine that with this kind of mode, to improve Construction quality.
Example IV:
The present embodiment is the improvement of above-described embodiment, is refinement of the above-described embodiment on coefficient of dispersion.Described in the present embodiment Coefficient of dispersionCVThreshold value be:0.3.
This gives the roughness sample length suitable for evaluating the pipeline hydraulic coefficient of friction resistance(Table 2, table 3), with And measurement result evaluation criterion(Coefficient of dispersion CV < 0.3).
Finally it should be noted that being merely illustrative of the technical solution of the present invention and unrestricted above, although with reference to preferable cloth Scheme is put to be described in detail the present invention, it will be understood by those within the art that, can be to technology of the invention Scheme(The type selecting of such as measuring instrument, the utilization of various formula, sequencing of step etc.)Modify or equivalent, Without deviating from the spirit and scope of technical solution of the present invention.

Claims (4)

1. a kind of fast appraisement method of pipeline hydraulic frictional resistance, it is characterised in that as follows the step of methods described:
The step of preresearch estimates:Estimated value for estimating the arithmetic average deviation of the roughness profile of tested pipelineR aEitimate 、 The estimated value of the maximum height of roughness profileR zEitimate
The step of choosing test and appraisal sample length and test and appraisal evaluation length:For basisR aEitimate R zEitimate Or roughness profile Arithmetic average deviation test and appraisal valueR ae With the test and appraisal value of the maximum height of roughness profileR ze Choose test and appraisal sample lengthlr e Or Test and appraisal measurement lengthln e
The step of test and appraisal detection:For according to test and appraisal sample lengthlr e Tested pipeline inwall is tested and assessed along tested pipeline axis Detection, obtains one groupR ae R ze Value;
Assessment sample length whether the step in proper range:For to the groupR ae WithR ze Numerical value be estimated, its numerical value Whether in test and appraisal sample lengthlr e Or test and appraisal measurement lengthln e The maximum height value of corresponding roughness profileR a And roughness The arithmetic average deviation value of profileR z In the range of value;The group is determined if "Yes"R ae WithR ze Numerical value be pre-selectionR a WithR z Value, and enter next step, returned to if "No" " the step of choosing test and appraisal sample length and test and appraisal evaluation length ";
The step of selecting shorter sample length:For choosing shorter test and appraisal sample lengthlr s , and one group of roughness profile of acquisition The shorter test and appraisal value of arithmetic average deviationR as With the shorter test and appraisal value of the maximum height of roughness profileR zs , described shorter survey Comment sample lengthlr s Refer to:Than in the most short survey selected by " the step of choosing test and appraisal sample length and test and appraisal evaluation length " Sample length is commented also to want short test and appraisal length;
The step of determining sample length:IfR as R zs Numerical value meet corresponding sampling scope, it is determined thatR as R zs Corresponding Shorter test and appraisal sample lengthlr s Sample length used in general measurelr, otherwise select " whether assessing sample length Step in proper range " pre-selectionR a R z The corresponding test and appraisal sample length of valuelr e Sampling used in general measure Lengthlr
The step of general measure:For using sample lengthlrCollection tested pipeline inwall multiple positionR ai WithR zi Value, wherein:i=1,2 ... ...,n,nIt is positive integer;
The step of calculating average value:For calculatingR ai WithR zi Value average value:
,
The step of calculating hydraulic simulation experiment:For by measureAskValue is substituted into Ke Shi formula:
,
The hydraulic simulation experiment λ of tested pipeline is calculated, and draws corresponding hydraulic calculation chart;
In formula:λ- hydraulic simulation experiment;D- internal diameter of the pipeline;k- equivalent roughness, usesInstead of;Re- Reynolds number, calculates Formula isνThe kinematic coefficient of viscosity of-water;V- pipeline mean flow rate;
DescribedR aEitimate R zEitimate R an R zn R ae R ze R as R zs R ai R zi It is the different phase in detection processR a WithR z , changed to show difference with footnote;Describedlr e ln e lr s It is the different phase in detection processlrln, use pin Mark changes to show difference.
2. method according to claim 1, it is characterised in that described sample lengthlrOr measurement lengthlnCorrespondingR a WithR z Scope it is as follows:
R a ≤ 10 microns,lrIt is 2.5 millimeters,lnIt is 12.5 millimeters;10 microns of <R a ≤ 80 microns,lrIt is 8 millimeters,lnIt is 40 millis Rice;
R z ≤ 50 microns,lrIt is 2.5 millimeters,lnIt is 12.5 millimeters;50 microns of <R z ≤ 200 microns,lrIt is 8 millimeters,lnIt is 40 millis Rice.
3. method according to claim 1 and 2, it is characterised in that set after the step of described calculating average value: Calculate standard variance the step of and judge pipeline quality step:
The step of calculating standard variance:
,
The step of judging pipeline quality:For being calculated by formulaR a WithR z Coefficient of dispersionCV Ra CV Rz
,
And judgeCV Ra CV Rz Whether coefficient of dispersion is less thanCVThreshold value, illustrated if "Yes" measured pipeline be it is uniform, it is no Then illustrate the non-material liner of the same race of measured pipe fitting, or the partial tube liner quality of production is not up to standard.
4. method according to claim 3, it is characterised in that described coefficient of dispersionCVThreshold value be:0.3.
CN201710072656.8A 2017-02-10 2017-02-10 A kind of fast appraisement method of pipeline hydraulic frictional resistance Expired - Fee Related CN106777830B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710072656.8A CN106777830B (en) 2017-02-10 2017-02-10 A kind of fast appraisement method of pipeline hydraulic frictional resistance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710072656.8A CN106777830B (en) 2017-02-10 2017-02-10 A kind of fast appraisement method of pipeline hydraulic frictional resistance

Publications (2)

Publication Number Publication Date
CN106777830A true CN106777830A (en) 2017-05-31
CN106777830B CN106777830B (en) 2017-12-19

Family

ID=58956224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710072656.8A Expired - Fee Related CN106777830B (en) 2017-02-10 2017-02-10 A kind of fast appraisement method of pipeline hydraulic frictional resistance

Country Status (1)

Country Link
CN (1) CN106777830B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108005746A (en) * 2017-12-29 2018-05-08 潍柴动力股份有限公司 A kind of engine main oil gallery pressure limiting valve and engine
CN108414113A (en) * 2018-03-15 2018-08-17 山东微感光电子有限公司 The fire alarm system and method for fiber optic temperature are predicted with multi-point temperature coefficient of dispersion
CN109779607A (en) * 2018-12-29 2019-05-21 长江大学 A kind of Frictional Drag of Casing Running coefficient based on hole diameter degree of irregularity determines method and system
CN110929458A (en) * 2019-11-15 2020-03-27 长江大学 Method for calculating annulus comprehensive hydraulic equivalent diameter of irregular borehole section
CN111141250A (en) * 2019-12-31 2020-05-12 中国航发动力股份有限公司 Method for detecting roughness of oil injection hole
CN111191402A (en) * 2020-02-24 2020-05-22 西安石油大学 Method for establishing drag reducer drag reduction effect evaluation and prediction model
CN113138077A (en) * 2021-03-19 2021-07-20 上海船舶电子设备研究所(中国船舶重工集团公司第七二六研究所) Equivalent length testing device and method and gas fire extinguishing system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101050958A (en) * 2007-04-30 2007-10-10 浙江建设职业技术学院 Measuring method for large size structure surface contour curve roughness coefficient
US20110315249A1 (en) * 2009-01-07 2011-12-29 Framo Engineering As Device for providing a controllable pressure reduction
CN203400635U (en) * 2013-07-15 2014-01-22 苏州东大仁智能科技有限公司 Concentrated water discharge device for reverse osmosis seawater desalination system
US20140081606A1 (en) * 2012-09-17 2014-03-20 Mohammed Omer System and method for modeling corrosion-based multiphase flow friction in pipes

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101050958A (en) * 2007-04-30 2007-10-10 浙江建设职业技术学院 Measuring method for large size structure surface contour curve roughness coefficient
US20110315249A1 (en) * 2009-01-07 2011-12-29 Framo Engineering As Device for providing a controllable pressure reduction
US20140081606A1 (en) * 2012-09-17 2014-03-20 Mohammed Omer System and method for modeling corrosion-based multiphase flow friction in pipes
CN203400635U (en) * 2013-07-15 2014-01-22 苏州东大仁智能科技有限公司 Concentrated water discharge device for reverse osmosis seawater desalination system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张福坤等: "输气管道水力摩阻因数计算", 《油气储运》 *
苑伟民: "显式Colebrook-White摩阻系数方程", 《天然气与石油》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108005746A (en) * 2017-12-29 2018-05-08 潍柴动力股份有限公司 A kind of engine main oil gallery pressure limiting valve and engine
CN108005746B (en) * 2017-12-29 2024-04-16 潍柴动力股份有限公司 Engine main oil duct pressure limiting valve and engine
CN108414113A (en) * 2018-03-15 2018-08-17 山东微感光电子有限公司 The fire alarm system and method for fiber optic temperature are predicted with multi-point temperature coefficient of dispersion
CN109779607A (en) * 2018-12-29 2019-05-21 长江大学 A kind of Frictional Drag of Casing Running coefficient based on hole diameter degree of irregularity determines method and system
CN110929458A (en) * 2019-11-15 2020-03-27 长江大学 Method for calculating annulus comprehensive hydraulic equivalent diameter of irregular borehole section
CN110929458B (en) * 2019-11-15 2023-08-18 长江大学 Annulus comprehensive hydraulic equivalent diameter calculation method for irregular well bore section
CN111141250A (en) * 2019-12-31 2020-05-12 中国航发动力股份有限公司 Method for detecting roughness of oil injection hole
CN111141250B (en) * 2019-12-31 2021-06-15 中国航发动力股份有限公司 Method for detecting roughness of oil injection hole
CN111191402A (en) * 2020-02-24 2020-05-22 西安石油大学 Method for establishing drag reducer drag reduction effect evaluation and prediction model
CN113138077A (en) * 2021-03-19 2021-07-20 上海船舶电子设备研究所(中国船舶重工集团公司第七二六研究所) Equivalent length testing device and method and gas fire extinguishing system

Also Published As

Publication number Publication date
CN106777830B (en) 2017-12-19

Similar Documents

Publication Publication Date Title
CN106777830B (en) A kind of fast appraisement method of pipeline hydraulic frictional resistance
CN108254032A (en) River ultrasonic wave time difference method method of calculating flux
Jonsson et al. Experiments on turbulent-flow phenomena in eccentric annular ducts
CN206656765U (en) A kind of combination type multifunction water leveling ruler
CN108151766A (en) Localization method, magnetic nail positioning navigation error correcting method and the positioning device of magnetic nail
CN106092137B (en) The outdoor calibrator (-ter) unit and method of a kind of vehicle-mounted three-dimensional laser pavement detection system
De Chiffre et al. Investigation on capability of the reaming process using minimal quantity lubrication
CN103643997A (en) Method for computing gas emission quantity of low-wind-speed coal tunneling face
CN106225654B (en) A kind of concrete surface roughness measuring device and its measurement method
CN110045144A (en) A kind of closed conduct flow rate of liquid distribution method for automatic measurement and device
CN101187538A (en) Screw gauge detection method
CN105115455B (en) A kind of alignment schemes of the different gauge head measurement data of triaxial coordinate measuring apparatus
CN113358048A (en) Method for recognizing beam structure deflection by graph multiplication based on distributed optical fiber
CN102538876A (en) Tester for parameter of flow field of gas static pressure restrictor
CN117518982B (en) Method and system for improving machining precision of machine tool
CN108507444A (en) The detection device and method of workpiece specific position Inner bore dias
CN204495229U (en) A kind of guide hole cubing
CN204611352U (en) A kind of Portable thermal hydraulic piping detection device
CN112113481A (en) Device and method for measuring dam body breach geometric shape in water tank test
CN107220513A (en) A kind of method and its application for calculating mud-rock flow mean flow rate
CN108871152A (en) A kind of internal diameter of the pipeline measuring device and method
CN106643443A (en) Coaxiality detection device and method for mixing cylinders of concrete mixer
CN205300473U (en) Mar measuring apparatu
CN101498672B (en) Method for detecting surface finish quality of dust-free cloth
CN105004385B (en) A kind of colliery lotion pipeline flow-measuring method

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20171219

Termination date: 20190210

CF01 Termination of patent right due to non-payment of annual fee