Scientific journal
Advances in current natural sciences
ISSN 1681-7494
"Перечень" ВАК
ИФ РИНЦ = 0,775

Features of measurements and dynamics modeling expense of a small waterway of a spring «A crystal spruce»

Mazurkin P.M., 1 Tarasova E.I. 1
1 Volga State University of Technology
In the summer of 2012 monitoring of a consumption of water on a small waterway was carried out. The measured vessel was accepted in the form of a ladle in capacity one liter. All measurements were carried out in the evening from 17-00 p.m. Therefore the current time undertakes the whole days. The model of dynamics has two components: the first component is the law of exponential growth, and the second wave indignation with variables amplitude and frequency of fluctuation. The modeling technique with procedures is shown: 1) identification of the permanent member; 2) on the remains from the permanent member, consistently complicating a design, wave function is identified; 3) the permanent member is combined with wave function; 4) the trend design to steady not the wave law becomes complicated.
spring
expense
measurement
law

Динамика – это изменение во времени [3–5]. Для родника динамическими показателями являются время наполнения емкости и расход водотока, вычисляемый как отношение объема мерного сосуда к времени его наполнения.

Эксперименты проводили в местности, которая богата родниками, в так называемом «родниковом крае». Действительно, Куженерский район можно назвать настоящим родниковым краем [2]. Проходящий по территории района Вятский увал создает уникальное движение подземных вод на своих склонах, из-за чего в этих местах много источников, ключей ручьев. На территории поселка Куженер находится немало мелких родников.

Одним из них является «Хрустальная ель», родник находится в долине реки Немда, под Якайсуртом в овраге. Его кристально-чистая вода по душе многим жителям.

Ключ бьёт из-под корней ели и образует несколько ручейков. Сам родник представляет собой небольшой ключ, который под напором выходит из земли в специальном углублении. Вместе с подземными выходами образуется некая система водотоков в роднике, которая собирается в деревянном лотке и вытекает в нескольких крупных струях.

Устойчивые законы. По принципу «от простого к сложному» в таблице 1 предложены «кирпичики» для построения, по ходу структурно-параметрической идентификации биотехнического закона [3, 5], искомой закономерности.

При моделировании временных рядов тренд нужно вначале искать по закону экспоненциального роста или гибели (спада). Тренды во времени всегда имеют предысторию, так как ничто не возникает в историческом поведении на пустом месте, то есть с нуля. Тенденции по структурной динамике вполне могут начаться с нулевого элемента системы, т.е. в новых условиях. Все шесть устойчивых законов распределения из таблицы 1 являются частными случаями биотехнического закона, показанного в нижней строке.

Уровни коэффициента корреляции. Теперь можно определиться с уровнями коэффициента корреляции. Это позволит затем ранжировать полученные по функциональным факторным связям по тесноте коррелятивной вариации (термин по Ч. Дарвину).

Таблица 1

Математические конструкты для построения статистической модели [5]

Фрагменты без предыстории

изучаемого явления или процесса

Фрагменты с предысторией

изучаемого явления или процесса

y = ax – закон линейного роста или спада (при отрицательном знаке перед правой стороной приведенной формулы)

y = a – закон не влияния принятой переменной на показатель, который имеет предысторию значений

