Да. ЭДС всегда фазная.
Чтобы не перепутать, правило: Z1=Z2=ZL-Zm, Z0=ZL+2*Zm (Zm - взаимоиндукция между фазами, ZL - собственное сопротивление фазы).
Тогда после всех подстановок и сокращений Ia=3*E1/(Z1+Z2+Z0) получим Ia=Ea/ZL - фазная ЭДС делить на сопротивление фазы (думаю, логичная формула :) ).

Вот из этой таблицы данные в эксель подставил, поделил одно на другое, подставил в формулы расчета сопротивлений и получил цифры в районе 40%.

Смысл в том, что параметры трансформатора (сквозное сопротивление, полное сопротивление) и стандартные формулы расчета ориентированы на внешние КЗ, при внутренних придется самостоятельно формулу выводить, и там (при точном расчете) придется использовать не сквозные или полные сопротивления, а какие-то параметры, которые из них косвенным образом выводятся. В тех формулах, которые привели вы (СТО Газппрома), видимо, берутся такие большие сопротивления резисторов, что ток получается преимущественно активным и определяется в основном им. Остается только убедиться, что сопротивление резистора действительно много больше всего остального.

Ветвь намагничивания бывает по прямой последовательности (это то, что можно получить из паспортных данных Ixx), а бывает по нулевой последовательности.
Для получения параметров ветви намагничивания по нулевой последовательности нужно делать специальный опыт холостого хода нулевой последовательности, и мерять Ixx0, там схема опыта не такая, как для просто Ixx. И нужно еще, по идее, делать опыт короткого замыкания нулевой последовательности, и мерять uk0% (там тоже своя схема опыта, не такая, как для uk%).
Опытов по нулевой последовательности производители тупо не делают, наверное, им лень или типа того.

В итоге для моделирования недостаточно данных. В симулинке, если взять модель трансформатора, есть специальный параметр zero-sequence excitation current when delta winding opened, это и есть тот самый Ixx0. Когда я в первый раз увидел этот параметр, я вообще был хз, чему его задавать и что это, потом только, когда разобрал модель трансформатора в матричной форме ( в виде магнитосвязанных катушек), то тогда стало понятно, что двумя опытами, приведенными в паспорте все параметры этих катушек не задать, будет неопределенность. В хелпе матлаба рекомендуют задавать Ixx0=100%, ссылка вот https://www.mathworks.com/help/physmod/ … hworks.com
Тогда получается, что Z0=100*Uном.ф^2/(Ixx0%*Sном.ф)=Uном.ф^2/Sном.ф=Uном^2/Sном, т.е. Z0 равно сопротивлению номинальной нагрузки трансформатора.
В книге Беляева, если обработать данные из таблицы для каких-то старых трансформаторов, получается Ixx0≈250%...

Да, именно это я и имел ввиду; видимо я не понял этого из ваших предшествующих сообщений. Что касается производителей тр-ров, то действующий ГОСТ не предполагает измерения этого сопротивление на заводе. Если правильно помню (опять же согласно ГОСТ на тр-ры) заказчик при заказе может попросить завод предоставить эти данные, и заводу ничего не останется, как измерить самому и предоставить. Возможно у них есть некоторая архивная база подобных измерений, нужно бы как-нибудь написать.
В матлабе дефолтная величина тока Ixx0% 100% мне кажется отвлечена от реальности и завышена. Собственно там важна и величина активных потерь х.х. по нулевой последовательности.
alexej p., что касается постоянной равной 1, то это как в школе: примем всю работу за единицу, тут же принимаем всю обмотку за единицу. Тогда сопротивление замкнувшейся части обмотки пропорционально числу витков от начала обмотки до места замыкания.

А реально ток ОЗЗ будет изменяться в зависимости от сопротивления до замкнувшейся части обмотки ВН трансформатора КТП соединенной в треугольник.
А в формуле СТО считается, что при ОЗЗ в середине обмотки статора ток обязательно снизится в половину.
http://rzia.ru/uploads/8706/thumbnail/B1EZr_vYuyai29tHSGw0.png

Да появится возможность измерим, очень интересно узнать и результатами конечно поделюсь.
Спасибо, всем за помощь коллеги.

Обычно используются (в почти 100% случаев) только две ступени с независимой выдержкой.
Одна (защита в режиме ОЗЗ, при отказе собственной защиты или выключателя):
Iсраб = Irn/Кн
Тсраб = Тдоп/Кн
Где Irn - номинальный ток резистора.
Тдоп - допустимое время пребивании резистора под ном. током (обычно 10 сек).
Кн = 1,3 - 1,5 (и для тока и времени)
Другая (защита от недопустимой длительной несимметрии):
Iсраб = (Кн/Кв).Ir.прод (продолжительно допустимый ток через резистор, нормально 10% номинала).
Тсраб сигн = 3,0 сек
Тсраб откл > 30-60 сек
Все это конечно недостаточно и осознается, но использование самых простых защит свело защита резистора до этого.
После того как сгорело несколько резисторов именно от длительной несимметрии стали думать, но есть ли решение не знаю...
У меня своя выдумка. В одном проекте (СН 6,3 кВ ТЭЦ на сборных шинах 2х63 МВт) кроме выше перечисленных добавил МТЗ с зависимой х-кой (типa LIT). Если хотите могу выслать расчет. Но предупреждаю - это я так делаю, а не официально принятое решение.
Тогда не пользовался фукцией 49 RMS (ее в реле просто не было).

Для защит присоединений используем трехтрансформаторный фильтр НП (от бубликов постепенно оказываемся). Для защиты резистора - один обычный ТТ в заземляющем проводе (встраивается после резистора, в смысле м/у резистором и землей).

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