y = axb – закон показательного роста (закон показательной гибели y = ax–b не является устойчивым, из-за бесконечности при нулевом значении объясняющей переменной

y = a exp(±cx) – закон Лапласа (Ципфа в биологии, Парето в экономике, Мандельброта в физике) экспоненциального роста или гибели, относительно которого создан метод операторных исчислений

y = axb(–cx) – биотехнический закон в упрощенной форме

y = a exp(±cxd) – закон экспоненциального роста или гибели, – по проф. П.М. Мазуркину

y = axb(–cxd) – биотехнический закон – проф. П.М. Мазуркин

Всего рассмотрим восемь уровней адекватности полученной после идентификации закономерности, начиная от коэффициента корреляции, равного 0,000. При этом для научно-технического творчества особо рекомендуются к вниманию зависимости с корреляцией от 0,300 до 0,500. Их считают слабыми связями и, не задумываясь, отбрасывают.

Но, в других условиях функционирования, именно слабые по результатам текущих наблюдений и прошлых экспериментов в будущем могут стать сильными.

Таблица 2

Коэффициент корреляции между факторами

Интервал коэффициента корреляции

Характер тесноты связи

Существующая классификация

Предлагаемая классификация

1,0000

сильная связь

однозначная

0,9000…1,0000

сильнейшая

0,7000…0,9000

сильная

0,500…0,700

слабая связь

средняя

0,300...0,500

слабоватая

0,100…0,300

нет связи

слабая

0,000…0,100

слабейшая

0,000

нет связи

Летом 2012 года нами проведен мониторинг расхода воды на малом водотоке в роднике «Хрустальная Ель». Все данные по результатам эксперимента представлены в табл. 3. Методика в основном та же, что и для предыдущего эксперимента. Отличие заключается в том, что мерный сосуд был ковшом в один литр.

Все измерения проводились вечером с 17-00 часов. Поэтому текущее время τ берется целыми сутками. Расход малого водотока q взяли в размерности мл/c, т.к. в размерности л/с значения объемного расхода были значимыми для крупного водотока и были слишком малыми по числам для мелкого по расходу водотока родника.

Для моделирования в программной среде CurveExpert соотношения чисел на осях абсцисс и ординат должны быть соразмерными. Продолжительность в 76 суток достаточна, чтобы выявлять к тренду волновые закономерности возмущения изучаемого признака.

Ранее измерения большого водотока получились неточными. Большой водоток требует быструю реакцию или автоматическую регистрацию. Не соблюдая такие условия, может возникнуть большая погрешность при измерении крупного водотока.

На рис. 1 показано фото места истечения одного малого водотока.

Расход в нем слишком мал, поэтому за мерное устройство приняли ковш объемом в один литр. И его было достаточно, чтобы мерный сосуд наполнялся не менее чем за 15 секунд (расчеты и обоснование погрешности измерений приведены в седьмой главе). Если для наполнения мерного сосуда времени требуется больше, то значит, при ручной фиксации времени секундомером результаты измерений будут также точнее.

Таблица 3

Результаты измерений расхода воды малого водотока из родника «Хрустальная ель»

№ п/п

Дата

Время измерения τ, сутки

Время наполнения T, с

Расход q, л/с

Расход

q, мл/с

№ п/п

Дата

Время измерения τ, сутки

Время наполнения T, с

Расход q, л/с

Расход q, мл/с

1

2

3

4

5

6

7

8

9

10

11

12

1

23.07.12

0

57.38

0,0174

17,4

66

 

12

70.43

0,0142

14,2

2

17-00

0

56.68

0,0176

17,6

67

 

12

69.75

0,0143

14,3

3

 

0

54.32

0,0184

18,4

68

 

12

68.11

0,0147

14,7

4

 

0

63.26

0,0158

15,8

69

 

12

74.40

0,0134

13,4

5

 

0

59.28

0,0169

16,9

70

 

12

72.45

0,0138

13,8

6

 

0

52.08

0,0192

19,2

71

6.08.12

14

70.59

0,0142

14,2

7

 

0

65.92

0,0152

15,2

72

 

14

68.79

0,0145

14,5

8

 

0

63.18

0,0158

15,8

73

 

14

72.10

0,0139

13,9

9

 

0

69.12

0,0145

14,5

74

 

14

76.38

0,0131

13,1

10

 

0

65.27

0,0153

15,3

75

 

14

72.67

0,0138

13,8

11

24.07.12

1

61.99

0,0161

16,1

76

 

14

70.12

0,0143

14,3

12

17-00

1

65.52

0,0153

15,3

77

 

14

74.93

0,0133

13,3

13

 

1

72.20

0,0139

13,9

78

 

14

75.35

0,0133

13,3

14

 

1

70.84

0,0141

14,1

79

 

14

70.43

0,0142

14,2

15

 

1

69.97

0,0143

14,3

80

 

14

72.02

0,0139

13,9

16

 

1

70.71

0,0141

14,1

81

15.08.12

23

77.34

0,0129

12,9

17

 

1

74.03

0,0135

13,5

82

 

23

86.22

0,0116

11,6

18

 

1

72.67

0,0138

13,8

83

 

23

83.31

0,0120

12

19

 

1

75.86

0,0132

13,2

84

 

23

85.79

0,0117

11,7

20

 

1

75.63

0,0132

13,2

85

 

23

81.20

0,0123

12,3

21

 

1

72.19

0,0139

13,9

86

 

23

86.96

0,0115

11,5

22

25.07.12

2

73.78

0,0136

13,6

87

 

23

80.43

0,0124

12,4

23

 

2

73.01

0,0137

13,7

88

 

23

78.40

0,0128

12,8

24

 

2

75.87

0,0132

13,2

89

 

23

82.25

0,0122

12,2

25

 

2

73.90

0,0135

13,5

90

 

23

89.20

0,0112

11,2

26

 

2

75.97

0,0132

13,2

91

13.09.12

52

105.00

0,0095

9,5

27

 

2

77.32

0,0129

12,9

92

 

52

105.73

0,0095

9,5

28

 

2

75.27

0,0133

13,3

93

 

52

101.49

0,0099

9,9

29

 

2

71.12

0,0141

14,1

94

 

52

103.47

0,0097

9,7

30

 

2

71.88

0,0139

13,9

95

 

52

100.87

0,0099

9,9

31

 

2

66.06

0,0151

15,1

96

 

52

106.29

0,0094

9,4

32

26.07.12

3

76.65

0,0130

13

97

 

52

104.50

0,0096

9,6

33

 

3

77.35

0,0129

12,9

98

 

52

104.41

0,0096

9,6

34

 

3

70.65

0,0142

14,2

99

 

52

102.51

0,0098

9,8

35

 

3

69.72

0,0143

14,3

100

 

52

103.22

0,0097

9,7

36

 

3

65.01

0,0154

15,4

101

14.09.12

53

104.59

0,0096

9,6

37

 

3

63.69

0,0157

15,7

102

 

53

102.02

0,0098

9,8

38

 

3

75.69

0,0132

13,2

103

 

53

100.43

0,0100

10

39

 

3

71.15

0,0141

14,1

104

 

53

105.35

0,0095

9,5

40

 

3

73.41

0,0136

13,6

105

 

53

100.93

0,0099

9,9

41

 

3

71.41

0,0140

14

106

 

53

102.43

0,0098

9,8

42

27.07.12

4

74.11

0,0135

13,5

107

 

53

103.25

0,0097

9,7

43

 

4

73.25

0,0137

13,7

108

 

53

103.78

0,0096

9,6

44

 

4

71.20

0,0140

14

109

 

53

105.02

0,0095

9,5

45

 

4

70.93

0,0141

14,1

110

 

53

104.12

0,0096

9,6

46

 

4

74.98

0,0133

13,3

111

6.10.12

75

106.35

0,0094

9,4

47

 

4

70.43

0,0142

14,2

112

 

75

108.40

0,0092

9,2

48

 

4

69.75

0,0143

14,3

113

 

75

108.79

0,0092

9,2

49

 

4

68.11

0,0147

14,7

114

 

75

109.10

0,0092

9,2

50

 

4

74.40

0,0134

13,4

115

 

75

107.64

0,0093

9,3

51

28.07.12

5

72.45

0,0138

13,8

116

 

75

108.21

0,0092

9,2

52

 

5

70.59

0,0142

14,2

117

 

75

108.07

0,0093

9,3

53

 

5

68.79

0,0145

14,5

118

 

75

109.15

0,0092

9,2

54

 

5

72.10

0,0139

13,9

119

 

75

107.19

0,0093

9,3

55

 

5

76.38

0,0131

13,1

120

 

75

107.29

0,0093

9,3

56

 

5

72.67

0,0138

13,8

121

7.10.12

76

108.49

0,0092

9,2

57

 

5

70.12

0,0143

14,3

122

 

76

108.25

0,0092

9,2

58

 

5

74.93

0,0133

13,3

123

 

76

108.48

0,0092

9,2

59

 

5

75.35

0,0133

13,3

124

 

76

109.03

0,0092

9,2

60

 

5

70.48

0,0142

14,2

125

 

76

108.72

0,0092

9,2

61

4.08.12

12

74.11

0,0135

13,5

126

 

76

108.81

0,0092

9,2

62

 

12

73.25

0,0137

13,7

127

 

76

107.97

0,0093

9,3

63

 

12

71.20

0,0140

14

128

 

76

107.95

0,0093

9,3

64

 

12

70.93

0,0141

14,1

129

 

76

108.81

0,0092

9,2

65

 

12

74.98

0,0133

13,3

130

 

76

109.12

0,0092

9,2

pic_42.tif

Рис. 1. Малый водоток родниковой воды

В табл. 4 представлены значения показателей вариации, рассчитанные по формулам описательной статистики. Данные обрабатывались в программной среде Ecxel.

Показатель расхода воды имеет относительно высокий коэффициент вариации, но всё же не соблюдается общеизвестное условие Eqn68.wmf.

Это означает, что есть достаточно высокая погрешность усреднения статистической выборки. При этом среднее арифметическое значение определено надежно и этому свидетельствует показатель точности P = 3,92 %.

Таблица 4

Статистические параметры расхода воды эксперимента

Статистические показатели

Формула

Результат расчета

Количество наблюдений, шт.

n

130

Размах R, мл/с

R = xmax – xmin

10

Минимум xmin, мл/с

Выбор из ряда данных

9,2

Максимум xmax , мл/с

Выбор из ряда данных

19,2

Среднее арифметическое Eqn59.wmf, мл/с

Eqn62.wmf

12,604

Ср. линейное отклонение Eqn60.wmf, мл/с

Eqn63.wmf

0,000154

Ср. квадратичное отклонение Eqn61.wmf

Eqn64.wmf

5,631

Коэффициент вариации Vσ, %

Eqn65.wmf

44,68

Ср. ошибка выборки mσ

Eqn66.wmf

0,49

Показатель точности P, %

Eqn67.wmf

3,92

В связи с этим следует выполнить статистическое моделирование, а по измеренным параметрам малого водотока можно выявлять биотехнические закономерности с коэффициентом корреляции выше 0,9.

По точкам в окне программной среды сразу же видно, что изменение времени наполнения T в зависимости от текущего времени τ с начала эксперимента происходит волнообразно. В таком случае нужно определить тренд (тенденцию). Эксперимент в 76 суток оказался недостаточным по продолжительности для моделирования экспоненциального закона.

В связи с этим остается идентифицировать наипростейшую формулу закона не влияния, то есть конструкцию вида Eqn69.wmf среднего арифметического значения (рис. 2 и 3).

pic_43.wmf

Рис. 2. График среднего арифметического значения времени наполнения

pic_44.wmf

Рис. 3. Остатки после модели (1)

Вначале была получена простейшая формула

T = 82,44115. (1)

Далее была идентифицирована простая формула косинуса с постоянными параметрами модели. Запись волнового возмущения приобретает вид формулы

Eqn71.wmf (2)

Отрицательный знак показывает кризисный характер колебательного возмущения. В знаменателе находится число 66,04238 половины периода колебательного значения. Весь период равен примерно 132 сутки. Частота колебания будет равна 132 или 1/132.

В этой модели половина амплитуды колебательного процесса равна 16,46117, а число π принимается символом pi в программе CurveExpert-1.3 с 18 знаками после запятой. Полупериод равен примерно 66 суток, то есть период колебания ожидается 132 сутки. Однако остатки после формулы (2) позволяют усложнить конструкцию статистической модели (2). Дополним постоянную амплитуду экспоненциальным законом, то есть введем функцию экспоненты.

Тогда можем дальше записать формулу так:

Eqn72.wmf (3)

После усложнения конструкции модели амплитуда изменяется по устойчивому закону экспоненциального роста, аналогично сложную зависимость может иметь и период колебательного изменения. Но начало волны по формуле (3) не совпадает с началом координат при условии τ = 0. Поэтому сдвиг начала волны колебательного возмущения учитывается дополнительным параметром модели, помещаемым внутри функции косинуса, и получимформулу:

Eqn73.wmf (4)

В этой формуле число 0,65128 показывает сдвиг начала волны относительно начала координат. Этот сдвиг измеряется в радианах.

Остатки после формулы (4) показали, что полупериод колебания является непостоянным. Поэтому усложняем конструкцию волновой функции. Повторная идентификация дает формулу

Eqn74.wmf (5)

Формулу (5) можно записать в другом виде:

Eqn75.wmf (5а)

Eqn76.wmf

Eqn78.wmf

где A – амплитуда (половина) колебания, с; p – период (половина) колебания, сутки.

Полупериод колебания в начале эксперимента, то есть 23.07.12, был равен 111,8 суток, но затем он уменьшается. Это происходит потому, что летом характерна стабильная погода, а к осени влажность и температура атмосферного воздуха сильно и часто колеблется (рис. 4).

pic_45.wmf

Рис. 4. Динамика времени наполнения

После следующего сеанса идентификации с усложнением конструкции программная среда дает аварийный останов. Поэтому окончательная модель приобретает вид формулы (5).

Теперь нужно соединить тренд (1) с волновой функцией (5) и получить общую закономерность:

T = T1 + T2; (6) (6)

T1 = 76,56800;

Eqn79.wmf

Eqn80.wmf

Eqn81.wmf

Далее усложним первую составляющую тем, что постоянный член дополним экспонентой. После идентификации получена формула

T = T1 + T2; (7)

Eqn82.wmf

Eqn83.wmf

Eqn84.wmf

Eqn85.wmf

Первая составляющая стала законом экспоненциального роста. Вторая составляющая осталась кризисной, но по сравнению с формулой (5) изменилась конструкция полупериода: она стала возрастающей. Частота колебания уменьшается, т.е. время наполнения успокаивается.

Далее формулу (7) доведем до полноты конструкции (рис. 5, 6) в виде

Eqn86.wmf (8)

pic_46.wmf

Рис. 5. Динамика времени наполнения ковшика в 1 л малым водотоком

pic_47.wmf

Рис. 6. Остатки от модели (8)

Готовая модель образует две составляющие, которые четко характеризуют волновое возмущение наполнения ковша водой родником. Первая составляющая являются известным законом экспоненциального роста при интенсивности роста 0,27262, а вторая – волновое возмущение с переменными амплитудой и частотой колебания.

По приведенной методике моделирования была приведена идентификация закономерностей по расчетному расходу малого водотока.

Кратко эта методика содержит следующие процедуры:

1) выявление постоянного члена, то есть тренда в виде среднего арифметического значения;

2) по остаткам от постоянного члена, усложняя конструкцию модели, идентифицируется волновая функция;

3) постоянный член совмещается с волновой функцией;

4) усложняется конструкция тренда (постоянного члена) до полной первой составляющей, т.е. до устойчивого не волнового закона.

На основе этой методики была получена (рис. 7 и 8) модель

Eqn87.wmf (9)

pic_48.wmf

Рис. 7. Динамика расхода водотока

pic_49.wmf

Рис. 8. Остатки от модели (9)

Оказалось, что первый член формулы (9) является постоянным. Это также объясняется относительно недостаточной продолжительностью эксперимента.

Коэффициент корреляции модели (9) равен 0,9241, что указывает на сильнейшую факторную связь. Чтобы объективно доказать волновой характер изменения расхода воды, необходимо добиться ламинарного потока жидкости. Крупный водоток сильно пульсирует, а малый водоток, как правило, является ламинарным.

Кроме того, время наполнения ковшика постоянной емкости, например, в один литр, значительное и это позволяет уменьшить погрешности измерения секундомером. Поэтому на одном малом водотоке нужно продолжить измерения в течение нескольких месяцев.

Для измерения расхода воды от водотоков разной мощности у одного и того же родника для ориентировочного анализа водотоков применяют мерные сосуды такой емкости, чтобы они заполнялись не быстрее пять секунд. Тогда, при реакции человека по 0,3 с на включение и отключение секундомера, погрешность измерений времени наполнения составит 0,6 с. Измерения будут проведены с общей погрешностью измерений не более 100∙0,6/5 = 12 %.

С учетом времени взятия пробы, сезонности и особенностей рельефа [1] можно будет составлять исходные данные для изучения влияния орографических факторов местности и количества осадков в виде дождя и снега на удельный расход воды родником. Это позволит выявлять биотехнические закономерности влияния лесов и растительного покрова в целом, а также ведения сельскохозяйственных и иных работ вокруг родников, на окружающую среду и учет антропогенного влияния на расход воды в сети родников.